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1  Introduction 


For  the  past  several  decades,  the  pavement  construction  industry  has 
moved  away  from  material  and  method  specifications,  toward  quality  assur¬ 
ance  specifications.  Material  and  method  specifications  required  complete 
direction  of  construction  by  a  representative  of  the  contractual  “owner.  ”  The 
contractor  was  directed  to  use  specific  materials  and  specific  methods  of 
placement.  This  type  of  specification  created  several  problems.  It  prevented 
the  contractor  from  using  innovative  materials  or  techniques  and  it  obligated 

the  owner  to  accept  the  final  product  regardless  of  quality. 

Quality  assurance  specifications  depend  primarily  on  end-result  criteria 
which  require  the  contractor  to  take  responsibility  for  supplying  a  product  or 
an  item  of  construction.  The  owner’s  responsibility  is  then  either  to  accept 
or  reject  the  final  product  or  to  apply  a  price  adjustment  that  compensates  for 
the  degree  of  compliance  with  the  specifications  (TRB  1996). 

The  end-result  criteria  used  for  a  quality  assurance  specification  must  be 
related  to  pavement  performance.  These  criteria  involve  testing  for  material 
and  structural  characteristics  that  are  known  to  correlate  with  pavement  per¬ 
formance  (e.g.  void  content  and  thickness  for  asphalt  concrete).  The  charac¬ 
teristics  used  for  end-result  criteria  are  often  not  fundamental  engineering 
properties  (e.g.  resilient  modulus  or  fatigue  resistance)  because  they  have  to 
be  highly  repeatable  and  amenable  to  timely  acceptance  testing  (TRB  1996). 

To  execute  a  quality  assurance  specification  properly,  the  contractor  must 
implement  a  quality  control  (QC)  program  and  the  owner  must  implement  a 
quality  assurance  (QA)  program.  A  QC  program  involves  sampling  and  test¬ 
ing  during  construction  to  control  the  level  of  quality  of  the  final  product.  A 
QA  program  involves  all  the  actions  necessary  to  provide  confidence  that  a 
product  or  facility  will  perform  satisfactorily  in  service.  The  QA  program 
must  include  acceptance  sampling  and  testing  to  determine  if  the  quality  of 
produced  material  or  construction  is  acceptable  in  terms  of  the  specifications 


The  results  of  QA  acceptance  tests  are  compared  to  acceptance  limits, 
which  usually  include  provisions  for  pay  adjustments.  These  acceptance 
limits  must  be  realistic,  particulary  when  nonconformance  results  in  reduced 
pay  for  a  contractor.  The  principal  obstacle  for  developing  realistic  limits  is 
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the  ability  to  predict  with  some  degree  of  accuracy  the  amount  of  variability 
that  should  be  expected  in  “satisfactory”  construction  (Nicotera  1974). 


Problem 


The  pavement  industry  has  not  maintained  or  published  extensive  variabil¬ 
ity  data  related  to  materials  and  construction,  as  have  other  fields  of  engi¬ 
neering  and  science.  Pavement  construction  has  long  been  regarded  as  an  art 
that  is  highly  dependent  on  the  skills  and  experience  of  the  field  engineer,  so 
decisions  have  often  been  based  on  engineering  judgement.  The  use  of  tradi¬ 
tional  material  and  method  specifications  for  pavement  projects  has  been 
accompanied  by  a  minimal  amount  of  conformance  testing.  Consequently, 
there  has  typically  not  been  an  organized  method  for  resolving  disputes  when 
materials  or  methods  used  for  pavement  construction  are  suspected  to  not 
conform  with  specification  requirements.  Other  industries,  in  contrast, 
involve  production  lines  in  a  controlled  environment  that  includes  confor¬ 
mance  testing.  Successful  automobile  manufacturing,  for  example,  has  relied 
on  a  statistical  approach  to  quality  control  for  over  40  years. 

The  pavement  industry  is  now  adopting  statistics-based  specifications  and 
reliability  assessments  for  several  reasons.  First,  the  pace  of  pavement  con¬ 
struction  projects  has  increased,  necessitating  that  the  parties  involved  share 
responsibilities.  Secondly,  the  duties  of  engineers  has  become  more  wide¬ 
spread,  requiring  that  they  spend  time  away  from  their  field  projects.  Fin¬ 
ally,  modem  construction  has  been  besieged  with  litigation.  Owners  must 
therefore  protect  themselves  by  writing  clear  specification  requirements  that 
are  based  on  fair  assumed  risks  by  all  contractual  parties. 


Scope 


This  report  has  two  purposes.  First  of  all,  it  describes  probabilistic  and 
statistical  tools  that  are  useful  for  analyzing  the  variability  of  pavement  mate¬ 
rials  and  pavement  structures.  Secondly,  this  report  presents  a  summary  of 
published  variability  data  pertaining  to  pavements.  These  tools  and  data  will 
be  useful  for  the  development  of  fair  and  enforceable  quality  assurance 
specifications. 

Basic  statistical  concepts  are  reviewed  in  Chapter  2  to  prepare  the  reader 
for  the  variability  analysis  tools  presented  in  Chapter  3.  Published  variability 
data  are  summarized  in  Chapter  4  and  are  presented  in  detail,  with  refer¬ 
ences,  in  the  appendices. 
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2  Statistical  Concepts 


This  chapter  provides  a  review  of  basic  probabilistic  and  statistical  concepts, 
which  will  be  used  for  analyses  of  material  variability  later  in  this  report. 
This  chapter  should  be  useful  to  those  who  develop,  write,  and  implement 
statistics-based  construction  specifications. 


Terminology 


Population  versus  sample 

A  population  is  a  set  of  data  that  includes  an  entire  entity  of  interest.  If 
an  Army  installation  is  conducting  inventory  on  its  highway  system,  the 
design  thicknesses  of  all  portland  cement  concrete  (PCC)  pavements  would 
be  considered  a  population  of  data.  A  sample  is  a  set  of  data  that  is  selected 
randomly  from  within  a  population,  typically  to  represent  the  population 
from  which  it  was  derived.  If  an  Army  installation  were  to  make  a  judge¬ 
ment  concerning  the  thicknesses  of  its  PCC  pavements  by  coring  100 
randomly-selected  slabs,  the  selected  slabs  would  constitute  a  sample  of  the 
installation’s  population  of  PCC  pavements. 


Unit,  observation,  and  variable 

A  set  of  data  is  composed  of  information  gathered  from  a  number  of 
experimental  units.  The  information  gathered  from  each  experimental  unit  is 
called  an  observation.  Each  observation  may  consist  of  one  or  more  pieces 
of  information,  called  variables. 

Example.  If  an  Army  installation  wishes  to  gather  information  about 
nearby  sources  of  asphalt  cement,  each  source  of  asphalt  cement  would  be  a 
unit.  The  observation  for  each  of  these  units  may  consist  of  several  vari¬ 
ables,  such  as  specific  gravity,  solubility,  and  viscosity  at  a  given 
temperature. 
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Types  of  variables 

Variables  may  be  continuous  or  discrete.  Continuous  variables  can  take 
on  any  value  within  a  given  interval. 

Example.  The  asphalt  cement  content  of  hot-mix  asphalt  (HMA),  which 
could  be  measured  from  extractions  performed  on  pavement  core  samples,  is 
a  continuous  variable.  Conceptually,  the  asphalt  cement  content  could  be 
any  value  between  and  including  0  percent  and  100  percent.  Similar  to 
almost  all  continuous  variables,  the  reported  measurements  of  asphalt  cement 
content  are  limited  to  discrete  increments  due  to  the  limited  precisions  of 
measurement  devices.  Reported  asphalt  cement  contents  are  typically  limited 
to  increments  of  one-tenth  of  a  percent. 

A  discrete  variable  must  assume  a  single  category  from  within  a  predeter¬ 
mined  list.  These  categories  may  be  based  on  a  nominal  scale  or  an  ordinal 
scale.  In  a  nominal  scale,  each  category  is  given  a  name.  The  category 
names  may  be  alphabetic,  numeric,  or  alphanumeric. 

Example.  If  an  Army  installation  is  conducting  inventory  on  its  roadway 
system,  they  may  classify  each  pavement  section  by  an  alphabetic  name  that 
reflects  the  nature  of  its  design:  rigid,  flexible,  or  unsurfaced.  As  an  exam¬ 
ple  of  numeric  categories  in  a  nominal  scale,  an  installation  may  classify  its 
roadways  according  to  their  number  of  lanes.  All  observations  for  the  num¬ 
ber  of  lanes  would  be  discrete;  commonly,  2  or  4.  In  this  case,  the  discrete 
variable  names  would  have  quantitative  implications. 

An  ordinal  scale  for  a  discrete  variable  involves  a  ranking  procedure. 
These  ranks  are  typically  given  integer  values  including  “1”  for  either  the 
lowest  rank  or  the  highest  rank. 

Example.  If  an  Army  installation  is  concerned  with  surface  scaling  dete¬ 
rioration  for  reinforced  concrete  bridge  decks,  they  may  inspect  their  bridge 
decks  visually  and  then  rank  the  appearances  of  the  bridge  decks  with  an 
ordinal  scale. 


Frequency  Distributions 

A  frequency  distribution  is  a  fundamental  method  for  summarizing  data. 

In  order  to  construct  this  distribution,  descriptive  categories  (or  class  inter¬ 
vals)  are  established  so  that  each  data  observation  falls  into  a  single  category. 
The  number  of  observations  that  fall  within  each  category  are  then  counted 
and  tabulated.  The  necessary  descriptive  categories  may  exist  naturally  for 
discrete  variables,  but  they  must  be  fabricated  for  continuous  variables.  The 
appropriate  number  of  class  intervals  is  usually  between  5  and  20,  depending 
on  the  number  of  observations  and  the  apparent  trends  in  the  data  (Spiegel 
1988).  The  frequency  of  occurrences  for  the  various  class  intervals  form  the 
basis  of  a  frequency  distribution. 
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Example.  Frequency  of  occurrences  for  a  continuous  variable,  California 
bearing  ratio  of  a  silty  gravel,  are  developed  in  Table  2-1 . 


Table  2-1 

California  Bearing  Ratio  for  a  Silty  Gravel 

Class  Interval 

Class 

Mark1 

Frequency 

Relative 

Frequency 

Cumulative 

Relative 

Frequency 

25  to  29 

27 

2 

30  to  34 

32 

4 

0.036 

0.055 

35  to  39 

37 

7 

0.064 

0.1 18 

40  to  44 

42 

1 1 

0.100 

0.218 

47 

15 

hruhihib 

0.354 

52 

18 

0.164 

m^i 

55  to  59 

57 

20 

0.182 

0.700 

60  to  64 

62 

14 

0.127 

0.827 

65  to  69 

67 

12 

0.936 

70  to  74 

72 

5 

0.046 

0.982 

75  to  79 

77 

2 

0.018 

1.000 

Total 

110 

1.000 

1  Midpoint  of  class  interval. 

Once  the  frequencies  of  occurrence  are  established,  they  can  be  converted 
to  relative  frequencies  by  expressing  the  number  of  observations  for  each 
category  as  a  proportion  of  the  total  number  of  observations  (Table  2-1). 
Each  relative  frequency,  which  takes  a  value  between  0.0  and  1.0,  represents 
an  estimated  probability  that  an  observation  will  take  on  a  value  within  its 
descriptive  data  category.  For  this  reason,  relative  frequencies  are  often 
referred  to  as  empirical  probabilities  (Freund  and  Wilson  1993).  Once  rela¬ 
tive  frequencies  are  established,  they  can  be  converted  to  cumulative  relative 
frequencies,  which  also  take  values  between  0.0  and  1.0  (Table  2-1)  Cumu¬ 
lative  relative  frequencies  are  often  referred  to  as  empirical  cumulative 
probabilities. 

A  histogram  is  a  graphical  representation  of  either  frequencies  of  occur¬ 
rence  or  relative  frequencies  (Figure  2-1).  These  histograms  are  often  refer¬ 
red  to  as  frequency  distributions  or  relative  frequency  distributions 
depending  on  the  type  of  ordinate  used.  These  histograms  provide  an  easy- 
to-interpret  view  of  both  the  range  of  measurements  and  the  shape  of  the  dis¬ 
tribution  of  measurements.  Important  characteristics  of  the  shape  include 
symmetry  and  data  dispersion.  Overlaying  the  histogram  with  a  cumulative 
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Figure  2-1 .  Relative  frequency  distribution  for  soil  CBR 

relative  frequency  plot  provides  the  viewer  with  an  additional  perspective  of 
the  same  data,  as  shown  in  Figure  2-1. 


Numerical  Descriptive  Statistics 

Although  frequency  distributions  provide  useful  descriptions  of  data, 
numerical  descriptors  are  often  needed  for  quantitative  comparisons  between 
data  sets.  When  numerical  descriptors  are  used  to  represent  frequency  dis¬ 
tributions,  however,  some  information  may  be  lost  (e.g.  the  shape  of  the  fre¬ 
quency  distribution).  For  this  reason,  numerical  descriptors  must  be  used 
with  care  and  with  knowledge  of  any  accompanying  assumptions. 


Central  tendency 

One  of  the  most  useful  single  characteristics  of  a  distribution  is  its  central 
tendency.  Central  tendency  can  be  calculated  several  ways,  including  those 
listed  below  (Freund  and  Wilson  1993).  These  methods  of  calculating  cen¬ 
tral  tendency  are  valid  for  both  populations  of  data  and  samples  of  data. 
Mean  is  the  most  commonly  used  measure  of  central  tendency. 
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a.  Mean.  The  mean  is  the  sum  of  all  values  (assuming  a  finite  population 
size)  divided  by  the  number  of  values.  A  mean  for  a  population  is 
commonly  denoted  fi  and  a  mean  for  a  sample  is  commonly  denoted  y 
(where  the  individual  sample  observations  are  denoted  y).  As  an 
example,  the  mean  of  the  relative  density  data  in  Table  2-2  is  97  33 
(  =  1557.2/16). 


b.  Median.  The  median  is  the  middle  value  when  the  measurements  are 
arranged  from  lowest  to  highest.  The  median  will  be  the  average  of 
two  values  if  the  number  of  measurements  is  even.  The  median  of  the 
relative  density  data  in  Table  2-2  is  97.4. 

c.  Mode.  The  mode  is  the  most  frequently  occurring  measurement.  The 
mode  will  not  be  a  unique  value  if  two  or  more  measurements  occur 
with  the  same  greatest  frequency.  The  mode  will  not  be  defined  if 
each  measurement  occurs  only  once.  The  mode  of  the  relative  density 
data  in  Table  2-2  is  97.5. 
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d.  Geometric  Mean.  The  geometric  mean  is  the  N*  root  of  the  product  of 
N  values.  The  geometric  mean  will  not  be  defined  if  any 
measurements  are  less  than  or  equal  to  zero.  The  geometric  mean  of 
the  relative  density  data  in  Table  2-2  is  97.31. 

e.  Midrange.  The  midrange  is  the  average  of  the  smallest  and  largest 
measured  values.  The  midrange  is  not  often  used  because  it  ignores 
most  of  the  information  provided  by  the  data.  The  midrange  of  the 
relative  density  data  in  Table  2-2  is  97.20. 


Dispersion 

The  variability  or  dispersion  of  a  data  set  is  also  important  and  should  be 
quantified  if  possible.  A  value  for  dispersion  will  provide  an  indication  of 
whether  a  frequency  distribution  is  “broad”  or  “narrow,”  however,  it  will 
not  provide  an  indication  of  symmetry  (or  skewness).  The  simplest  measure 
of  dispersion  is  the  range,  which  is  defined  as  the  difference  between  the 
largest  and  smallest  observed  values.  Similar  to  the  calculation  for  mid¬ 
range,  the  calculation  for  range  ignores  most  of  the  values  in  a  data  set  and 
therefore  has  limited  usefulness. 

Another  measure  of  dispersion,  called  variance,  is  commonly  used.  Vari¬ 
ance  is  represented  by  two  symbols,  depending  on  its  application  (Steel  and 
Torrie  1980):  o2  is  used  if  the  data  set  is  a  population  and  s2  is  used  if  the 
data  set  is  a  sample.  The  population  variance  (assuming  discrete  observa¬ 
tions)  is  defined  as  the  sum  of  squared  deviations  from  the  mean,  divided  by 
the  total  number  of  observations,  N: 


£( yt- 

,2  _  <=! _ 


N 


(1) 


where 


y.  =  value  for  each  observation 
H  =  population  mean 
N  =  number  of  units  in  the  population 

The  numerator  in  Equation  1  is  often  referred  to  as  the  “sum  of  squared 
deviations”  or  the  “sum  of  squares.”  In  order  to  simplify  computations,  the 
sum  of  squares  equation  can  be  transformed  to  a  working  formula: 


N  N 

Efy-/*)2  =  'Ey  i 


<  N  N 

E* 

V  M  ) 


N 


(2) 
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The  calculation  of  sample  variance,  s2,  is  similar  to  the  calculation  of 
population  variance  with  the  exception  that  the  divisor  is  (n-1),  rather  than  N: 


E  (y,-y)2 

2  =  ±1 -  (3) 

n-l 


where 


y,  —  value  for  each  observation 
y  =  sample  mean 
n  =  number  of  units  in  the  sample 

The  necessity  of  using  (n-l)  as  the  divisor  for  the  sample  variance  (Equa¬ 
tion  3),  rather  than  n,  is  a  consequence  of  using  a  sample  statistic,  y,  as  an 
estimate  of  a  population  parameter,  n,  in  its  calculation.  While  population 
parameters  are  fixed  values,  sample  statistics  are  random  variables.  This 
idea  will  be  expanded  upon  in  the  section  tided,  “Sampling  Distributions.” 

For  the  relative  density  data  in  Table  2-2,  variance  is  calculated  as: 


2  0.99  _  n  2 

s  ~  — j-j-  =  0.066  percent 2 


In  addition  to  range  and  variance,  the  dispersion  of  a  data  set  may  be 
quantified  by  its  standard  deviation.  Standard  deviation  is  simply  the  square 
root  of  the  variance,  as  shown  below.  At  times,  the  standard  deviation  is 
convenient  because  its  units  are  the  same  as  the  units  for  the  data  from  which 
it  is  calculated. 


o  =  ^  and  (5) 

s  =  (6) 

where 

a  =  standard  deviation  for  a  population 
s  =  standard  deviation  for  a  sample 

For  the  relative  density  data  in  Table  2-2,  standard  deviation  is  calculated 
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s  =  JP  =  \f 0.066  =  0.257  percent 


(7) 


Since  sample  standard  deviations  and  variances  are  themselves  variables, 
several  independent  estimates  for  these  statistics  may  be  available.  The  pro¬ 
per  method  for  combining  these  estimates  is  to  calculate  a  weighted  average 
for  variance,  with  the  weights  proportional  to  the  respective  degrees  of  free¬ 
dom  [Granley  1969,  ASTM  (1995b)  E  178].  The  calculation  of  weighted 
average  for  variance  is  often  called  “pooling”  the  variance  estimates: 

£  Ks?)(n-l)] 

pooled  sample  variance  =  — -  (8) 

£  ni~M 

;=i 


where 


n{  =  number  of  tests  in  project  I 
M  =  number  of  projects 

Pooling  is  essentially  a  weighted  average  with  the  degrees  of  freedom  in  each 
project  (nrl)  as  the  weighting  factor.  If  the  number  of  tests  within  each  pro¬ 
ject  are  the  same,  the  pooled  value  is  equal  to  the  straight  average. 

Finally,  the  dispersion  of  a  data  set  can  be  quantified  by  a  coefficient  of 
variation  (CV),  as  shown  below.  The  coefficient  of  variation  can  be 
described  as  the  standard  deviation  expressed  as  a  percentage  of  the  mean. 
This  measure  of  dispersion  has  advantages  at  times  because  the  magnitude  of 
the  mean  and  the  units  of  measure  are  factored  out.  Using  CV,  distributions 
representing  data  with  different  units  and  distributions  representing  data  with 
different  means  can  all  be  compared  in  terms  of  their  dispersion.  Equations 
for  computing  CV  for  both  populations  and  samples  are  shown  below. 

CV  =  —  X  100%  or  (9) 


CV  =  -4-  x  100%  (10) 

y 

For  the  relative  density  data  in  Table  2-2,  coefficient  of  variation  is  cal¬ 
culated  as: 


CV  =  °-25—  x  100%  =  0.26% 
97.33 


(ID 
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Skewness  is  a  measure  of  the  degree  of  symmetry  for  a  frequency  distri¬ 
bution.  If  the  right  tail  of  a  frequency  distribution  extends  farther  from  the 
central  maximum  than  the  left  tail,  the  distribution  is  said  to  be  skewed 
positively  (Figure  2-2).  If  the  left  tail  of  a  frequency  distribution  extends 
farther  from  the  central  maximum  than  the  right  tail,  the  distribution  is  said 
to  be  skewed  negatively  (Spiegel  1990).  The  moment  coefficient  of  skew¬ 
ness  for  a  population  (Pj)  is  calculated  as: 


ECy,--/*)3 


The  moment  coefficient  of  skewness  for  a  sample  (b,)  is  calculated  as 
(SAS  1988): 


b 


l 


Ety-y)3 

-  x  — - 

(«-!)(« -2)  s3 


(13) 


If  a  distribution  is  positively  skewed,  its  moment  coefficient  of  skewness 
will  be  positive.  If  a  distribution  is  negatively  skewed,  its  moment  coeffi¬ 
cient  of  skewness  will  be  negative.  If  a  distribution  is  symmetric,  its  moment 
coefficient  of  skewness  will  be  equal  to  zero.  By  raising  the  standard  devi¬ 
ation  in  the  denominator  to  a  power  of  three,  the  moment  coefficient  of 
skewness  becomes  unitless. 

For  the  relative  density  data  in  Table  2-2,  coefficient  of  skew  is  calculated 
as: 


bi 


16  x  -0.18 
15-14  (0.257)3 


-0.83 


(14) 


Discrete  Probability  Distributions 


Relative  frequency  distributions  were  described  previously  as  empirical 
probability  distributions.  Relative  frequency  distributions  are  considered 
empirical  because  they  are  usually  constructed  from  sample  values  Relative 
frequency  distributions  can  also  be  called  “discrete”  distributions  because 

they  are  constructed  from  data  that  has  been  categorized  into  discrete  class 
intervals. 
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Relative  Frequency  Relative  Frequency  Relative  Frequency 


Figure  2-2 


Class  Mark 


Class  Mark 


Relative  frequency  distributions  with  various  skew: 

(a)  symmetric,  (b)  positive  skew,  and  (c)  negative  skew 
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Contrary  to  relative  frequency  distributions,  which  are  constructed  from 
sample  data,  the  construction  of  probability  distributions  is  based  on  assump¬ 
tions  related  to  either  the  shape  of  the  distribution  or  the  probability  of  occur¬ 
rence  for  each  discrete  value.  Probability  distributions  are  therefore 
theoretical,  rather  than  empirical,  and  represent  theoretical  characteristics  of 
a  population  of  data. 


Binomial  distribution 

A  relatively  simple  and  useful  discrete  probability  distribution  is  the  bino¬ 
mial  distribution.  This  distribution  is  applicable  for  statistical  problems  that 
include  “repeated  trials.”  In  these  situations,  we  are  interested  in  the  prob¬ 
ability  of  getting  x  successes  in  n  trials.  In  order  for  the  binomial  distribu¬ 
tion  to  apply,  the  following  assumptions  must  also  be  true  (Johnson  1994). 

a.  There  are  only  two  possible  outcomes  for  each  trial  (usually  referred  to 
as  “success”  and  “failure”). 

b.  The  probability  of  success  (p)  is  the  same  for  each  trial. 

c.  The  n  trials  are  independent. 

Let  Y  be  the  random  variable  that  represents  a  set  of  all  possible  successes 
m  n  trials.  The  probability  of  the  occurrence  of  y  successes,  which  is  repre¬ 
sented  as  p(y),  can  be  calculated  for  each  possible  value  of  y.  The  binomial 
distribution  can  then  be  constructed  by  plotting  p(y)  as  a  function  of  the  dis¬ 
crete  values  of  y.  In  order  to  emphasize  the  fact  that  p(y)  is  calculated  for  a 
binomial  distribution,  p(y)  is  often  denoted  by  b(y;  n,  p). 

b(yi  n,  p)  =  P>(1  -pf-y  for  y  =  0,1,2 . ,«  (15) 


where 


M  =  -  »! 

\3v  y!  (n-y)!  (16) 

Cumulative  probabilities  of  the  binomial  distribution  are  often  needed  for 
statistical  applications.  The  cumulative  binomial  probability  B  can  be 
calculated  as  follows: 


k 

B(k;  n,p)  =  £  b(y\  n,  p)  for  k  =  0,1,2 ,...,« 
y=0 


(17) 
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The  mean  and  the  variance  for  the  binomial  random  variable  y  can  be 
calculated  without  constructing  the  probability  distribution,  provided  that  the 
two  distribution  parameters  n  and  p  are  known. 


1 u  =  np 

and 

(18) 

a2  =  np{  1 

~P) 

(19) 

If  n  is  large,  calculating  the  binomial  probabilities  and  the  cumulative 
binomial  probabilities  can  be  cumbersome.  Statistical  software  packages  are 
often  equipped  to  handle  these  calculations.  In  addition,  many  textbooks 
include  an  appendix  with  cumulative  binomial  probabilities  tabulated  for  a 
range  of  values  for  each  of  y,  n,  and  p. 

Example.  Assume  that  a  contractor  is  removing  concrete  cores  from  a 
rigid  pavement.  The  contractor  needs  four  cores  that  remain  intact  and  are 
of  sufficient  length  for  compression  testing  (successes).  He/she  knows  from 
experience  that  only  about  one-half  of  the  cores  attempted  are  of  sufficient 
quality  (successes),  so  p=0.5.  If  the  contractor’s  crew  is  cutting  12  cores 
(n=12),  his/her  binomial  distribution  for  this  situation  is  shown  in  Fig¬ 
ure  2-2.  The  probability  that  the  contractor  will  obtain  three  or  fewer  quality 
cores  is  only  7  percent.  The  probability  that  the  contractor  will  obtain  the 
necessary  four  quality  cores,  or  even  more,  is  approximately  93  percent. 


Continuous  Probability  Distributions 


Probability  distributions  can  also  be  constructed  for  continuous  random 
variables,  which  can  assume  an  infinite  number  of  different  values  (any  value 
within  an  interval).  Continuous  probability  distributions  for  these  continuous 
variables  have  the  following  characteristics  (Freund  and  Wilson  1993). 

a.  The  graph  of  the  distribution  is  usually  a  smooth  curve.  This  is  in  con¬ 
trast  to  the  histogram  type  of  distribution  for  discrete  variables. 

b.  The  curve  is  described  by  an  equation, /(y),  called  the  distribution 
function.  This  function  corresponds  to  p(y)  for  discrete  variable 
distributions. 

c.  The  total  area  under  the  curve  is  one.  This  corresponds  to  the  sum  of 
all  possible  probabilities  for  discrete  variables. 

d.  The  probability  of  a  random  variable  taking  a  value  within  an  interval 
(a,  b)  is  equal  to  the  area  between  the  distribution  curve  and  the  hori¬ 
zontal  axis  within  the  interval  (a,  b).  This  corresponds  to  adding  prob¬ 
abilities  for  discrete  variables. 
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Normal  distribution 


The  normal  distribution  is  a  continuous  distribution  that  has  proved  to  be 
useful  in  many  applications.  The  normal  distribution  was  developed  in  the 
eighteenth  century  when  scientists  observed  regularity  in  errors  of  measure¬ 
ment,  caused  by  laws  of  chance.  They  called  the  family  of  distributions  the 
normal  curve  of  errors”  and  fit  the  following  continuous  equation,  often 
referred  to  as  the  normal  probability  density  function  (Johnson  1994). 


fly)  =  — - —  exp 
a  JTk 


2  o2 


(20) 


where 


-ao  <  y  <  CO 


The  position  and  breadth  of  the  resulting  normal  probability  density  func¬ 
tion  is  dependent  on  only  two  parameters:  mean  (p)  and  standard  deviation 
(a).  The  curve  is  always  bell-shaped  and  symmetric  about  the  mean  so  its 
degree  of  skewness  is  always  zero.  The  breadth  of  the  curve  increases  with 
increasing  a.  This  effect  is  shown  in  Figure  2-3  for  functions  with  a  mean  of 
0  and  standard  deviations  of  1/2,  1,  and  2.  While  the  breadth  of  normal 
distributions  can  change,  the  following  rules  of  thumb  remain  applicable  for 
data  dispersion  (Willenbrock  1974b). 

a.  Approximately  68.0  percent  of  the  data  falls  within  plus  or  minus  one 
standard  deviation  from  the  mean  (fi  ±  a). 

b.  Approximately  95.5  percent  of  the  data  falls  within  plus  or  minus  two 
standard  deviations  from  the  mean  (fi  ±  2a). 

c.  More  than  99.7  percent  of  the  data  falls  within  plus  or  minus  three 
standard  deviations  from  the  mean  (fi  ±  3a). 


Standard  normal  distribution 


The  calculus  required  to  integrate  under  normal  distribution  curves  can  be 
quite  complex,  so  probabilities  related  to  areas  under  these  density  functions 
are  seldom  calculated  without  the  aid  of  computer  software.  If  software  is 
not  convenient  for  a  particular  application,  tables  representing  a  standard 
normal  distribution  can  be  used  to  facilitate  manual  calculations.  The  stan¬ 
dard  normal  distribution  has  a  mean  of  zero  and  a  standard  deviation  of  one, 
resulting  in  the  probability  density  function  shown  below. 


fly) 
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Figure  2-3.  Discrete  binomial  distribution 


where 

-oo  <  y  <  oo 

For  all  normal  distributions,  the  degree  of  skewness  is  zero. 

The  random  variable  associated  with  the  standard  normal  distribution  is 
usually  represented  by  the  letter  Z  [={z}].  Areas  under  a  normal  distribution 
for  any  variable  Y  [={y}]  can  be  determined  by  converting  a  specific  value  y 
to  a  value  z  according  to  the  equation  below.  The  calculated  value  z  can  then 
be  used  as  input  for  standard  normal  distribution  tables,  such  as  that  shown 
in  Table  A1  (Appendix  A),  to  determine  the  probability  of  finding  measure¬ 
ments  less  than  or  greater  than  y. 


o 

Example.  Assume  that  a  specific  mixture  of  portland  cement  concrete  is 
required  to  have  a  compressive  strength  greater  than  or  equal  to  34  MPa  and 
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that  the  concrete  supplier  involved  has  historically  reported  coefficients  of 
variation  for  strength  in  the  range  of  20  percent.  If  the  average  compressive 
strength  for  the  delivered  concrete  is  40  MPa,  we  can  use  the  “Z”  tables  to 
calculate  the  probability  that  a  measured  compressive  strength  will  fall  below 
the  specified  34  MPa.  From  the  given  values  for  coefficient  of  variation  and 
mean,  the  standard  deviation  of  compressive  strengths  is  calculated  to  be 
8  MPa  (0.2  x  40  MPa).  The  z-value  that  corresponds  to  a  y-value  of  34  MPa 
equals  -0.75.  Since  the  standard  normal  distribution  is  symmetric,  we  know 
that  the  probability  of  obtaining  a  value  of  z  less  than  -0.75  is  equal  to  the 
probability  of  obtaining  a  value  of  z  greater  than  +0.75.  Therefore,  using 
the  standard  Z  table  (Table  Al)  and  interpolation,  we  can  determine  that  the 
probability  of  measuring  a  compressive  strength  smaller  than  the  specified 
34  MPa  equals  0.2266  or  23  percent.  The  assumed  normal  distribution  for 
compressive  strength,  along  with  a  cumulative  probability  distribution,  are 
plotted  in  Figure  2-4.  As  expected,  the  cumulative  probability  at  34  MPa  is 
approximately  0.23. 


Normal  Variate  Y 


Figure  2-4.  Effect  of  standard  deviation  on  the  breadth  of  a  normal 
distribution 
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Beta  distribution 


The  normal  distribution  is  not  the  most  appropriate  representation  for  all 
applications.  Some  random  variables  are  represented  more  accurately  by 
distributions  with  finite  limits.  For  example,  California  bearing  ratio  (CBR) 
measurements  for  a  high-plasticity  clay  would  have  a  mean  close  to  zero  and 
a  relatively  high  standard  deviation.  However,  the  real-life  measurements 
could  never  have  a  value  less  than  zero.  In  cases  such  as  this  and  for  skewed 
distributions,  a  beta  distribution  may  be  more  applicable  than  a  normal  distri¬ 
bution.  Beta  probability  distributions  are  continuous  and  they  have  a  finite 
range  from  y=a  to  y=b.  These  probability  distributions  can  be  described  by 
the  density  function  shown  below  (Harr  1987). 


fly)  =  C(y-aT(b-yf 


(23) 


where 


a  >  -1,  P  >  -1,  and  (24) 


C  =  — (a  +  P  +  !)• —  (G  and  p  are  integers)  (25) 

a!p!(6-a)“+p+1 

Given  a,  b,  a,  and  fl  for  a  particular  beta  distribution,  the  mean  and  the  vari¬ 
ance  of  the  distribution  can  be  calculated  as  shown  below. 


(X  +  l  /u  \ 

/*  =  «+ - iT-T(Z,_a) 

a  +  p  +2 


(26) 


a2  =  (fr-Q)2(a  +  l)(P  +  l) 

(a  +P  +  2)2(a  +  P  +  3) 


(27) 


Given  a  range,  a  mean,  and  a  standard  deviation  for  a  random  variable  Y, 
values  of  a  and  p  for  the  appropriate  beta  distribution  can  be  calculated  as 
shown  below. 


a  =  — (1  -P)  -  (1 +P)  and 

Q2 


(28) 


P  =  £li-(a+2) 


(29) 
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where 


P  =  ^~fl) 
(b-a) 


and 


(30) 


Q 


g 

(b-a) 


(31) 


TJeA.becudiStribUt!°n  Can  assume  many  shaPes’  depending  on  the  values  of 
a  and  fj.  Shapes  include  uniform,  triangular,  skewed  left  or  right  and  even 
a  shape  resembling  a  normal  distribution  with  finite  limits,  as  shown  in 
Figure  2-5.  The  beta  distribution  is  similar  to  the  normal  distribution  when  a 
J**  3'  ,Due  t0  ^  complexities  of  beta  distributions,  they  are  not  easily 
tabulated  into  a  standard  form  such  as  the  “Z"  tables  for  the  normal  distribu- 
tion.  However,  computer  software  is  available  to  help  with  the  construction 
ot  beta  distributions  and  for  calculating  probabilities  associated  with  these 
distributions.  Harr  (1987)  provides  the  source  code  for  this  type  of  software. 


Compressive  Strength  (MPa) 
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Figure  2-5.  Concrete  compressive  strength  with  an  assumed  normal 
distribution 
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Example.  Assume  that  a  pavement  engineer  is  interested  in  load  transfer 
between  slabs  that  rely  solely  on  aggregate  interlock.  The  engineer  learns 
from  literature  that  for  his/her  slab  configuration  the  average  load  transfer  is 
approximately  20  percent.  He/she  also  learns  that  the  standard  deviation  for 
load  transfer  is  typically  about  10  percent.  Due  to  the  fact  that  measurements 
cannot  possibly  fall  below  zero  percent  or  above  50  percent  and  due  to  the 
fact  that  the  distribution  has  been  reported  to  be  skewed,  the  engineer  may 
decide  that  a  beta  distribution  would  provide  a  better  representation  of  the 
load  transfer  measurements  than  a  normal  distribution.  The  beta  distribution 
and  the  cumulative  probability  distribution  for  this  problem  were  calculated 
with  Harr’s  (1987)  software,  as  shown  in  Figure  2-6.  If  the  engineer  was 
interested  in  the  likelihood  that  load  transfer  would  be  less  than  15  percent,  a 
quick  measurement  taken  from  the  cumulative  probability  plot  would  indicate 
that  this  should  occur  approximately  35  percent  of  the  time. 


Figure  2-6.  Shapes  of  the  beta  distribution  (Harr  1987) 


Sampling  Distributions 


The  example  problems  used  up  to  this  point  have  included  data  sets  that 
were  considered  to  be  “populations.”  When  dealing  with  populations,  the 
distribution  parameters  can  be  considered  to  be  fixed  values.  When  taking 
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samples  from  a  population,  however,  the  statistics  that  describe  the  distribu¬ 
tion  of  sample  values  cannot  be  considered  to  be  fixed.  Different  samples 
taken  from  the  same  population  can  generate  different  statistics.  A  statistic 
computed  from  a  random  sample  is  therefore  a  random  variable.  This  inher¬ 
ent  variability  for  sample  statistics  must  be  considered  whenever  samples  are 
used  to  make  inferences  about  a  population. 

All  the  possible  values  of  a  sample  statistic  can  be  described  by  a  prob¬ 
ability  distribution  for  the  statistic,  often  called  a  sampling  distribution. 
Characteristics  of  a  sampling  distribution  can  be  related  to  characteristics  of 
the  population  from  which  the  samples  were  drawn. 

Assume  samples  that  include  n  observations  each  are  drawn  from  a  popu¬ 
lation,  Y  t-{y}],  with  mean  fi  and  variance  o2.  As  the  number  of  samples 
drawn  approaches  infinity,  the  distribution  of  sample  means,  f„  will  have  a 
mean  that  approaches  n  and  a  variance  that  approaches  rf/n,  as  shown 
below. 


Ky()  =  f1  and 


(32) 


o2(y,) 


(33) 


This  statement  makes  intuitive  sense  as  one  would  expect  the  sample 
means  to  cluster  around  the  population  mean  and  one  would  expect  the  vari¬ 
ance  of  sample  means  to  be  less  than  the  variance  of  individual  observations 
Sample  means  are  referred  to  as  “unbiased”  estimates  of  the  population 
mean. 


Distribution  for  sample  means 

An  ^additional  important  characteristic  of  the  distribution  of  sample  means 
is  the  central  limit  theorem.  ”  This  theorem  states  that  the  distribution  of 
sample  means  can  be  closely  approximated  by  the  normal  distribution, 
regardless  of  the  population  from  which  the  samples  are  drawn.  The  size  of 
the  samples  required  to  validate  this  theorem  is  dependent  on  the  shape  of  the 
parent  population.  If  the  population  resembles  normality,  sample  sizes  of 
10  or  more  should  be  sufficient.  Sample  sizes  of  30  or  more  should  be  suffi¬ 
cient  for  populations  of  any  other  shape  (Freund  and  Wilson  1993). 

Grant  and  Leavenworth  (1972)  stated  that  even  if  sample  size  (n)  is  small 
the  distribution  of  the  means  of  the  samples  can  be  very  close  to  normal  if 
the  number  of  samples  is  sufficiently  large.  This  theory  holds  true  even  if 
the  parent  population  is  far  from  normal.  Grant  and  Leavenworth  (1972) 
reported  on  a  study  in  which  1,000  samples  of  size  n  =  4  were  taken  from 
two  bowls:  one  containing  numbered  tags  from  a  rectangular  distribution 
population  and  one  containing  numbered  tags  from  a  triangular  distribution 
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population.  Neither  of  the  original  populations  resembled  a  normal  distribu¬ 
tion,  however,  the  distribution  of  sample  means  in  each  case  was  normal. 
They  reported: 

“The  great  practical  importance  of  the  normal  curve  arises  even  more 
from  its  use  in  sampling  theory  than  from  the  fact  that  some  observed  distri¬ 
butions  are  described  by  it  well  enough  for  practical  purposes.  Of  great 
practical  significance  is  the  fact  that  distributions  of  averages  of  samples  tend 
to  be  approximately  normal  even  though  the  samples  are  drawn  from  non¬ 
normal  populations.” 

When  the  distribution  of  sample  means  can  be  assumed  to  be  normal  in 
shape,  the  standard  normal  variate  Z  can  be  used  as  a  problem-solving  tool. 
Recall  the  transformation  of  a  random  variable  from  a  population  to  a  stan¬ 
dard  normal  variable. 


z 


=  1UL 
o 


(34) 


Similarly,  a  sample  mean  can  be  transformed  from  a  population  to  a 
standard  normal  variable: 


ziyi) 


(35) 


The  denominator  in  Equation  30  represents  the  standard  deviation  for  the 
distribution  of  sample  means,  which  can  be  calculated  from  the  variance  of 
sample  means  (Equation  28),  as  shown  below.  The  standard  deviation  of 
sample  means  is  often  referred  to  as  the  “standard  error  of  the  mean.” 


o(y,) 


(36) 


Example.  Suppose  that  an  engineer  is  planning  to  construct  a  parking  lot 
in  an  area  that  has  a  variable  water  table  due  to  the  presence  of  various  types 
of  soil.  Suppose  also  that  he/she  knows  from  historical  data  that  the  depth  to 
water  around  the  site  should  follow  a  near-normal  distribution.  The  mean 
and  standard  deviation  of  the  depth  to  water  is  2  meters  and  0.5  meters, 
respectively.  If  the  engineer  is  interested  in  the  probability  that  a  sample  of 
16  measurements  will  have  a  mean  less  than  1.75  meters,  he/she  can  use  the 
Z-table  as  shown  below. 


22 


Chapter  2  Statistical  Concepts 


(37) 


«©  -W-* 


nfr° “  standard  z  table  (Appendix  A),  the  engineer  could  see  that  the 

getting  a  z‘value  Iess  than  -2  is  approximately  2  percent 
0.0228).  The  engineer  could  then  feel  confident  that  the  mean  of  a 
16-rephcate  sample  would  very  rarely  be  less  than  1.75  meters. 

rp,  F(*  this  same  problem,  assume  the  sample  consists  of  only  4  replicates. 
Then  the  corresponding  £-value  is  as  shown  below. 


*(?,)  =  -25~2  =  -i 
0.5 


From  the  standard  z-tables,  the  engineer  could  see  that  the  probability  of 
getting  a  z-value  less  than  -1  is  approximately  16  percent  (0.1587)  With  the 
decrease  in  sample  size  from  n=16  to  n=4,  the  standard  deviation  for  the 
distribution  of  sample  means  doubled.  As  a  result,  the  estimated  probability 
tor  obtaining  a  given  deviation  from  the  mean  increased. 

The  effects  of  sample  size  on  the  probability  distribution  for  sample 
means  is  illustrated  in  Figure  2-7.  In  this  figure,  the  distributions  for  sample 
means  are  shown  to  be  more  narrow  (or  tighter)  than  the  standard  normal 
distribution  from  which  they  were  derived.  As  sample  size  increases,  the 
distributions  for  sample  means  became  progressively  more  narrow. 

Distribution  for  sample  variances 

When  samples  of  n  observations  each  are  drawn  from  a  population  with 
mean  M  and  variance  a2,  the  sample  variances  s2  can  be  used  to  construct  a 
^?UwCy  dlst"^Utl0n:  The  theoretical  relationship  between  the  mean  of  this 
distribution,  n(s  ),  and  the  true  population  variance,  o2,  is  shown  below. 

Ksh  =  —a2 


For  small  sample  sizes,  the  mean  of  the  distribution  of  sample  variances 
would  be  less  than  the  true  population  variance.  However,  as  the  sample 
size  71  increases,  the  mean  of  the  distribution  of  sample  variances  becomes  a 
good  estimate  for  the  true  population  variance.  This  inequality  between  u(s2) 
and  o2  causes  sample  variances  to  be  termed  “biased”  estimates  of  the  popu¬ 
lation  variance.  In  order  to  compensate  for  this  bias,  sample  variances  are 
calculated  with  ( n-1 )  as  the  divisor,  rather  than  n  (as  shown  below).  As  the 
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Load-Transfer  (percent) 


Figure  2-7.  Load-transfer  between  pavement  slabs  with  an  assumed  beta 
distribution 


sample  size  ( n )  increases,  the  effect  of  the  bias,  along  with  the  magnitude  of 
the  bias  correction,  decreases. 

unbiased  s2  =  — —  ( biased  s2)  (40) 

n- 1 

Recall  the  difference  between  the  calculations  for  population  variance  o2 
and  sample  variance  s2,  as  presented  previously.  While  the  divisor  for  popu¬ 
lation  variance  was  N,  the  divisor  for  sample  variance  was  n-1 ,  as  shown  in 
Equations  41  and  42,  respectively. 

n 

5D  (yrri2  (4i) 

=  /-! _ 
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(42) 


E  (y,--y)2 

2  _  _M _ 

n- 1 


The  term  /i  1  in  these  calculations  for  sample  variance  can  be  referred  to 
as  the  degrees  of  freedom  for  the  sample  variance  statistic.  In  general,  the 
degrees  of  freedom  for  a  statistic  ( v)  is  defined  as  the  number  of  independent 
observations  (n)  minus  the  number  of  population  parameters  (m)  that  must  be 
estimated  from  sample  observations  in  order  to  calculate  the  statistic  In 
symbols,  v  =  n-m  (Spiegel  1990). 

As  shown  previously ,  the  calculation  for  sample  variance  requires  an  esti¬ 
mate  for  the  population  mean  (fx),  which  is  provided  by  the  sample  mean  (y). 
Therefore,  the  degrees  of  freedom  for  sample  variance  can  be  calculated  as 
v  =  n-m  =  n-1. 


T  distribution 

In  previous  applications  of  the  Z  tables,  the  standard  deviation  of  the 
population  (a)  was  assumed  to  be  known.  In  most  practical  situations,  how¬ 
ever,  a  is  not  known  and  the  only  available  measure  of  the  population  stan¬ 
dard  deviation  is  the  sample  standard  deviation  (s).  In  order  to  make 
continued  use  of  the  standard  normal  tables  for  sample  mean  distributions 
one  would  like  to  substitute  5  for  a  in  the  calculation  of  Z.  This  substitution 
can  be  made,  but  the  resulting  distribution  is  not  exactly  normal  and  has  to  be 
handled  differently  than  the  Z  distribution.  If  s  substitutes  for  o,  the  distri¬ 
bution  is  referred  to  as  a  t  distribution  (or  Student’s  t).  The  calculation  of  t 
is  performed  as  shown  in  Equation  43.  The  t  distribution  is  similar  to  normal 
m  that  it  is  bell-shaped  and  symmetric  about  its  mean.  Relative  to  normal, 
however,  the  t  distribution  is  broader  (has  fatter  tails).  This  breadth  reflects 
increased  dispersion  for  the  variate,  which  is  caused  by  the  uncertainty  of 
using  a  sample  statistic  to  estimate  the  population  standard  deviation.  As  the 
sample  size  increases,  confidence  in  the  estimate  of  the  population  standard 
deviation  increases,  causing  the  t  distribution  to  become  more  narrow  and  to 
approach  normality  (Baecher  1987),  as  shown  in  Figure  2-8. 


t 


=  y-A* 


5 


(43) 


The  probability  density  function  that  describes  the  t  distribution  is  tedious 
to  integrate  manually,  so  problems  are  typically  solved  with  the  help  of  com¬ 
puter  software  or  standard  t  tables,  such  as  that  shown  in  Table  A2 
(Appendix  A). 
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Figure  2-8.  Effect  of  sample  size  on  the  probability  distribution  for 
sample  means 


Chi-squared  (x2)  distribution 

The  distribution  provides  a  method  for  comparing  sample  variances 
with  population  variances.  If  is  calculated  for  one  or  more  samples  taken 
randomly  from  a  parent  normal  population,  then  also  becomes  a  random 
variable.  Chi-square  can  be  defined  for  each  sample  as  shown  below. 


£  ( yry )2 


Y2  =  _  =  ±i _ 

2 

a  a 


where 

yi  =  each  value  within  a  sample 
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y  -  the  calculated  mean  for  the  sample 
SS  =  sum  of  squared  deviations  for  the  sample  data 
o2  =  the  variance  for  the  parent  population 
«  =  the  number  of  values  obtained  for  the  sample 

The  relationship  between  sample  variance  and  population  variance  Equa¬ 
tion  3),  contained  within  the  variable,  can  be  seen  more  clearly  with  a  dif¬ 
ferent  expansion  on  the  sum-of-squares  term.  Recall  the  definition  for  the 
unbiased  estimate  of  population  variance,  calculated  from  a  sample  of  data- 
s  -  SS/(n-l).  Substitution  of  this  relationship  for  the  sum-of-squares  in 
Equation  44  provides  the  following  formula  for  / 


(h-1)s2 


(45) 


The  probability  density  function  that  describes  a  f  distribution  is  tedious 
to  integrate  manually,  so  software  or  published  tables,  such  as  that  shown  in 
Table  A3  (Appendix  A),  are  often  used.  Similar  to  the  variances  from  which 
they  are  calculated,  values  cannot  be  negative.  The  jf  distribution  can  be 
other  than  normal  and  its  shape  is  dependent  on  the  number  of  values 
obtained  for  the  sample  ( n ),  as  shown  in  Figure  2-9. 

A  few  important  characteristics  of  the  chi-square  distribution  are  listed 
below  (Freund  and  Wilson  1993). 

a.  The  mean  of  the  /  distribution  is  equal  to  the  degrees  of  freedom,  v 
(—n-1),  and  the  variance  is  equal  to  2v. 

b.  For  large  values  of  v  (usually  greater  than  30),  the  /  distribution  can 
be  approximated  by  the  normal  distribution  with  a  mean  and  variance 
as  described  in  “a.”  Thus,  in  these  cases,  the  standard  normal 
distribution  can  be  used  as  shown  below. 


i/2\  (46) 


c.  The  accuracy  with  which  /  represents  SS/o2  for  samples  is  not  signifi¬ 
cantly  affected  when  the  parent  population  deviates  slightly  from  nor¬ 
mal.  Severe  departures  from  normality,  however,  can  affect  the 
reasonableness  of  using 

Example.  Suppose  that  an  agency  is  monitoring  asphalt  cement  content 
for  a  paving  job.  If  a  sample  of  10  measurements  had  a  variance  of 
0.04  percent2  and  the  variance  had  historically  been  0.03  percent2,  what 
would  be  the  probability  of  obtaining  a  sample  variance  of  0.04  percent  or 
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variable  t  or  variable  z 


Figure  2-9.  Effect  of  sample  size  on  the  probability  distribution  for 
student's  t 


higher?  The  calculation  is  illustrated  in  Equation  42.  Degrees  of 
freedom  is  nine  (10-1). 

2  __  (10-1)004  =  12 
*  0.03 

From  the  distribution  (Table  A3)  and  linear  interpolation,  the  prob¬ 
ability  of  obtaining  a  /  value  of  12  or  higher  is  found  to  be  approximately 
22  percent  (0.2227).  The  evidence  should  suggest  to  the  agency  that  this 
sample  variance  is  not  highly  unlikely.  Based  on  this  data  alone,  there  is  no 
reason  to  judge  that  the  manufacturing  process  for  asphalt  concrete  is  out  of 
control. 
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F  distribution 


The  ^-statistic  is  defined  as  the  ratio  of  two  independent  j?  random  vari¬ 
ables, ^each  divided  by  their  respective  degrees  of  freedom.  Assume  that  s2 
and  s2  are  the  variances  of  two  independent  samples  of  size  n,  and  n2 
respectively.  If  these  samples  are  taken  from  the  same  population  or  from 
two  different  populations  with  variances  that  are  expected  to  be  equivalent 
then  F  can  be  calculated  as  shown  below. 


F  = 


(48) 


The  choice  for  which  variance  estimate  to  place  in  the  numerator  is  some¬ 
what  arbitrary.  Typically,  reference  tables  for  the  distribution  of  the 
F-statistic  assume  that  the  larger  variance  is  used  as  the  numerator.  Since  the 
samples  from  which  variances  are  calculated  are  selected  randomly,  F  is  also 
a  random  variable.  Similar  to  the  /  probability  distribution,  the  f’ distribu¬ 
tion  is  defined  only  for  positive  values.  The  shape  of  the  F  distribution  is 
dependent  on  the  degrees  of  freedom  for  each  of  the  sample  sets  from  which 
variance  was  calculated. 


The  probability  density  function  for  the  F  distribution  is  complex  and 
would  be  difficult  to  integrate  manually.  Similar  to  other  sampling  distribu¬ 
tions,  the  use  of  software  or  reference  tables,  such  as  that  shown  in  Table  A4 
(Appendix  A),  are  essential. 

The  F  distribution  is  nonsymmetric,  particularly  when  one  or  both  of  the 
distribution  parameters  has  a  low  value.  Figure  2-10  shows  F  distributions 
when  the  number  of  replicates  for  each  sample  are  equal  and  consist  of 
values  ranging  from  5  to  50.  Figure  2-1 1  shows  the  effect  on  the  F  distribu¬ 
tion  when  the  number  of  replicates  for  one  of  the  samples  drops  from  50 

The  F  distribution  is  used  to  conduct  analysis-of-variance  (ANOVA)  pro¬ 
cedures,  which  will  be  demonstrated  in  the  next  chapter.  The  F  distribution 
can  also  be  used  to  make  simple  inferences  concerning  the  equality  of  vari¬ 
ances  for  two  normal  populations. 

Example.  Suppose  that  two  construction  crews  (A  and  B)  are  placing 
asphalt  concrete  surface  course  mixtures.  Suppose  further  that  21  randomly 
located  cores  were  removed  from  crew  A’s  pavement  while  11  randomly 
located  cores  were  removed  from  crew  B’s  pavement.  The  standard  devi¬ 
ation  for  surface  course  thickness  was  determined  to  be  5  mm  and  3  mm, 
respectively.  If  the  two  crews  should  be  paving  with  equal  uniformity  in 
thickness,  what  is  the  probability  of  finding  a  difference  in  standard  devia¬ 
tions  as  large  or  larger  than  that  measured? 
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0  5  10  15  20  25  30 

Chi  Square 


Figure  2-10.  Effect  of  sample  size  on  the  chi-square  distribution 


The  F-statistic  is  calculated  as  shown  in  Equation  49.  The  value  2.8  is 
compared  to  values  provided  in  the  standard  F  distribution  table,  with 
degrees  of  freedom  equal  to  20  and  10,  F  (20,10).  Using  linear  interpolation 
within  Table  A4,  the  probability  of  obtaining  a  ratio  of  variances  this  large  is 
found  to  be  less  than  5  percent  (0.0488).  The  variances  can  therefore  be 
considered  significantly  different  with  a  high  degree  of  confidence.  The 
project  manager  can  feel  confident  that  the  pavement  placed  by  crew  B  has  a 
more  uniform  thickness  than  the  pavement  placed  by  crew  A. 


^  2 

F  =  —  =  —  =  2.8  (49) 

si  (3)2 
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Figure  2-1 1 .  Effect  of  sample  sizes  on  the  F  distribution 
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Probability 


Figure  2-12.  Effect  of  unequal  sample  sizes  on  the  F  distribution 
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3  Analyzing  Pavement 
Materials  Data 


This  chapter  provides  a  review  of  techniques  needed  for  studying  the 
variability  of  construction  data. 


Distribution  Shapes 


Smooth  curve  approximations  for  the  frequency  distribution  histograms 
tor  experimental  data  can  take  various  forms  (Willenbrock  1974),  such  as 
those  shown  in  Figure  3-1.  The  symmetrical,  or  bell-shaped  curve,  is  char¬ 
acteristic  of  the  normal  distribution.  Knowledge  of  whether  frequency  distri¬ 
bution  histograms  are  skewed  is  important  because  many  statistical  tests  used 
tor  the  interpretation  of  data  are  based  on  the  assumption  of  symmetry,  or 
even  normality.  The  J-shaped  and  U-shaped  distributions  shown  in  Fig¬ 
ure  3-1  are  not  common  in  cases  for  construction  materials  unless  data  has 
been  manipulated.  For  example,  compressive  strength  data  for  concrete  can 
be  reverse  J-shaped  if  all  strengths  less  than  the  specificied  strength  are 
discarded. 

Symmetric  and  skewed  distributions  are  examples  of  unimodal  (single 
peak)  distributions.  Distributions  can  also  be  bimodal  or  multimodal  if  more 
an  one  peak  exists.  These  characteristics  have  serious  implications  for  sta¬ 
tistical  analyses  because  the  sample  of  data  has  probably  not  been  obtained 
from  a  single  population.  In  a  bimodal  distribution,  each  of  the  two  peaks 
may  represent  the  mean  of  a  different  population. 

The  overriding  pattern  for  the  frequency  distribution  histograms  for  con¬ 
struction  materials  (e.g.  quality  control  data)  is  that  they  tend  toward  bell- 
curves>  rather  than  J-shaped,  U-shaped,  or  multimodal  curves 
(Willenbrock  1974).  When  variability  among  data  is  caused  by  the  cumula¬ 
tive  effect  of  a  large  number  of  small  perturbations  or  errors,  the  resulting 
frequencies  of  observations  typically  exhibit  a  near-normal  distribution 
which  is  bell-shaped  (Baecher  1987,  Benjamin  and  Cornell  1970)  Even 
when  the  material  property  of  interest  does  not  have  a  normal  frequency 
distribution,  the  frequency  distribution  of  sample  means  is  approximately 
normal  (Baecher  1987).  This  idea  is  an  element  of  the  central  limit  theorem 
as  discussed  in  the  previous  chapter. 
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Figure  3-1 .  Characteristic  shapes  of  smooth-curve  approximations  for 
frequency  distributions  {Spiegel  1988) 


Departures  From  Normality 

For  the  case  where  the  frequency  distribution  for  pavement  materials  data 
departs  form  the  “normal”  shape,  the  cause  is  probably  one  of  the  following 
(Rethati  1983): 

a.  Physical  limits. 

b.  Errors  caused  by  sampling  and  testing. 

Physical  limits  means  that  the  possible  numerical  values  for  a  measured 
material  property  have  either  an  upper  limit,  a  lower  limit,  or  both  (Rethati 
1983).  For  example,  unconfined  compressive  strength  for  soil  cannot  be  less 
than  zero,  so  the  results  for  low-strength  soils  would  tend  to  have  a  positive 
skew.  Physical  limits  that  are  not  near  the  mean  of  a  material  property 
would  tend  not  to  induce  skewness.  The  unconfined  compressive  strength  of 
concrete  cannot  be  less  than  zero,  but  average  strengths  are  typically  so  far 
removed  from  zero  that  skewness  is  not  induced  by  this  numerical  boundary. 
An  example  of  potential  negative  skew  occurs  for  measured  degree  of  satu¬ 
ration  for  wet  soil.  If  the  mean  degree  of  saturation  is  relatively  large 
(greater  than  60  percent),  the  distribution  would  tend  to  end  abruptly  at  the 
physical  limit  of  100  percent. 
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To  attain  normality,  the  data  forming  a  frequency  distribution  must  typi¬ 
cally  come  from  a  single  population  (Rethati  1983).  If  the  material  is  con¬ 
crete,  all  the  concrete  must  be  from  the  same  mixture  design.  If  the  material 
is  soil,  horizontal  and  vertical  survey  boundaries  must  be  employed  so  that 
the  soil  can  be  considered  a  single  deposit.  If  two  different  types  of  soil  are 
analyzed  as  a  single  population,  a  sampling  error  would  have  been  commit¬ 
ted  and  the  distribution  for  a  material  property  may  likely  appear  bimodal. 

If  the  property  is  plasticity  index  (PI)  and  the  two  soil  types  are  clays  of  dif¬ 
ferent  plasticity,  the  measurements  of  one  soil  type  may  cluster  around  a  PI 
of  30  (for  example),  while  the  measurements  of  the  second  soil  type  may 
cluster  around  a  PI  of  20.  In  order  to  evaluate  this  data  properly,  the  two 
soil  types  should  be  treated  as  two  different  populations. 

Additional  sampling  errors  can  be  caused  by  deviations  in  the  time  of 
sampling.  Time  can  permit  changes  in  material  properties,  which  can  invali¬ 
date  assumptions  concerning  the  material’s  homogeneity,  including  the 
assumption  that  all  the  samples  come  from  the  same  population  (Rethati 
1983).  If  the  material  is  a  silt  soil  and  the  measured  property  is  moisture 
content,  a  day’s  difference  in  sampling  time  provides  an  opportunity  for 
moisture  content  to  change.  If  the  material  is  an  asphalt  concrete  and  the 
measured  property  is  temperature,  measurements  just  after  loading  a  truck 
should  be  considered  a  separate  population  from  measurements  upon  arrival 
at  the  construction  site. 

Errors  caused  by  testing  can  skew  results  (Rethati  1983).  For  soil,  the 
level  of  disturbance  imposed  on  a  sensitive  clay  during  sampling  will  affect 
shear  strength  measurements.  If  some  of  the  samples  are  severely  disturbed 
they  will  exhibit  lower  strengths,  which  will  cause  the  distribution  to  skew  in 
the  negative  direction.  If  some  concrete  cylinders  are  allowed  to  dry  before 
measuring  compressive  strength,  their  measured  strength  will  be  high,  which 
will  cause  the  distribution  to  skew  in  the  positive  direction. 


Tests  for  Normality 


Many  common  statistical  procedures  assume  that  the  frequency  distribu¬ 
tion  of  data  is  Gaussian,  or  normal.  It  is  important  to  know  when  incorrect 
assumptions  for  normality  can  lead  to  gross  errors  and  to  know  how  to  check 
for  normality.  Incorrect  assumptions  concerning  the  shape  of  probability 
distributions  are  most  critical  when  the  statistical  analysis  is  concerned  with 
the  “tails”  of  the  distribution,  or  the  very  low  and  very  high  values  (Coleman 
and  Steele  1989).  For  example,  the  use  of  statical  methods  for  the  deter¬ 
mination  of  outliers  in  a  data  set  involves  the  tails  specifically.  These 
methods  are  sensitive  to  the  shape  of  the  probability  distribution.  Statistical 
methods  that  use  the  entire  body  of  data,  such  as  methods  for  comparing 
sample  means,  are  less  sensitive  to  the  shape  of  the  probability  distributions. 


The  simplest  method  for  determining  whether  or  not  a  sample  frequency 
distribution  is  normal  involves  visual  inspection  of  the  shape  of  its  histogram. 
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If  the  histogram  is  symmetric  and  bell-shaped,  the  assumption  of  normality  is 
probably  acceptable. 

The  next  level  of  sophistication  for  testing  normality  is  another  graphical 
alternative,  involving  the  production  of  a  “normal  probability  plot.”  Several 
types  of  normal  probability  plots  exist,  but  they  all  share  a  common  char¬ 
acteristic:  a  distribution  that  is  normal  or  nearly  normal  will  plot  as  a 
straight  line,  extending  from  the  lower  left  quadrant  to  the  upper  right  quad¬ 
rant  (Willenbrock  1974c).  Examples  of  these  plots  will  be  shown  as  part  of 
an  example  problem  later  in  this  chapter.  A  popular  style  used  by  statistical 
software  packages  is  to  format  the  x-  and  y-axes  as  unitless  linear  scales. 

The  x-axis  is  used  for  representing  all  real  data  in  terms  of  standard  devia¬ 
tions  from  the  mean.  The  y-axis  is  used  for  representing  standard  normal 
values  (Z),  obtained  using  the  cumulative  percent  frequencies  for  the  data. 

For  example,  if  an  observation  is  greater  than  50  percent  of  the  entire  data 
set,  its  expected  normal  Z  would  be  0.0  (Table  Al).  If  an  observation  is 
greater  than  80  percent  of  the  entire  data  set,  its  expected  normal  Z  would  be 
approximately  0.84.  If  the  data  is  normal,  its  calculated  standard  deviations 
from  the  mean  should  correlate  positively  with  the  standard  normal  devia¬ 
tions  from  the  mean  (Norusis  1993).  These  plots  are  typically  referred  to  as 
Z-plots.  An  example  problem  is  used  later  in  this  chapter  to  demonstrate  the 
construction  of  a  Z-plot. 

Another  version  of  the  normal  probability  plot  uses  linear  scales  for  the  x- 
and  y-axes,  with  percent  as  their  units.  The  x-axis  is  used  to  represent 
cumulative  percent  frequencies  for  the  actual  data.  Each  data  point  is  con¬ 
verted  to  standard  deviations  from  the  mean  and  this  value  is  used  in  the 
standard  normal  Z  table  to  obtain  an  expected  cumulative  percent  frequency, 
which  is  plotted  on  the  y-axis.  If  the  data  is  normal,  its  calculated  cumula¬ 
tive  percentages  should  correlate  positively  with  the  standard  normal  cumula¬ 
tive  percentages  (Norusis  1993).  These  plots  are  typically  referred  to  as 
P-plots.  An  example  problem  is  used  later  in  this  chapter  to  demonstrate  the 
construction  of  a  P-plot. 

Several  methods  exist  for  purely  quantitative  testing  of  normality.  Collec¬ 
tively,  these  methods  are  referred  to  as  “goodness-of-fi”  tests.  In  each 
method,  the  decision  of  normality  is  reduced  to  statistical  testing  of  the  null 
hypothesis,  which  states  that  the  distribution  is  normal  in  shape.  Rejection  or 
acceptance  of  this  hypothesis  involves  a  comparison  between  a  calculated 
statistic  and  a  test  statistic.  These  methods  include  the  chi-square  (jf)  test, 
the  Kolmogorov-Smimov  test,  the  Lilliefors  test,  and  the  Shapiro-Wilks  test. 
These  methods  vary  in  terms  of  their  leniency  towards  rejecting  the  null 
hypothesis  of  normality.  It  is  important  to  realize  that  as  data  sets  become 
large,  each  of  these  goodness-of-fit  tests  becomes  hard  to  satisfy  (Norusis 
1993).  Frequency  distributions  for  real-world  data  will  rarely  exhibit  perfect 
normality,  even  if  they  appear  to  be  symmetric  and  bell-shaped. 

In  order  to  account  for  the  inherent  strictness  of  normality  testing  for 
large  datasets,  the  probability  for  falsely  rejecting  normality  should  be 
decreased,  essentially  making  the  normality  test  more  lenient.  The  prob¬ 
ability  of  falsely  rejecting  normality  is  controlled  by  setting  a  level  of  test 


36 


Chapter  3  Analyzing  Pavement  Materials  Data 


significance  (a).  Practitioners  who  test  for  normality  by  these  quantitative 
methods  should  supplement  their  results  with  normal  probability  plots  This 
practice  will  enable  the  practitioners  to  develop  a  sense  of  the  relationship 
between  dataset  size,  normality,  and  appropriate  level  of  test  significance 


Example.  The  “goodness-of-fit”  test  is  demonstrated  in  this  example 
because  it  is  commonly  used  and  it  is  simple.  The  data  are  moisture  content 
determinations  for  a  residual,  high-plasticity  clay  (Table  3-1).  There  are 
229  measurements,  with  a  mean  of  31.9  percent  and  a  standard  deviation  of 
12.8  percent.  To  begin  a  chi-square  test,  the  data  are  categorized  into  inter¬ 
vals,  as  if  producing  a  frequency  distribution  plot.  Distances  of  the  intervals 
from  the  mean  are  determined  from  the  upper  end  of  the  intervals  and  are 
expressed  in  terms  of  standard  deviations.  This  is  necessary  so  that  the 
actual  data  can  be  compared  to  expected  values  obtained  from  the  standard 
normal  table  (Table  Al).  The  standard  normal  table  provides  the  probability 
of  obtaining  a  sample  value  within  each  class  interval,  if  the  data  was  nor¬ 
mally  distributed.  Expected  observation  frequencies  for  each  class  interval, 
assuming  normality,  can  then  be  obtained  by  multiplying  each  class  interval 
probability  by  the  total  number  of  observations.  The  contribution  to  chi- 
square  from  each  class  interval  is  calculated  as: 


X2  contribution  =  k?^serve^  expected )2 

expected 


The  total  is  then  the  sum  of  all  the  class  interval  contributions.  For  this 
example,  the  total  chi-square  (sum  of  column  F)  was  10.52. 

The  degrees  of  freedom  for  chi-square  is  the  total  number  of  class  inter¬ 
vals  (k)  minus  one,  minus  the  number  of  independent  population  parameters 
estimated  ( m ).  The  subtraction  of  one  from  k,  in  addition  to  subtracting  m, 
accounts  for  a  constraint  condition.  The  constraint  states  that  if  k-1  expected 
frequencies  are  known,  the  remaining  frequency  can  be  determined  (Spiegel 
1990).  For  this  example,  there  were  14  class  intervals  and  we  estimated  two 
population  parameters:  mean  (  =  31.9  percent)  and  standard  deviation 
(  =  12.8  percent).  Therefore,  we  had  11  (14-1-2)  degrees  of  freedom.  From 
Table  A3,  we  can  determine  the  probability  of  obtaining  a  chi-square  larger 
than  10.62  if  the  population  was  normal.  Using  linear  interpolation  we  find 
this  probability  to  be  0.48  or  48  percent.  Therefore  for  this  example  we  can¬ 
not  reject  the  assumption  of  normality.  The  assumption  of  normality  should 
be  rejected  only  if  the  probability  of  obtaining  a  larger  than  the  calculated 
value  is  less  than  the  established  level  of  significance.  For  this  example  if 
the  level  of  significance  had  been  set  at  5  percent,  the  assumption  of  nor¬ 
mality  could  only  have  been  rejected  if  the  probability  of  obtaining  a  r2 
larger  than  the  calculated  value  was  less  than  5  percent. 

Using  the  same  set  of  data,  values  for  both  a  Z-plot  and  a  P-plot  were  cal¬ 
culated,  as  shown  in  Table  3-2.  Values  for  these  plots  can  be  obtained  using 
eitfier  individual  data  points  or  class  intervals.  Statistical  software  packages 
will  typically  use  all  the  data  points.  For  simplicity,  the  class  intervals  were 
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Table  3-1 

Goodness-of-Fit  for  Moisture  Contents  (Percent)  of  Residual, 
Hiah-Plasticity  Clay  (adapted  from  Steel  and  Torrie  1980) 

Class  Interval 

A 

B 

C 

D 

E 

F 

0.6  to  5.5 

7 

-26.43 

-2.065 

0.0194 

4.44 

1.48 

5.6  to  10.5 

5 

-21.43 

-1.674 

0.0277 

6.34 

0.28 

10.6  to  15.5 

7 

-16.43 

-1.284 

0.0525 

12.0 

2.08 

15.6  to  20.5 

18 

-11.43 

-0.893 

0.0863 

19.8 

0.16 

20.6  to  25.5 

32 

-6.43 

-0.502 

0.1219 

27.9 

0.60 

25.6  to  30.5 

41 

-1.43 

-0.112 

0.1477 

33.8 

1.53 

30.6  to  35.5 

37 

3.57 

0.279 

0.1545 

35.4 

0.07 

35.6  to  40.5 

25 

8.57 

0.670 

0.1386 

31.7 

1.42 

40.6  to  45.5 

22 

13.57 

1.060 

0.1068 

24.5 

0.26 

45.6  to  50.5 

19 

18.57 

1.451 

0.0712 

16.3 

0.45 

50.6  to  55.5 

6 

23.57 

1.841 

0.0405 

9.27 

1.15 

55.6  to  60.5 

6 

28.57 

2.232 

0.0201 

4.60 

0.43 

60.5  to  65.5 

3 

33.57 

2.623 

0.0084 

1.92 

0.61 

II  65.6  to  70.5 

1 

38.57 

3.010 

0.0044 

1.01 

0.00 

A  =  Observed  frequency. 

B  =  Deviation  of  high  endpoint  from  mean. 

C  =  Standard  deviations  from  mean. 

D  =  Normal  probability  of  obtaining  a  value  in  the  class  interval  (Table  A1 ). 

E  =  Expected  frequency,  based  on  normality. 

F  -  Contribution  to  chi-square  [(observed-expected)2/expected). 

used  in  this  example  problem.  The  plots  are  shown  in  Figures  3-2  and  3-3. 
Based  on  the  results  from  the  x2  analysis,  one  would  expect  the  plots  to  show 
a  strong  positive  correlation  between  the  data  and  the  standard  normal 
values.  The  plots  support  this  expectation  and  show  that  the  data  are  at  least 
nearly  normal. 


Transformations 

If  a  data  set  is  found  not  to  be  normally  distributed,  two  options  exist  for 
transforming  the  data  into  a  new  data  set  that  may  be  normally  distributed. 
The  first  option  is  to  obtain  or  calculate  new  data  that  have  similar  engineer¬ 
ing  value.  For  example,  California  bearing  ratio  measurements  for  in-situ 
clay  may  not  be  normally  distributed,  but  plate-bearing  load  test  results  for 
the  same  clay  may  be  normal  in  shape. 
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Table  3-2 

u"PL°Diand-  P’Plot  Calcu,ati°ns  for  Moisture  Contents  of  Residual, 


|  Class  Interval 

Observed 

Frequency 

Z-Plot 

b  ume  ijjo  \J)  ] 

P-Plot 

x-axis1 

y-axis2 

x-axis3 

y-axis4 

0.6  to  5.5 

7 

-2.065 

-1.866 

3.1 

1  9 

5.6  to  10.5 

5 

-1.674 

-1.626 

5.2 

4  7 

10.6  to  15.5 

7 

-1.284 

-1.385 

8.3 

10.0 

15.6  to  20.5 

18 

-0.893 

-0.986 

16.2 

18  6 

20.6  to  25.5 

32 

-0.502 

-0.519 

30.1 

30  8 

25.6  to  30.5 

41 

-0.112 

-0.050 

48.0 

45.6 

30.6  to  35.5 

37 

0.279 

0.364 

64.2 

61  0 

35.6  to  40.5 

25 

0.670 

0.678 

75.1 

74.9 

40.6  to  45.5 

22 

1.060 

1.024 

84.7 

85  5 

45.6  to  50.5 

19 

1.451 

1.476 

93.0 

92.7 

50.6  to  55.5 

6 

1.841 

1.706 

95.6 

96  7 

55.6  to  60.5 

6 

2.232 

2.120 

98.3 

98  7 

60.5  to  65.5 

3 

2.623 

2.650 

99.6 

99  6 

65.6  to  70.5 

1 

3.010 

N/A 

100 

99.9 

■  ===^ J 

2  !n'arS®  'r:™!®  ^Presented  as  standard  deviation  from  the  mean  (real  data), 
frequency  dafd  deV,at,0n  from  the  mean  <based  on  real  cumulative  percent 
|  Cumulative  percent  frequency  (real  data). 

meNa°n)mal  CUmU'ative  frec>uency  aa  a  Percent  (based  on  real  standard  deviation  from  the 

N/A  Not  applicable. 

The  second  option  for  transforming  data  into  a  normal  distribution 
involves  recomputing  all  data  with  either  a  square  root,  a  logarithm,  or  a 
trigonometric  function.  If  a  data  set  is  lognormal  in  shape,  taking  the  loga¬ 
rithm  of  all  data  points  should  create  a  normal  distribution.  Conversely  tak¬ 
ing  the  logarithm  of  normally  distributed  data  can  cause  the  distribution  to 
deviate  from  normality.  While  soil  porosities  (p)  and  void  ratios  (e)  are 

?Q^n°rmalIy  distributed>  lof(P)  and  log(e)  are  typically  not  (Schultze 
1972).  As  an  example  of  a  trigonometric  transformation,  researchers  found 
that  the  angle  of  internal  friction  (<j>)  for  sand  could  be  transformed  into  a 
normal  distribution  by  applying  the  cotangent  function  (Schultze  1972). 

Example.  Suppose  an  airfield  runway  needs  to  be  evaluated  for  its  allow¬ 
able  passes  by  a  fully-loaded  F-15  aircraft.  Suppose  further  that  119  esti¬ 
mates  for  allowable  passes  have  been  obtained  from  non-destructive 
evaluations  using  a  falling-weight  deflectometer.  An  example  data  set  orga¬ 
nized  by  class  intervals,  is  shown  in  Table  3-3.  These  data  were  evaluated 
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Observed  Standard  Deviations  from  Mean 


Figure  3-2.  Z-Plot  for  a  example  problem  concerning  the  moisture 
content  of  soil 


for  normality  using  both  the  test  and  the  Kolmogorov-Smirnov  test.  The 

test  did  not  reject  normality,  but  the  Kolmogorov-Smirnov  test  rejected 
normality  at  a  significance  level  of  5  percent.  These  data,  which  are  shown 
in  Figure  3-6,  appear  to  have  a  positive  skew.  This  was  confirmed  by  calcu¬ 
lating  their  coefficient  of  skew  (=  0.62). 

Distributions  for  allowable  passes  are  often  lognormal  in  shape,  so  trans¬ 
formations  using  the  logarithmic  function  can  often  improve  normality.  The 
transformed  allowable  passes  data  are  also  shown  in  Table  3-3.  Neither  the 
£  test  or  the  Kolmogorov-Smirnov  test  rejected  normality  for  these  trans¬ 
formed  data.  These  data,  which  are  shown  in  Figure  3-7,  do  not  appear  to 
have  any  significant  skew.  This  was  confirmed  by  calculating  their  coeffi¬ 
cient  of  skew  (=  -0.07). 
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Observed  Cumulative  Percent  Frequency 


Figure  3-3.  P-Plot  for  an  example  problem  concerning  the  moisture 
content  of  soil 


Handling  Potential  Outliers 


In  handling  large  quantities  of  data,  some  results  may  appear  to  be 
isolated  and  much  different  than  the  rest  of  the  population.  If  the  physical 
reason  for  the  odd  results  is  known,  these  data  may  be  (Rethati  1983): 

a.  Discarded. 


b.  Discarded  and  replaced. 


c.  Corrected  on  a  physical  basis. 

Corrections  on  a  physical  basis  involve  the  use  of  physical  laws  or  known 
physical  relationships  to  correct  sampling  and  testing  errors. 
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Table  3-3  1 

Transformation  of  Allowable  Passes  Data  for  an  Airfield  Feature  | 

Allowable  Passes 

Log  (Allowable  Passes)  |j 

Class  Intervals 

Observations 

Class  Intervals 

Observations 

181000-200000 

3 

5.26-5.28 

1 

201000-220000 

2 

5.29-5.31 

2 

221000-240000 

5 

5.32-5.34 

2 

241000-260000 

6 

5.35-5.37 

4 

261000-280000 

9 

5.38-5.40 

5 

281000-300000 

14 

5.41-5.43 

7 

301000-320000 

16 

5.44-5.46 

12 

321000-340000 

21 

5.47-5.49 

16 

341000-360000 

10 

5.50-5.52 

21 

361000-380000 

11 

5.53-5.55 

16 

381000-400000 

5 

5.56-5.58 

11 

401000-420000 

5 

5.59-6.61 

8 

421000-440000 

4 

5.62-6.64 

6 

441000-460000 

2 

5.65-5.67 

3 

461000-480000 

2 

5.68-5.70 

3 

481000-500000 

1 

5.71-5.73 

1 

501000-520000 

1 

5.74-5.76 

1 

521000-540000 

1 

541000-560000 

1 

Example.  Suppose  a  concrete  ready-mix  plant  performs  quality  control 
testing  with  152x305-mm  cylinders.  Suppose  further  that  a  few  data  points 
in  a  month’s  production  appear  to  be  outliers  (or  extreme  values),  well  above 
the  mean.  The  concrete  supplier  may  discover  that  these  data  were  obtained 
with  76xl52-mm  cylinders.  The  supplier  may  choose  to  factor  down  the 
strength  measurements  with  known  relationships  between  the  strengths  mea¬ 
sured  by  the  two  cylinder  sizes.  Whether  discarding  or  correcting,  all 
activities  should  be  documented. 

If  the  physical  reason  for  outliers  is  not  known,  the  outliers  may  only  be 
removed  from  the  data  set  with  proper  statistical  procedures.  Even  with 
these  tools,  data  removal  should  be  performed  sparingly  and  should  be  docu¬ 
mented.  Removing  data  can  be  counterproductive  if  the  true  dispersion  or 
shape  of  a  distribution  is  lost.  If  the  apparent  outliers  are  actually 
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Figure  3-4.  Variance  components  for  a  construction  site  study 


reflections  of  a  data  set’s  skewness, 
provide  a  false  image  of  symmetry. 


then  the  removal  of  the  outliers  would 


When  considering  removal  of  outliers  by  statistical  methods,  all  possible 
physical  reasons  for  the  outlying  data  should  be  investigated.  The  sampling 
procedure  should  be  reviewed  to  ensure  that  the  outlying  data  came  from  the 
same  population  as  the  other  data.  Test  procedures  should  be  reviewed  to 
ensure  that  the  outlying  data  was  tested  in  the  same  manner  as  the  rest  of  the 
observations.  Finally,  data  analyses  should  be  reviewed  to  ensure  that  there 
were  no  calculation  errors.  Once  these  efforts  have  been  exhausted,  the  next 
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Figure  3-7.  Frequency  distribution  for  transformed  (log)  allowable  passes 


h?JS  f  WhCther  **  distribution  ^  normal  or  other  shape.  The 

f  f  Can  mfluence  ^  choices  for  oulier  testing.  Published 

™  a„rd  rhermmm?  °UtherS  mClUde  **  T~statistic’  ^  Tietjen-Moore 
£2 f  ’ r  2T  ?0t  S  cnten0n-  The  T-statistic,  which  is  calculated  dif- 
ferendy  than  Student  s  t-statistic  discussed  in  Chapter  2,  can  only  be  applied 

ormal  populations  and  can  only  handle  a  single  outlier.  The  Tietien- 
Moore  statistic  can  also  only  be  applied  to  normal  populations  but  it  can  test 
or  multiple  low  and  high  outliers.  Chauvenot’s  criterion  can  detect  multiple 
outliers  and  it  can  be  applied  to  distributions  that  are  other  than  normal. 

nn Jht6  rStJ1StiC’  f°r  USe  With  a  single  outlier-  squires  conversion  of  the 
potential  outlier  to  a  test  criterion  (Tcalculated): 


calculated 


yt-y 


(51) 
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where 


yi  -  the  potential  outlier 
y  =  the  sample  mean 
s  =  standard  deviation  for  the  sample. 

This  test  criterion,  Tcalculaled,  is  compared  to  standard  T-values  (Table  A5)  to 
determine  the  probability  of  obtaining  a  deviation  of  this  magnitude,  relative 
to  the  sample  mean.  When  using  the  standard  T-values,  the  engineer  would 
select  a  level  of  significance,  which  will  define  the  probability  of  falsely  dis¬ 
carding  a  data  point  as  an  outlier.  According  to  guidelines  in  ASTM  (1995b) 
E  178,  levels  of  significance  of  5  percent  or  less  are  generally  advisable  for 
testing  outliers.  If  TcalculaIed  is  larger  than  the  standard  T-value  at  the  selected 
level  of  significance,  the  potential  outlier  may  be  discarded  (Schiff  and 
D’Agostino  1996). 

Example.  Assume  that  the  following  data  represents  extracted  asphalt 
cement  contents  from  a  single  lot:  5.68,  5.70,  5.70,  5.70,  5.72,  5.72,  5.72, 
5.78,  5.84,  and  5.96.  The  mean  and  standard  deviation  for  this  data  is  5.75 
percent  and  0.09  percent,  respectively.  If  we  question  the  validity  of  the 
measured  value  5.96,  we  would  calculate  T  as  (5.96-5.75)/0.09  =  2.33. 
Inspection  of  Table  A5  reveals  that  the  probability  of  occurrence  for  a 
T-value  this  high  is  approximately  2  percent.  If  we  had  chosen  a  significance 
level  of  5  percent,  we  would  conclude  that  the  extreme  value  did  not  come 
from  the  same  population  as  the  rest  of  the  data  and  it  would  be  discarded. 

The  Tietjen-Moore  (T-M)  statistic  can  be  applied  in  cases  where  multiple 
outliers  exist  on  either  the  low  side  of  the  data,  the  high  side  of  the  data,  or 
on  both  sides  of  the  data  (Tietjen  and  Moore  1972).  Let  the  sample  values 
be  yh  y2,  y3,  ...,  yn  and  compute  the  sample  mean,  y.  Then  compute  the  n 
absolute  residuals: 


/-,  =  lyj-yl.  '2  =  fri-yl  ••••  =  ly.-yl  (52> 

Now  relabel  the  original  observations,  y„  y2,  y3 . as  z’s  in  such  a 

manner  that  z,  is  that  observation  whose  absolute  residual  is  the  i,h  largest. 
Relabeled  observation  Z/  should  have  the  smallest  residual  and  relabeled 
observation  z„  should  have  the  largest  residual.  The  Tietjen-Moore  statistic 
for  testing  the  significance  of  the  k  largest  residuals  is  then: 

n-k 

£  (zrz*)2 

Ek  =  ^ -  (53) 

*  n 

£(z,-z)2 
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where 


n-k 


h  =  £ 

1=1 


(n-k) 


(54) 


The  term  zk  represents  the  mean  of  the  (n-k)  least  extreme  observations 
and  the  term  z  is  the  mean  of  the  full  sample.  If  the  calculated  value  of  E„  is 
smaller  than  the  critical  value  of  Ek,  shown  in  Table  A6  (Appendix  A)  then 
the  outlier  can  be  discarded.  As  with  the  T-statistic,  the  engineer  will ’have 
to  select  a  level  of  significance  (a),  which  will  define  the  probability  of 
falsely  discarding  a  potential  outlier. 


For  the  example  data  shown  in  Table  3-4,  a  T-M  statistic  can  be  calcu¬ 
lated  for  each  of  three  different  cases:  one  outlier  (E,),  two  outliers  (£,) 

and  three  outliers  (E3).  The  calculated  T-M  statistic  for  each  case  is  shown 
m  Table  3-5. 


Table  3-4 

Percent  of  Particles  Finer  Than  the  No.  200  Sieve  for  Asphalt 

1  Concrete  Aqqregate 

Observation  (Percent) 

Absolute  Residual  (r.) 

(2,Z)2 

6.1 

0.21 

0.0441 

6.0 

0.11 

0.0121 

5.8 

0.09 

0.0081 

5.7 

0.19 

0.0361 

5.9 

0.01 

0.0001 

5.5 

0.39 

0.1521 

5.4 

0.49 

0.2401 

5.6 

0.29 

0.0841 

6.0 

0.11 

0.0121 

6.9 

1.01 

1.0201 

Average  observation  (z)  =  5.890  percent. 

£  U-,-z)2  =  1 .609  percent2. 

The  critical  T-M  statistics  for  this  problem  can  be  obtained  from 
Table  A6.  Assuming  that  the  chosen  level  of  significance  is  5  percent,  the 
T-M  statistic  for  10  observations  for  one,  two,  and  three  outliers  is  0.353 
0.172,  and  0.083,  respectively.  An  observation  can  be  considered  an  outlier 
if  the  calculated  T-M  statistic  is  less  than  the  critical  T-M  statistic.  There¬ 
fore,  the  data  set  has  only  one  outlier.  This  outlier  would  be  the  observation 
with  the  largest  absolute  residual,  which  is  6.9  percent. 
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Table  3-5 

Execution  of  the  Tietjen-Moore  Statistic  for  Determining  of 

Outliers 

j  Number  of  Outliers 

Ek 

1 

5.778 

0.476/1.609  =  0.296 

2 

5.825 

0.315/1.609  =  0.196 

1  3 

5.871 

0.194/1.609  =  0.121 

Chauvenof s  criterion  was  originally  developed  for  applications  related  to 
the  standard  normal  distribution  (Coleman  and  Steel  1989).  However,  as 
described  below,  it  can  be  applied  to  any  known  or  assumed  probability  den¬ 
sity  function.  The  criterion  permits  rejection  of  all  points  that  fall  outside  of 
a  designated  band  around  the  sample  mean.  The  width  of  the  band  cor¬ 
responds  to  an  area  under  the  probability  density  function  equal  to  l-(l/(2n)) 
percent,  where  n  equals  the  number  of  measurements  in  the  data  set.  For 
sample  sizes  ranging  from  5  to  50,  the  area  for  the  acceptance  region  ranges 
from  0.90  to  0.99,  corresponding  to  areas  for  rejection  regions  from  0.10  to 
0.01,  respectively,  as  shown  in  Table  3-6.  Chauvenot’s  criterion  is  always 
treated  as  a  two-tailed  significance  test:  one-half  of  the  rejection  region  lies 
at  the  upper  end  of  the  probability  density  function  and  one-half  of  the  rejec¬ 
tion  region  lies  at  the  lower  end  of  the  probability  density  function. 

If  the  distribution  of  data  is  assumed  to  be  normal,  potential  outliers  can 
be  compared  to  the  rejection  band  limits  by  converting  their  measured  values 
to  a  Chauvenot  Criterion  as  follows: 


x 


calculated 


(55) 


where 


y,  =  the  potential  outlier 

y  =  the  sample  mean 

s  =  standard  deviation  for  the  sample 

This  conversion  is  the  same  as  that  used  for  the  T-statistic.  The  values  for 
Tcalculmed  can  then  be  compared  to  values  for  vrejection,  which  represent  the 
boundaries  of  the  two-tail  rejection  region  shown  in  Table  3-5.  Maintaining 
the  assumption  of  normality,  the  values  for  'crejection  can  be  obtained  from  any 
standard  normal  table,  such  as  that  shown  in  Table  Al.  Each  value  for 
Rejection  is  obtained  using  the  appropriate  two-tailed  rejection  region. 

Using  the  10  observations  from  the  T-statistic  example,  determination  of 
r  would  require  the  use  of  0.025  as  the  two-tailed  rejection  region,  as 

r€J€CtlQJ I  T 
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Table  3-6 
Chauvenot's  ( 
Steele  1989) 

Criteria  for  Determination  of  Outliers  (Coleman  and 

Number  of 

Samples 

Area  of 

Acceptance 

Total  Area  of 
Rejection 

One-Half  Rejection 

Area1 

5 

0.90 

0.100 

0.050 

10 

0.95 

0.050 

0.025 

20 

0.975 

0.025 

0.0125 

50 

r 

0.99 

0.010 

0.005 

1  To  be  used  for  two-tailed  tests 

defined  in  Table  3-6.  From  Table  Al,  rrejeaion  is  found  to  equal  1.96.  The 
Va  ^ calculated  for  the  highest  measured  asphalt  cement  content  of  5  96 

would  be  2.33,  which  is  greater  than  trejeclion.  This  measurement  happens  to 
be  the  only  data  point  with  vcalculated  greater  than  vreJection.  Similar  to  the 
T-statistic  conclusion,  Chauvenot’s  criterion  indicates  that  5.96  should  be 
discarded  as  an  outlier. 


Using  the  data  from  the  Tietjen-Moore  (T-M)  example,  which  also 
included  10  observations,  rreJeclion  would  be  defined  similar  to  the  T-statistic 
example  (1.96).  Using  the  sample  standard  deviation  of  0.42  percent,  the 
lowest  and  highest  measurements  for  percent  passing  would  transform  to 
values  of  1.17  and  2.40,  respectively.  The  highest  measurement 

haS  \  TTT  Val.Ue  larger  than  of  L96-  This  mea¬ 

surement  happens  to  be  the  only  data  point  with  r  greater  than  r 

Therefore  similar  to  the  T-M  statistic  conclusion"  Chauvenot’s  criterion""0" 
indicates  that  the  measurement  of  6.9  percent  passing  the  No.  200  sieve 
should  be  discarded  as  an  outlier. 


If  the  distribution  is  known  to  conform  to  a  beta  shape,  rather  than  nor¬ 
mal,  the  concept  of  Chauvenot’s  criterion  can  still  be  applied  However 
instead  of  comparing  the  data  to  a  standard  normal  distribution,  the  data  is 
compared  to  a  theoretical  beta  distribution.  The  area  of  the  acceptance 
region  can  still  be  calculated  as  l-l/(2n)  percent  and  the  area  of  the  rejection 
region  as  l/(2n)  percent.  The  rejection  region  can  still  be  considered  two- 
tailed.  The  band  limits  on  the  distribution,  to  be  used  for  rejection  criteria 
can  then  be  determined  using  Harr’s  program  BETA  (Harr  1987).  The  user 
of  BETA  only  needs  to  know  the  distribution's  minimum  value,  maximum 
value,  mean,  and  coefficient  of  variation.  If  the  rejection  area  has  been 
determined  to  be  0.10,  then  the  upper  band  limit  corresponds  to  the  value 
that  has  only  a  5  percent  chance  of  being  exceeded  by  randomly  selected 
values.  The  lower  band  limit  corresponds  to  the  value  that  has  a  95  percent 
chance  of  being  exceeded  by  randomly  selected  values. 

Regardless  of  the  distribution  shape  or  the  method  of  testing  for  outliers 
the  test  should  only  be  performed  one  time  for  a  set  of  data.  In  other  words, 
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once  one  or  more  points  have  been  rejected,  the  tests  for  outliers  should  not 
be  repeated  on  the  new,  smaller  data  set.  Repeated  application  of  outlier 
tests  on  the  same  data  set  would  result  in  an  erroneous  overall  significance 
level  (a). 


Components  of  Variability 

When  studying  the  variability  of  materials,  there  are  two  types  of  analyses 
that  are  of  primary  interest  to  the  construction  industry.  The  first  type  of 
analysis  is  intended  to  study  the  variability  of  a  test  method  or  process  speci¬ 
fically.  In  this  case,  the  material  to  be  subjected  to  the  test  is  collected  and 
homogenized  at  a  single  location.  Typically,  the  material  is  then  sent  to 
several  laboratories  for  testing  by  a  single  operator  from  each  laboratory. 
These  analyses  are  often  referred  to  as  “round  robin”  studies  because  sam¬ 
ples  from  the  same  homogeneous  material  are  provided  to  all  the  parties 
involved. 

The  second  type  of  anlysis  is  intended  to  differentiate  between  the  vari¬ 
ability  associated  with  sampling  and  testing  a  heterogeneous  material 
obtained  from  a  construction  site,  rock  quarry,  asphalt  plant,  etc.  In  order  to 
keep  the  components  of  variability  reasonable  in  this  type  of  analysis,  a 
single  laboratory  should  be  involved  in  the  testing.  These  analyses  are  often 
referred  to  as  “construction  site”  studies  because  the  results  pertain  to  a 
particular  project  and  the  results  are  useful  for  developing  reasonable  con¬ 
tract  specification  criteria. 


Round  robin  studies 

Round  robin  studies  are  often  used  for  developing  precision  statements  for 
test  methods.  In  accordance  with  ASTM  (1995a)  C  670,  two  basic  elements 
of  precision  are  needed:  single-operator  precision  and  multilaboratory  preci¬ 
sion.  Single  operator  precision  provides  an  estimate  of  the  difference  that 
may  be  expected  between  duplicate  measurements  made  on  the  same  material 
in  the  same  laboratory  by  the  same  operator  using  the  same  apparatus  within 
a  time  span  of  a  few  days  [ASTM  (1995a)  C  802].  This  within-laboratory 
precision  is  often  referred  to  as  the  repeatability  of  the  test.  Multilaboratory 
precision  provides  an  estimate  of  the  difference  that  may  be  expected 
between  measurements  made  on  the  same  material  in  two  different  labora¬ 
tories  [ASTM  (1995a)  C  802].  This  between-laboratory  precision  is  often 
referred  to  as  the  reproducibility  of  the  test  method. 

At  least  ten  different  laboratories  should  be  included  in  a  round  robin 
study  [ASTM  (1995a)  C  802],  The  number  of  different  types  of  a  material 
that  should  be  included  in  this  type  of  study  depend  on  many  factors,  such  as 
the  range  of  material  types  to  which  the  test  method  should  apply,  the 
expense  of  distributing  multiple  samples,  and  the  test  duration.  In  general,  a 
minimum  of  three  materials  should  be  included  [ASTM  (1995a)  C  802].  The 
number  of  replicates  necessary  for  each  test  on  each  material  is  a  function  of 
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the  number  of  laboratories  involved.  If  15  laboratories  are  participating 
then  three  replicates  is  adequate.  If  less  than  15  laboratories  are  involved, 
the  number  of  replicates  should  be  increased  according  to:  (30/p)  + 1  where 
p  =  the  number  of  laboratories.  If  more  than  15  laboratories  are  involved 
the  number  of  replicates  can  be  reduced  to  two  [ASTM  (1995a)  C  802], 

Example.  Assume  we  are  interested  in  the  repeatability  and  reproduci¬ 
bility  associated  with  determining  plasticity  index  (PI)  for  soil.  The  fabri¬ 
cated  data  for  this  study,  which  included  15  laboratories,  3  types  of  soil  and 
3  test  replicates,  are  shown  in  Table  3-7.  Within-laboratory  averages  and 
withm-laboratory  variances  are  shown  in  Table  3-8.  The  overall  mean  and 
the  pooled  variance  for  PI  for  each  soil  type  are  also  shown  in  Table  3-8. 


Table  3-7 

Raw  Data  for  a  Fabricated  Round  Robin  Experiment  Concerning 

Plasticity  Index  of  Soils  a 

1  ~  "  '  ■' - - 

Laboratory 

High-Plasticity  Soil 

Medium-Plasticity  Soil 

Low-Plasticity  Soil 

1 

30 

31 

30 

13 

11 

11 

4 

3 

3 

2 

32 

31 

30 

10 

11 

12 

8 

5 

7 

3 

35 

35 

34 

14 

13 

14 

5 

6 

5 

4 

30 

30 

31 

7 

6 

5 

2 

4 

2 

5 

36 

35 

37 

1 1 

12 

10 

5 

3 

4 

6 

36 

35 

36 

10 

11 

11 

1 

2 

3 

7 

35 

36 

37 

13 

11 

12 

1 

2 

2 

8 

31 

32 

32 

12 

13 

14 

6 

5 

7 

9 

36 

37 

38 

11 

10 

1 1 

8 

7 

6 

10 

31 

30 

30 

8 

8 

9 

4 

5 

4 

11 

31 

30 

32 

1 1 

10 

12 

5 

7 

8 

12 

31 

32 

32 

12 

14 

13 

5 

6 

7 

13 

30 

31 

30 

7 

6 

5 

2 

4 

2 

14 

32 

33 

33 

1 1 

13 

11 

3 

4 

3 

15 

35 

36 

36 

11 

10 

11 

2 

1 

3 

Overall 

variance 

6.6 

5.5 

4.1 

Before  proceeding  with  a  statistical  analysis  of  the  round  robin  data,  we 
need  to  verify  two  assumptions:  homogeneity  of  within-laboratory  variances 
and  absence  of  substantial  interactions  between  laboratory  and  soil  type.  If  a 
single  laboratory  has  an  erratic  variance  compared  to  the  others,  it  may  be 
eliminated.  If  all  the  variances  are  erratic,  however,  there  is  a  problem  with 
the  test  method.  Substantial  interactions  between  laboratories  and  materials 
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Table  3-8 

Summary  of  Mean  and  Within-Laboratory  Variance 


Hiqh-Plasticity  Soil 

Medium-Plasticity  Soil 

Low-Pla 

sticity  Soil 

Laboratory 

Mean 

Variance 

Mean 

Variance 

Mean 

Variance 

1 

30.3 

0.33 

11.7 

1.33 

3.3 

0.33 

2 

31.0 

1.00 

11.0 

1.00 

6.7 

2.33 

3 

34.7 

0.33 

13.7 

0.33 

5.3 

0.33 

4 

30.3 

0.33 

6.0 

1.00 

2.7 

1.33 

5 

36.0 

1.00 

11.0 

1.00 

4.0 

1.00 

6 

35.7 

0.33 

10.7 

0.33 

2.0 

1.00 

7 

36.0 

1.00 

12.0 

1.00 

1.7 

0.33 

8 

31.7 

0.33 

13.0 

1.00 

6.0 

1.00 

9 

37.0 

1.00 

10.7 

0.33 

7.0 

1.00 

10 

30.3 

0.33 

8.3 

0.33 

4.3 

0.33 

11 

— 

31.0 

1.00 

11.0 

1.00 

6.7 

2.33 

12 

31.7 

0.33 

13.0 

1.00 

6.0 

1.00 

13 

30.3 

0.33 

6.0 

1.00 

2.7 

1.33 

14 

32.7 

0.33 

11.7 

1.33 

3.3 

0.33 

15 

35.7 

0.33 

10.7 

0.33 

2.0 

1.00 

Pooled12 

33.0 

0.55 

10.7 

0.82 

4.2 

1.00 

1  Overall  mean  [  =  (£within-lab  means)/number  of  labs]. 

2  Pooled  within-laboratory  variance  [  =  (V within-lab  variances)/number  of  labs]. 


would  also  indicate  either  a  problem  with  a  laboratory  or  a  problem  with  the 
test  method.  In  an  extreme  case  of  interaction,  the  laboratories  may  not 
report  the  same  ranking  of  materials,  in  terms  of  the  magnitude  of  test 
results. 

Homogeneity  of  variances  for  each  material  may  be  tested  with  the 
Hartley  F-max  test  (Freund  and  Wilson  1993).  This  method  simply  involves 
finding  the  largest  ratio  of  within-laboratory  variances  and  then  comparing 
this  ratio  to  the  F-max  distribution  data  shown  in  Table  A7  (Appendix  A). 
The  null  hypothesis  for  this  test  states  that  the  variances  are  homogeneous. 
The  alternative  hypothesis  states  that  at  least  two  variances  are  not  equal.  If 
the  calculated  ratio  is  larger  than  the  critical  value  shown  in  Table  A7,  then 
the  null  hypothesis  would  be  rejected.  The  largest  ratio  of  variances  found 
within  high-plasticity  soil,  medium-plasticity  soil,  and  low-plasticity  soil  was 
3.0  (1.0/0.33),  4.0  (1.33/0.33),  and  7.0  (2.33/0.33),  respectively.  For 
15  laboratories,  3  replicates  per  laboratory,  and  a  level  of  significance  of 
5  percent,  the  critical  F-max  value  is  855.  The  critical  F-max  value  is  high 
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because  our  estimates  for  within-laboratory  variability  are  poor  (only  3  repli¬ 
cates).  We  have  no  reason  to  reject  the  null  hypothesis  of  equal  variances 


Substantial  interactions  can  be  checked  by  inspecting  the  reported  means 
for  plasticity  index  in  Table  3-8.  If  we  ranked  the  reported  means  from  each 
laboratory,  we  would  find  that  each  laboratory  reported  the  lowest  PI  for  the 
same  material  and  the  highest  PI  for  the  same  material.  There  does  not  seem 
to  be  any  laboratory  that  reported  a  different  trend  in  results,  relative  to  the 
other  laboratories.  We  can  state  that  there  are  no  substantial  interactions. 


Having  verified  the  assumptions  of  homogeneous  variances  and  the 
absence  of  substantial  interaction,  a  statistical  computer  program  such  as 
bigmaStat  or  SAS  can  be  used  to  perform  an  analysis  of  variance  (ANOVA) 
procedure  for  each  soil  type.  This  permits  estimation  of  the  between- 
laboratory  variances  and  verification  of  the  pooled  within-laboratory  vari¬ 
ances.  Since  the  number  of  replicates  was  the  same  for  each  laboratory  the 
pooled  within-laboratory  variance  in  this  case  is  calculated  as  a  simple  a’ver- 
age.  The  analysis  is  similar  to  a  completely  randomized  experimental  design 
where  laboratories  are  the  treatment  factors  and  the  variation  within  labora 
tones  reflects  experimental  error  and  sampling  error.  Summary  ANOVA 

output  tables  are  shown  in  Tables  3-9  through  3-11. 


Table  3-9 

Round  Robin  Analysis  of  Variance  for  Plasticity  Index  Using  Hiqh- 
Plasticity  Soil  a 

-  - 

Source  of  Variability 

Degrees  of 
Freedom1 

Sum  of 

Squares 

Mean 

Square 

Expected 

Mean  Square2 

Laboratory 

II 

Q. 

273.2 

19.52 

o  2  +  r(OH2) 

Error 

P(r-1)  =  30 

16.7 

0.55 

oj 

Corrected  total 

pr-1  =  44 

289.9 

N/A 

N/A 

’  P  ~  Number  of  labs,  r  =  number  of  tests  per  lab. 

°b  ,=  Variance  between  labs,  o„2  =  variance  within  labs. 

N/A  No  calculation. 

The  ANOVA  tables  list  two  individual  sources  of  variability:  error  and 
laboratory.  The  error  component  represents  within-laboratory  variability. 
The  laboratory  component  includes  both  within-laboratory  variability  and 
between-laboratory  variability.  The  ANOVA  results  include  the  calculation 
of  a  sum  of  squares  term  and  a  mean  square  term  for  each  source  of  vari¬ 
ability  as  shown.  A  sum  of  squares  is  simply  the  sum  of  squared  deviations 
from  the  mean,  as  discussed  in  Chapter  2.  A  mean  square  is  the  sum  of 
squares  divided  by  the  degrees  of  freedom  for  the  particular  source  of  vari¬ 
ability  The  mean  square  terms  can  be  compared  to  the  expected  mean 
square  formulas,  which  are  known  because  they  comprise  the  mathematical 
basis  of  the  ANOVA.  Collectively,  these  comparisons  provide  two  equations 
and  two  unknowns  (within-laboratory  variance  and  between-laboratory 
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variance).  After  calculating  these  variances,  they  can  be  added  to  estimate 
the  overall  variance.  These  computed  results  are  summarized  in  Table  3-12. 


Table  3-10 

Round  Robin  Analysis  of  Variance  for  Plasticity  Index  Using  J 

Medium-Plasticity  Soil  1 

Source  of  Variability 

Degrees  of 
Freedom1 

Sum  of  Squares 

Mean  Square 

Expected 

Mean  Square2 

Laboratory 

p-1  =  14 

219.0 

15.64 

aj  +  r(oh2) 

Error 

p(r-1 )  =  30 

24.7 

0.82 

_2b? _ 

Corrected  total 

pr-1  =  44 

243.6 

N/A 

N/A 

1  p  =  Number  of  labs,  r  =  number  of  tests  per  lab. 

2  ob2  =  Variance  between  labs,  ow2  =  variance  within  labs. 

N/A  No  calculation. 

Table  3-1 1 

Round  Robin  Analysis  of  Variance  for  Plasticity  Index  Using  Low- 
Plasticity  Soil 

Degrees  of 

Sum  of 

Expected  Mean 

Source  of  Variability 

Freedom1 

Squares 

Mean  Square 

Square2 

l 

p-1  =  14 

150.3 

10.74 

o„2  +  r(o„2) 

Error 

p(r-1 )  =  30 

30.0 

1.00 

oj 

Corrected  total 

pr-1  =  44 

180.3 

||  2  0h2  _  Variance  between  labs,  ow2  = 

variance  within  labs. 

Table  3-12 

Summary  of  Round  Robin  PI  Data  by  Soil  Type 


. 7  w 

Soil  Type 

Mean, 

Percent 

Components  of  Variance 

Within- 

Laboratory 

Between-Laboratory 

Overall 

Variance1 

Low-plasticity 

4.2 

1.00 

3.25 

4.25 

Medium-plasticity 

10.7 

0.82 

4.94 

5.76 

Hiah-Dlasticity 

33.0 

0.55 

6.32 

6.87  1 

1  Overall  variance  -  within-laboratory  variance  +  between-laboratory  variance. 
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A  comparison  between  Tables  3-8  and  3-12  reveals  that  the  ANOVA- 
derived  withm-laboratory  variances  match  the  hand-calculated  pooled 
variances  exactly.  Inspection  of  Table  3-12  also  reveals  that  the  between- 
aboratory  variances  were  approximately  three  to  eleven  times  larger  than  the 
within-laboratory  variances.  The  estimated  variances  and  the  calculated 
means  can  be  used  to  calculate  standard  deviations  and  coefficients  of  vari¬ 
ation,  as  shown  in  Table  3-13.  The  soil  types  have  been  organized  in  the 
table  in  order  of  increasing  mean  to  permit  a  search  for  trends  in  standard 
deviation  and  coefficient  of  variation.  These  trends  will  determine  the  most 
appropriate  form  for  the  precision  statements. 


Table  3-13 

Standard  Deviations  and  Coefficients  of  Variation  for  PI  Data 

Soil  Type 

Mean, 

Percent 

Standard  Deviation, 
Percent 

Coefficient  of  Variation, 
Percent 

Within- 

Laboratorv 

Between- 

Laboratorv 

Within- 

Laboratorv 

Between- 

Laboratory 

Low-plasticity 

4.2 

1.00 

1.80 

24 

43 

Medium-plasticitv 

10.7 

0.91 

2.22 

8.5 

21 

High-plasticity 

33.0 

0.74 

2.51 

2.2 

7.6 

The  form  of  precision  statements  will  depend  on  the  relationship  between 
standard  deviation  (within-laboratory  and  between-laboratory)  and  mean 
material  property.  Generally,  this  relationship  will  fall  into  one  of  two  cate¬ 
gories:  either  the  standard  deviation  remains  approximately  constant  over  the 
range  of  mean  values  or  the  coefficient  of  variation  remains  approximately 
constant  over  the  range  of  mean  values.  In  the  case  of  constant  standard 
deviation,  the  coefficient  of  variation  would  have  to  decrease  as  mean 
increases.  In  the  case  of  constant  coefficient  of  variation,  the  standard  devi¬ 
ation  would  have  to  increase  as  mean  increases.  Inspection  of  Table  3-13 
reveals  that  for  this  round  robin  experiment,  standard  deviation  remained 
approximately  constant  across  the  range  of  mean  PI,  for  both  within- 
laboratory  and  between-laboratory  components  of  variability.  The  precision 
statement  would  therefore  be  most  useful  if  written  in  terms  of  standard  devi¬ 
ation,  rather  than  in  terms  of  coefficient  of  variation. 

The  precision  statement  developed  from  this  fabricated  round  robin 
experiment  is  shown  in  Table  3-14.  The  standard  deviations  reported  in 
Table  3-14  were  calculated  by  averaging  over  the  three  soil  types  included  in 
the  study.  Assuming  that  the  test  method  specification  will  require  two  repli¬ 
cates,  the  precision  statement  will  include  an  allowable  difference  for  two 
measurements,  which  is  calculated  by  multiplying  the  standard  deviation 
by  2.8.  According  to  ASTM  (1995a)  C  670,  this  difference  between  two 
measurements  represents  that  which  should  be  exceeded,  on  the  average 
only  5  percent  of  the  time.  There  are  many  other  forms  of  precision  state¬ 
ments,  including  those  that  are  based  on  coefficient  of  variation.  Also, 
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Table  3-14 

Precision  Statement  for  Plasticity  Index 

Type  Index 

Standard  Deviation,  Percent 

Acceptable  Range  of  Two 
Test  Results,  Percent 

Sinqle-operator 

0.88 

2.51 

Multilaboratory 

2.18 

6. 12  1 

1  Each  test  result  is  one  c 

2  Each  test  result  is  a  sine 

f  two  replicates  from  a  single  laboratory. 

}le  replicate  and  each  replicate  is  from  a  different  laboratory. 

precision  statements  can  be  written  for  specifications  that  require  more  than 
two  replicates.  Detailed  information  on  developing  precision  statements  is 
provided  in  ASTM  (1995a)  C  670  and  ASTM  (1995b)  E  177. 


Construction  site  studies 

The  variabilities  of  concern  for  construction  site  operations  typically 
involve  a  single  testing  laboratory.  Multiple  laboratories  would  not  normally 
be  involved  in  a  single  project’s  quality  control  or  quality  assurance.  Typi¬ 
cally,  a  contractors’  performance  and  subsequent  payment  are  judged  for 
each  “lot”  of  material  separately.  A  lot  can  be  defined  as  an  isolated  quan¬ 
tity  of  material  or  process  accumulated  under  conditions  that  are  considered 
uniform  for  sampling  purposes  [ASTM  (1995b)  E  105,  AASHTO  (1995) 

R  10].  Examples  of  lots  include  1,800  metric  tons  of  asphalt  concrete,  8  hrs 
worth  of  production  of  asphalt  concrete,  a  single  lift  of  select  fill  considered 
to  be  constructed  with  a  uniform  material  and  under  uniform  conditions,  or  a 
1 ,000-ft  section  of  completed  airfield  taxiway. 

The  concepts  presented  in  the  following  discussion  are  applicable  for  sam¬ 
pling  and  testing  all  construction  materials,  but  portland  cement  concrete  will 
be  used  as  the  example.  A  single  lot  of  concrete  can  consist  of  many  truck- 
loads,  each  of  which  can  be  considered  a  “batch.”  The  overall  variation  of 
concrete  in  a  lot  can  therefore  be  broken  down  into  between-batch  variation 
and  within-batch  variation.  Between-batch  variation  represents  the  difference 
in  test  results  from  one  batch  to  another  for  a  material  that  is  supposedly  pro¬ 
duced  consistently. 

The  within-batch  variation,  which  represents  the  variabilities  that  exist 
within  any  single  batch,  can  be  broken  down  into  sampling  error  and  experi¬ 
mental  error.  Sampling  error  is  caused  by  some  inconsistency  in  the  sam¬ 
pling  procedure  that  disturbs  the  random  selection  process.  A  popular 
example  is  sampling  aggregate  from  a  segregated  stockpile.  Experimental 
error  can  be  further  broken  down  into  two  components:  testing  error  and 
inherent  variability.  Testing  errors  are  caused  by  inconsistencies  that  occur 
while  performing  the  same  test  on  multiple  samples.  Testing  variation  is 
affected  by  the  repeatability  of  a  test,  technician  skills,  the  condition  of 
equipment,  and  the  methods  used  for  reducing  samples  to  a  size  that  is 
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appropriate  for  testing  (Nicotera  1974).  Testing  variation  can  be  quantified 
y  repeatmg  the  same  test  on  similar  samples  of  material.  Inherent  vari¬ 
ability  is  governed  by  the  laws  of  chance  and  it  is  unavoidable.  Its  magni¬ 
tude  will  vary  depending  on  the  material  and  the  test.  For  a  particular 
±e  influence  of  inherent  variability  on  an  analysis  of  variance 

wN.?Vmfn  0nly  b£  reduced  by  increasing  the  number  of  test  replicates 
(Waller  1966). 


Between-batch  variation  typically  exceeds  within-batch  variation.  This 
relationship  is  similar  to  the  case  where  between-laboratory  variation  typi¬ 
cally  exceeds  within-laboratory  variation  for  round  robin  experiments 
(Nicotera  1974). 

The  components  of  variance  for  a  construction  site  material  study  are 
additive.  The  Pythagorean  theorem  provides  a  visual  tool  for  relating  vari¬ 
ance  components  back  to  standard  deviations  (Newlon  1966,  Waller  1966). 

If  the  length  of  each  orthogonal  side  of  a  right  triangle  is  equal  to  a  compo¬ 
nent  standard  deviation,  the  variance  associated  with  the  hypotenuse  equals 
the  sum  of  the  orthogonal  sides  squared,  as  shown  in  Figure  3-4.  Overall 
variance  equals  between-batch  variance  plus  within-batch  variance.  Within- 
batch  variance  equals  sampling  variance  plus  variance  associated  with  experi¬ 
mental  error.  Experimental  variance  is  composed  of  testing  variability  and 
inherent  material  variability,  but  these  components  are  difficult  to  separate 
(Nicotera  1974). 


V2  =  c2 

overall  a 


between-batch 


+  S 


within-batch 


(56) 


within-batch 


=  S 2 


sampling 


+  j" 


experimental 


(57) 


Between-batch  variance  reflects  the  “control”  of  a  construction  operation 
so  it  is  typically  the  quantity  needed  forjudging  contractor  performance.  If’ 
the  batch-to-batch  variance  is  large  relative  to  the  within-batch  variance  the 
operation  can  be  judged  “out  of  control”  (Waller  1966).  Accurate  determi¬ 
nation  of  material  variances  is  facilitated  by  two  procedures:  (1)  reducing 
sampling  and  testing  variabilities  through  standardization  of  sampling  and 
testing  procedures  and  (2)  increasing  the  number  of  test  replicates  within  a 
sampling  plan  (Newlon  1966). 


Similar  to  the  round  robin  experiment,  different  components  of  variance 
can  be  isolated  with  an  anlaysis  of  variance  (ANOVA)  procedure.  The  pro- 
cedure  used  for  the  ANOVA  will  depend  on  the  number  and  nature  of  the 
different  sources  of  variation. 

Example.  Brown  (1966)  used  an  analysis  of  variance  procedure  to  isolate 
sources  of  variation  for  the  compressive  strengths  of  steam  cured  concrete 
He  was  concerned  with  the  relative  proportions  of  variation  that  existed 
between  concrete  batches  and  within  concrete  batches:  a2I0[al  =  a2  + 

°  within-  A  large  variation  between  groups,  relative  to  the  variation  within 
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groups,  would  indicate  the  existence  of  a  material  control  problem.  Example 
data  is  shown  in  Table  3-15. 


Table  3-15 

Compressive  Strengths  of  Steam-Cured  Concrete,  MPa  (after 

Brown  1966) 

1  Batch  No. 

Sample  No.  1 

Sample  No.  2 

Total 

1 

30.3 

30.1 

60.4 

2 

31.2 

31.4 

62.6 

3 

28.1 

29.0 

57.1 

4 

31.7 

32.3 

64.0 

5 

31.2 

31.7 

62.9 

Overall  total  (for  N  = 

10  samples) 

307.0 

———————— 

Overall  average 

30.7 

Due  to  the  few  sources  of  variation  and  the  small  dataset,  Brown  (1966) 
was  able  to  demonstrate  the  ANOVA  with  hand  calculations,  as  shown  in 
Table  3-16.  This  analysis  is  similar  to  a  completely  randomized  experi¬ 
mental  design  where  concrete  batches  are  the  treatment  factors  and  variation 
within  batches  reflects  experimental  error  and  sampling.  The  between-batch 
component  of  variance  accounted  for  a  significant  proportion  (93  percent)  of 
the  overall  variance. 

If  this  ANOVA  had  been  performed  with  commercial  software,  such  as 
SAS®,  the  output  would  look  similar  to  Table  3-17.  Similar  to  the  ANOVA 
in  the  round  robin  experiment,  the  calculated  mean  squares  can  be  compared 
to  the  expected  mean  squares  to  obtain  estimates  for  two  sources  of  variabil¬ 
ity.  In  this  case,  the  two  sources  are  within-batch  variance  and  between- 
batch  variance.  Results  from  calculations  using  the  ANOVA  output  are 
summarized  in  Table  3-18.  Similar  to  the  hand  calculations,  the  between- 
batch  variance  accounted  for  92  percent  of  the  total  variance.  The  large  ratio 
of  between-batch  variance  to  within-batch  variance  indicated  to  Brown  (1966) 
that  he  had  a  concrete  production  control  problem. 

Example.  In  a  study  with  a  similar  experimental  design,  Hughes  and 
Anday  (1966)  used  analysis  of  variance  techniques  to  optimize  the  efficiency 
of  nuclear  density  gage  quality  control  testing.  In  order  to  optimize  effi¬ 
ciency,  they  needed  to  differentiate  between  two  components  of  variance 
within  a  typical  lot  of  soil:  between-test  variability  and  between-test-site 
variability.  Between-test  variability  referred  to  variance  between  repeated 
tests  performed  at  the  same  location.  Between-test-site  variability  referred  to 
variance  between  the  test  results  at  different  locations  within  the  same  lot. 
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Table  3-16 

Example  of  Computational  Procedures  for  Variance  Analysis 
(after  Brown  1966) 

Computational  Steps 

-  ===========-= J 

Between- 

Batches 

Within- 

Batches 

Total 

a.  Enter  squared  values 

60. 42  +  62. 62  + 
...  +  62. 92 

N/A 

30. 32  +  30.1 2  + 

...  +  31. 72 

b.  Sum  of  "a" 

18,879.7 

N/A 

9,440.6 

c.  No.  of  tests  per  each  entry 
in  "a" 

2 

N/A 

1 

d.  Crude  sum  of  squares, 

9439.9 

N/A 

9440.6 

e.  Correction  factor, 

C  =  (qrand  total)2/N 

9424.9 

N/A 

9424.9 

f.  Corrected  sum  of  squares 

(d,*-e) 

15.0 

-  f,) 

0.7 

(d3  -  e) 

15.7 

g.  Degrees  of  freedom 

(#  batch  -  1) 

4 

(93  ‘  9,) 

5 

(N-  1) 

9 

h.  Mean  square,  "f"/"q” 

3.75 

0.14 

N/A 

1.  Component  of  variance 

(h,  -  h2)/c, 

1.81 

h2 

0.14 

(■l  +  <2> 

1.95 

j.  Component  of  variance,  % 

93 

7 

100 

— - — '  ===s=====sas=s=ss=========;===j 

N/A  No  calculation. 

Subscripts  represent  column  numbers  under  "Sources  of  Variation  * 

Table  3-17 

|  Analysis  of  Variance  for  Concrete  Compressive  Strenath 

Source  of  Variability 

Degrees  of 
Freedom1 

Sum  of 

Squares 

Mean 

Square 

Expected  Mean 
Square2 

Batches 

P-1  =  4 

14.97 

3.74 

oJ  +  r(oK2) 

Error 

p(r-1 )  =  5 

0.75 

0.15 

oJ 

Corrected  total 

1"" 

pr-1  =  9 

15.72 

N/A 

N/A 

l  p  ~  Number  of  batches,  r  =  number  of  tests  per  batch. 

m/a6  J  Va,ria?ce  between  batches,  ow2  =  variance  within  batches. 

N/A  No  calculation. 

They  selected  five  sections  of  roadway,  each  with  a  different  type  of  soil 
Mid  they  analyzed  the  variance  components  for  each  section  separately. 
Eight  sites  for  testing  were  selected  randomly  within  each  section  and  two 
tests  (nuclear  gage  determinations)  were  performed  at  each  site.  The 
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Table  3-18 

Summary  of  ANOVA  Calculation  Results  for  Concrete 
Compressive  Strength 

Com 

ponents  of  Variance 

Between- 

Overall 

Soil  Type 

Within-Laboratory 

Laboratory 

Variance1 

Component  of  variance  (MPa2) 

0.15 

1.80 

1.95 

Percent  of  variance 

8 

92 

100 

H  '  Overall  variance  -  within-laboratory  variance  +  between-laboratory  variance. 

between-test-site  variation  provided  an  indication  of  the  variabilities  associ¬ 
ated  with  the  material  and  the  construction  process.  The  between-test  vari¬ 
ation  provided  an  indication  of  the  repeatability  of  the  nuclear  gage  apparatus 
when  used  on  the  same  material.  Hughes  and  Anday  (1966)  found  that  the 
variation  between  tests  accounted  for  a  very  small  proportion  of  the  total 
variation.  The  variation  between  test-site  means,  however,  was  statistically 
significant.  They  concluded  that  if  sixteen  tests  were  to  be  performed  within 
a  section  (or  lot),  the  standard  error  for  section  means  would  be  minimized 
by  selecting  sixteen  different  sites  and  running  only  one  test  per  site.  Statis¬ 
ticians  refer  to  this  as  maximizing  experimental  precision  (Petersen  1985). 

Hughes  and  Anday  (1966)  were  not  concerned  with  variance  caused  by 
sampling  procedures  because  the  repeated  tests  were  performed  in  the  same 
access  hole.  In  the  concrete  production  example,  Brown  (1966)  also  did  not 
feel  it  necessary  to  differentiate  between  the  components  of  within-batch 
variation.  In  some  cases,  however,  an  engineer  may  want  to  break  down 
within-batch  variation  into  sampling  error  and  experimental  error.  These 
types  of  analyses  become  slightly  more  complex. 

In  order  to  differentiate  between  sampling  variance  and  experimental  vari¬ 
ance,  one  needs  a  sampling  plan  such  as  that  shown  in  Figure  3-5.  Multiple 
batches  are  randomly  selected  from  a  lot  of  material  and  then  at  least  two 
samples  are  randomly  selected  from  each  batch.  Finally,  material  for  at  least 
two  individual  test  replicates  are  split  out  of  each  sample.  The  number  of 
batches  selected  from  each  lot  would  be  decided  as  part  of  a  construction 
specification  acceptance  plan.  This  decision  is  based  on  a  number  of  factors, 
including  predetermined  party  risks,  anticipated  material  variabilities,  and 
type  of  payment  adjustment  scheme.  It  is  important  to  “split  out  the  test 
replicates  from  the  samples.  The  analysis  must  be  able  to  assume  that  the 
test  replicates  are  similar  representations  of  the  same  sample. 

The  ANOVA  compares  the  results  of  individual  test  portions  to  define  an 
experimental  (testing  and  inherent)  variation.  By  mathematically  combining 
the  test  portions  within  each  sample,  the  ANOVA  can  compare  samples  to 
estimate  the  sampling  variation.  Each  set  of  sample  increments  is  then 
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mathematically  combined  into  an  average  batch  value  to  permit  determination 
or  batch-to-batch  variation  (Nicotera  1974). 

Example.  Data  presented  by  Nicotera  (1974)  will  be  used  to  demonstrate 
this  type  of  anlysis  of  variance.  The  data  in  Table  3-19  are  values  for  per¬ 
cent  passing  the  19-mm  sieve  for  aggregates  obtained  from  asphalt  concrete 
extractions.  Two  samples  (A  and  B)  were  obtained  from  40  different 
batches.  Each  sample  had  enough  material  for  two  extraction/gradation  tests. 

Hand  calculations  for  this  level  of  ANOVA  would  be  tedious.  Results 
provided  by  SAS®  software,  however,  are  shown  in  Table  3-20.  The  experi¬ 
mental  design  can  be  considered  as  a  completely  randomized  treatment 
design  with  nested  subsampling.  The  batches  are  the  treatment  factors  and 
the  samples  are  nested  within  batches.  The  variability  within  samples  repre¬ 
sents  experimental  error.  Similar  to  the  previous  examples,  components  of 
variance  can  be  estimated  using  the  expected  mean  square  formulas  and  the 
calculated  mean  squares  shown  in  Table  3-20.  Results  for  these  calculations 
are  shown  in  Table  3-21.  Inspection  of  the  components  of  variance  reveals 
that  variation  between  batches  accounted  for  36  percent  of  the  overall  vari¬ 
ance.  Testing  and  inherent  variation  accounted  for  45  percent  of  the  overall 
variance,  which  is  a  relatively  large  proportion.  Sampling  accounted  for 
only  19  percent  of  the  overall  variance.  At  first  glance,  the  high  proportion 
of  variance  within  batches  (45  +  19=64  percent)  would  indicate  that  the  pro¬ 
cedure  for  measuring  percent  passing  the  19-mm  sieve  was  not  very  consis¬ 
tent.  However,  upon  more  thorough  study,  this  high  proportion  of  the 
variance  is  attributable  to  a  low  overall  variance.  The  overall  coefficient  of 
variation  for  percent  passing  the  19-mm  sieve  was  only  2.7  percent  The 
control  of  coarse  aggregates  for  the  production  of  asphalt  concrete  appears  to 
be  m  control.  The  repeatability  of  asphalt  extraction/gradation  tests  appears 
to  be  favorable. 

As  a  method  for  checking  the  ANOVA  results,  the  overall  variance  can  be 
calculated  by  spreadsheet  and  compared  with  the  estimated  value  shown  in 
Table  3-21.  Overall  variance  is  simple  to  calculate  because  it  includes  all  the 
data  (160  values),  without  regard  for  batch,  sample,  or  test  designations. 
Overall  variance  in  this  case  was  calculated  to  be  6.08  percent2  which 
agreed  with  the  ANOVA  results. 
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4  Historical  Variability  Data 


As  part  of  this  study,  a  literature  review  was  performed  to  collect  vari¬ 
ability  data  related  to  pavement  materials  and  pavement  structures.  These 
data  are  necessary  for  the  proper  development  of  specification  criteria, 
including  pay  adjustment  factors.  Details  concerning  the  data  found  in  liter¬ 
ature,  including  references,  are  presented  in  Appendices  B  through  F.  This 
chapter  serves  as  a  summary  for  the  appendices. 

The  material  variabilities  reported  in  this  chapter  involved  a  single  labora¬ 
tory.  Variabilities  associated  with  multiple  laboratories  would  typically  not 
be  an  issue  for  project  quality  testing.  Since  contractor  pay  is  typically  cal¬ 
culated  for  individual  lots  on  a  paving  project,  the  variabilities  reported  in 
this  chapter  are  intended  to  reflect  total  within-lot  variabilities.  Total  within- 
lot  variabilities  include  between-batch  variabilities,  between-sample  variabil¬ 
ities,  and  experimental  variation  (testing  and  inherent). 

When  several  estimates  for  variability  were  provided  by  the  same  source, 
variances  were  “pooled”  into  a  single  estimate  as  described  in  Chapter  2.  If 
the  number  of  replicates  for  variance  estimates  were  not  given,  pooling  of 
variances  was  performed  as  a  straight  average  (replicates  were  assumed  to  all 
be  equal). 

Material  variabilities  are  summarized  in  this  chapter  with  both  reported 
standard  deviations  (SDs)  and  coefficients  of  variation  (CVs).  Reported 
variabilities  for  any  particular  material  property  followed  one  of  three  trends 
with  respect  to  the  property  mean: 

a.  As  the  mean  increased,  the  standard  deviation  remained  approximately 
constant  while  the  coefficient  of  variation  decreased. 

b.  As  the  mean  increased,  the  coefficient  of  variation  remained 
approximately  constant  while  the  standard  deviation  increased. 

c.  The  range  of  reported  means  was  too  small  to  decipher  these  trends. 

If  SD  remained  approximately  constant,  SD  is  summarized  as  a  single  value 
and  CV  is  summarized  as  a  range.  If  CV  remained  approximately  constant, 
CV  is  summarized  as  a  single  value  and  SD  is  summarized  as  a  range.  If  the 
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range  of  reported  means  was  too  small  to  decipher  these  trends  SD  and  CV 
are  both  summarized  as  single  values. 


Any  comments  found  in  literature  concerning  the  normality  of  data  are 
included  m  this  chapter.  In  most  cases,  these  comments  are  just  perceptions 
related  to  the  shape  of  frequency  distributions.  However,  in  some  cases  a 
normality  test  was  performed.  All  reports  for  normality  address  individual 

test  results.  This  chapter  does  not  address  any  distributions  for  sample 
means.  F 


Residual  Fine-Grained  Soil  Deposits 


Krahn  and  Fredlund  (1983)  reported  on  several  soil  deposits  in  Canada. 
They  noted  a  surprising  consistency  with  respect  to  the  variability  of  individ¬ 
ual  soil  properties  between  deposits.  This  finding  supports  the  logic  behind 
collecting  and  summarizing  variability  data  for  use  in  future  construction 
projects. 

The  measured  properties  within  a  “lot”  (uniform  deposit)  of  residual  soil 
typically  exhibit  coefficients  of  variation  in  the  range  of  10  to  30  percent 
(Table  4-1).  Measured  properties  with  relatively  low  variability  include 
specific  gravity  of  soil  particles  and  dry  density.  Measured  properties  with 
relatively  high  variabilities  include  permeability,  plasticity  index,  unconfined 
compressive  strength,  and  shear  box  strength  results.  When  soil  strength 
tests  require  field  sampling  and  transportation  of  these  samples  to  a  labora¬ 
tory,  the  variability  of  measured  strength  is  highly  dependent  on  sampling 
technique  (Kay  and  Krizek  1972). 

Some  of  the  measured  properties  were  reported  to  be  non-normal:  natural 
water  content  (percent),  degree  of  saturation  (percent),  maximum  dry  density 
(kg/m  ),  and  unconfined  compressive  strength  (kPa).  Two  of  these  proper¬ 
ties  can  most  likely  attribute  their  deviations  from  normality  to  limiting 
numerical  boundaries.  Natural  water  contents  (w)  cannot  fall  below  0  per¬ 
cent,  so  the  distribution  for  a  soil  with  a  low  mean  w  would  tend  to  exhibit  a 
positive  skew.  Degree  of  saturation  (S)  cannot  exceed  100  percent  so  the 
distribution  for  a  soil  with  a  high  mean  S  would  tend  to  exhibit  a  negative 
skew.  Unconfined  compressive  strength  can  most  likely  attribute  its  devi¬ 
ation  from  normality  to  heterogeneity  induced  by  sampling.  Sampling  would 

have  the  potential  to  decrease  measured  strength,  causing  the  distribution  to 
exhibit  a  negative  skew. 


Engineered  Fill 


The  measured  properties  within  a  “lot”  of  engineered  fill  typically  exhibit 
coefficients  of  variation  in  the  range  of  10  to  30  percent  (Tables  4-2  and 
4-3).  Relative  compaction  had  a  low  variability  (approximately  5  percent), 
possibly  because  it  is  often  used  for  payment.  Contractors  typically  have  ’ 
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I  Table  4-1 

Summary  of  Reported  Variabilities  for  Residual  Fine-Grained  Soil 

Deposits 

I 

Standard 

Coefficient  of 

Property 

Deviation 

Variation,  Percent 

Normality 

Specific  qravity 

N/A 

2 

normal’ 

Natural  water  content, 
percent 

4 

10-25 

conflicts2 

Porosity 

0.05 

10-20 

normal1 

Void  ratio 

0.15 

20 

normal’ 

Degree  of 
saturation,  percent 

8 

7-15 

not  normal’ 

Dry  density,  kg/m3 

105 

5 

not  reported3 

Coefficient  of 

not  reported3 

permeability,  percent 
-  80  percent  saturated 

N/A 

90 

-  1 00  percent  saturated 

N/A 

700 

Liquid  limit 

4-15 

20 

normal’ 

Plastic  limit 

2-5 

15 

not  reported3  j 

Plasticity  index 

3-10 

40 

not  reported3 

Particles  finer,  percent 
-  No.  40  sieve 

15 

15 

not  reported3 

-  No.  200  sieve 

10 

20 

-  2  A/m  sieve 

10 

35 

Compression  index 

0.030- 

0.070 

30 

not  reported3 

Swell  index 

0.020 

25-55 

not  reported3 

Optimum  moisture 

Content,  percent 

5 

20 

normal4 

Maximum  dry  density, 
kq/m3 

130 

8 

not  normal4 

|  (Continued) 

1  According  to  the  chi-square 

"goodness-of-fit"  test. 

1  2  Reports  concerning  normality  conflicted. 

I  3  Reports  on  variability  did  not  address  normality. 

1  4  Based  on  the  shapes  of  histograms. 

U  N/A  Limited  data  or  not  reported. 
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Table  4-1  (Concluded) 

" - — -  - 

Property 

Standard 

Deviation 

Coefficient  of 
Variation,  Percent 

Normality 

Unconfined  compression 
strength,  kPa 

7-190 

45 

not  normal1 

California  bearing  ratio, 
percent 

3 

25 

not  reported3 

Field  vane,  kPa 

8 

30 

not  reported3 

Undrained  triax.  Test 

-  c 

-  tan(f) 

N/A 

N/A 

25 

20 

not  reported3 

Drained  triaxial  test 

-  c 

-  tanij) 

N/A 

N/A 

15 

2 

not  reported3 

Drained  shear  box 

-  c 

-  tanc]) 

N/A 

N/A 

100 

30 

not  reported3 

j  According  to  the  chi-square  "goodness-of-fit"  test. 

Reports  on  variability  did  not  address  normality. 

N/A  Limited  data  or  not  reported 

Table  4-2 

Summary  of  Reported  Variabilities  for  Compacted  Fine-Grained 

Soil  Deposits 

Property 

Standard 

Deviation 

Coefficient  of 
Variation,  Percent 

Normality 

Controlled  moisture 
content,  percent 

2 

20 

conflicts1 

Relative 

compaction,  percent 

4 

5 

conflicts1 

Resilient 
modulus,  MPa 

25 

25 

not  reported2 

California  bearing  ratio, 
percent 

1-15 

30 

not  reported2 

Modulus  of  subgrade 
reaction,  MPa/mm 

N/A 

35 

not  reported2 

Benkelman  beam 
deflections,  mm 

0.35 

25 

not  reported2 

Dynaflect  tests,  MPa 

N/A 

35 

not  reported2 

1  Reports  concerning  normality  conflicted. 

Reports  on  variability  did  not  address  normality. 

N/A  Limited  data  or  not  reported 
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Table  4-3 

Summary  of  Reported  Variabilities  for  Lime-Soil  Mixtures  and 
Soil-Cement 

Standard 

Coefficient  of 

Property 

Deviation 

Variation,  Percent 

Normality 

Unconfined  compression 
strenqth,  Mpa 

N/A 

10 

not  reported1 

Splitting-tensile  strength, 

MPa 

N/A 

10 

not  reported1 

Flexural  strength,  MPa 

N/A 

10 

not  reported1 

1  Reports  on  variability  did  not  address  normality. 

N/A  Limited  data  or  not  reported. 

financial  interest  in  meeting  density  requirements  and  this  objective  is  facili¬ 
tated  by  keeping  density  variabilities  as  low  as  possible.  The  CV  for  both 
controlled  moisture  content  and  California  bearing  ratio  are  in  the  same  range 
as  the  corresponding  values  for  residual  soil.  The  CV  for  unconfined  com¬ 
pressive  strength,  however,  was  lower  for  engineered  fill  than  for  residual 
soil. 

Inspection  of  Tables  4-2  and  4-3  also  reveals  that  conflicts  existed  for 
some  of  the  reports  on  normality.  These  conflicts  were  present  for  both  con¬ 
trolled  moisture  content  (percent)  and  relative  compaction  (percent).  The 
tendency  for  these  distributions  to  be  skewed  may  be  affected  by  mean  value 
and  variability.  Controlled  moisture  content  would  tend  to  exhibit  a  normal 
distribution  when  the  mean  is  far  removed  from  the  numerical  limit  of  zero. 
Relative  compaction  would  tend  to  exhibit  a  normal  distribution  when  the 
variability  is  low.  A  high  variability  for  relative  compaction  with  a  mean 
value  close  to  100  percent  would  tend  to  induce  a  negative  skew. 


Subbase  and  Base  Course  Materials 


The  liquid  limit  and  plasticity  index  for  untreated  subbase  and  base  course 
materials  have  SDs  that  are  approximately  the  same  as  for  the  fine-grained 
soils.  The  lower  mean  values  for  these  properties  for  subbase  and  base 
courses,  however,  cause  the  CVs  to  be  higher  than  for  the  fine-grained  soil 
(Table  4-4).  Similar  to  fine-grained  soil,  permeability  measurements  for  sub¬ 
base  and  base  course  materials  are  highly  variable.  Standard  deviations 
(SDs)  for  percent  passing  individual  sieve  sizes  range  from  approximately 
1  percent  to  5  percent,  while  the  CVs  for  percent  passing  individual  sieve 
sizes  range  from  approximately  5  percent  to  50  percent.  Measurement  with 
fine  sieves  tend  to  have  the  highest  CVs  because  they  have  the  lowest  mean 
values  for  percent  passing.  Finer  sieve  sizes  also  tend  to  be  susceptible  to 
errors  induced  by  sampling  method,  screen  cleanliness,  and  skill  of  the 
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Table  4-4 

Summary  of  Reported  Variabilities  for  Subbase  and  Base  Course 
Materials 

Property 

Standard 

Coefficient  of 

Deviation 

Variation,  Percent 

Normality 

Controlled  moisture  content, 
percent 

1 

10 

not  reported' 

Liquid  limit,  percent 

5 

10-80 

not  reported' 

Plasticity  index,  percent 

2 

25-65 

not  normal2 

Coefficient  of  permeability, 
percent 
-  silty  sand 

N/A 

240 

not  reported' 

Particles  finer,  percent 
-  1 9-mm  sieve 

3 

3-10 

not  reported' 

-  1 3.2-mm  sieve 

2 

3 

-  9.5-mm  sieve 

5 

3-20 

-  No.  4  sieve 

5 

4-20 

-  No.  8  sieve 

3 

5-15 

-  No.  1 6  sieve 

4 

20-25 

-  No.  30  sieve 

3 

10 

-  No.  50  sieve 

3 

5-25 

-  No.  100  sieve 

2 

30-45 

-  No.  200  sieve 

1 

5-25 

Sand  equivalence,  percent 

5 

15 

not  reported' 

Magnesium  sulfate 
soundness,  percent 

0.20-1.5 

3 

not  reported1 

Lift  thickness,  mm 
-  subbase 

8-55 

35 

not  reported' 

-  base 

5-30 

15 

Relative 

compaction,  percent 
-  subbase 

3 

3 

not  reported' 

-  base  ; 

2 

2 

Angle  of  internal  friction,  <t> 

N/A 

10 

not  reported' 

Unconfined  compression 
strength,  kPa 

N/A 

10 

not  reported1 

California  bearing  ratio 
(field),  percent 

8-35 

25 

not  reported' 

Benkelman  beam 
deflections,  mm 

0.25 

20 

not  reported' 

1  Reports  on  variability  did  not  address  normality. 

Based  on  the  shapes  of  histograms. 

N/A  Limited  data  or  not  reported. 
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technician  (Michigan  State  Department  of  Highways  1966).  Similar  to  the 
fine-grained  soils,  relative  compaction  had  a  low  variability. 

Tests  for  mechanical  properties  for  subbase  and  base  course  materials 
tend  to  have  high  variabilities,  as  shown  in  Tables  4-5  and  4-6.  Laboratory 
tests  for  modulus,  Poisson’s  ratio,  and  fatigue  properties  had  particularly 
high  variabilities.  The  plate-load  test,  which  is  performed  in  the  field,  also 
had  a  particularly  high  variability. 


Table  4-5 

Summary  of  Reported  Variabilities  for  Lime-  and  C 
Stabilized  Subbase  and  Base  Course  Materials 

ement- 

============— 

Property 

Standard 

Deviation 

Coefficient  of 
Variation,  Percent 

Normality 

Unconfined  compression 
strenqth,  kPa 

840 

15 

not  reported’ 

Compression  modulus,  MPa 

2,400 

60 

not  reported’ 

Indirect  tensile  strength,  kPa 

340 

35 

normal2 

Tensile  modulus,  MPa 

420 

70 

not  reported’ 

California  bearing  ratio 
(field),  percent 

90 

30 

not  reported’ 

Plate-load  tests,  MPa/mm 

N/A 

70 

not  reported' 

Dvnafiect  tests,  MPa 

N/A 

20 

not  reported’ 

1  Reports  on  variability  did  not  address  normality. 

2  Based  on  the  shapes  of  histograms. 

N/A  Limited  data  or  not  reported. 

Table  4-6 

Summary  of  Reported  Variabilities  for  Asphalt-Stabilized  Base 
Course  Materials 

Property 

Standard 

Deviation 

Coefficient  of 
Variation,  Percent 

Normality 

Indirect  tensile  strength,  kPa 

150 

20 

normal’ 

Tensile  modulus,  MPa 

140 

35 

not  reported2 

Indirect  poisson’s  ratio 

0.15 

50 

not  reported2 

Flexural  stiffness,  Mpa 

130-1,090 

25 

not  reported2 

Fatigue  life  of  beams 

N/A 

50 

not  reported2 

'  Based  on  the  shapes  of  histograms. 

2  Reports  on  variability  did  not  address  normality. 

N/A  Limited  data  or  not  reported. 

70 


Chapter  4  Historical  Variability  Data 


Asphalt  Concrete 


Penetration  measurements  for  asphalt  cement  had  lower  CVs  than 
capillary  viscosity  measurements,  as  shown  in  Table  4-7. 

Standard  deviations  for  percent  passing  individual  sieve  sizes  for  asphalt 
concrete  aggregates,  as  shown  in  Table  4-8,  were  very  similar  to  those  for 
subbase  and  base  course  materials.  The  CVs  for  bulk  specific  gravity  maxi¬ 
mum  theoretical  gravity,  and  field  density  were  all  relatively  low.  The  CV 
for  voids  total  mix  was  high  due  to  the  small  mean  values  (typically  3  to  6). 
Similar  to  stabilized  subbase  and  base  course  materials,  measurements  of 
modulus,  Poisson  s  ratio,  and  fatigue  had  high  variabilities. 

P°r/ield  tests’  deflections  measured  by  a  falling-weight  deflectometer 
(  WD)  were  generally  less  variable  than  deflections  measured  by  a  Benkel- 
man  beam  (Table  4-9).  Backcalculated  moduli  for  FWD  were  least  variable 
ror  the  subgrade  and  most  variable  for  asphalt  concrete. 


Table  4-7 

1 

L  Summary  of  Reported  Variabilities  for  Asphalt  Cement  II 

Property 

Standard 

Deviation 

Coefficient  of 
Variation,  Percent 

Normality 

Penetration  at  25°C, 

0.01  mm 

2-20 

25 

— _l 

not  reported'  j 

Viscosity  at  60°C.  Pa-s 

r  =====, 

240-1270 

55 

^ePor*s  on  variability  did  not  address  normality 

Portland  Cement  Concrete 


The  SDs  for  percent  passing  large  aggregate  sieve  sizes  is  typically  larger 

COnC,ret£  (PCC)  ^  for  asphalt  concrete  (AC),  as  shown 
m  lable  4-10.  This  trend  may  be  a  reflection  of  the  less  strict  aggregate 

control  requirements  for  PCC  construction,  as  compared  to  AC  construction. 
Ready-mix  PCC  plants  typically  include  only  two  aggregate  feed  bins,  while 
asphalt  plants  typically  include  four  bins. 


The  average  reported  CV  for  compressive  strength  was  15  percent  as 
shown  m  Table  4-11.  This  value  reflects  favorable  quality  control  according 
to  Baker  and  McMahon  (1969).  In  cases  of  poor  quality  control,  CV  for 
compressive  strength  exceeds  25  percent.  The  average  reported  CV  for  flex¬ 
ural  strength  was  7  percent,  which  reflects  favorable  quality  control  accord¬ 
ing  to  Witczak  et  al.  (1983).  In  cases  of  poor  quality  control  CV  for 
flexural  strength  exceeds  20  percent. 
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Table  4-8 

Summary  of  Reported  Variabilities  for  Asphalt  Concrete  Materials 

Standard 

Coefficient  of 

Property 

Deviation 

Variation,  Percent 

Normality 

1  Mixture  temperature  in  the 

8 

6 

not  reported1 

“■ h 

5' 

3- 

o 

Extracted  asphalt  cement 

0.25 

4 

not  reported1 

content,  percent 

Particles  finer,  percent 

not  reported1 

-  25-mm  sieve 

5 

5 

-  1 9-mm  sieve 

2 

2 

-  1 3.2-mm  sieve 

2 

2 

-  9.5-mm  sieve 

2 

2-3 

-  No.  4  sieve 

4 

1-10 

-  No.  8  sieve 

3 

5-10 

-  No.  1 6  sieve 

3 

5-10 

-  No.  30  sieve 

2 

5-10 

-  No.  50  sieve 

2 

10 

-  No.  100  sieve 

1 

10 

-  No.  200  sieve 

1 

10-20 

Bulk  specific  qravity 

0.030 

1 

not  reported1 

Theoretical  maximum 

0.010 

1 

not  reported1 

specific  qravity 

Voids  total  mix,  percent 

1.0 

10-65 

not  reported1 

Voids  filled,  percent 

5 

6 

not  reported1 

Field  density,  percent 

1.4 

2 

not  reported1 

relative  to  laboratory 

Marshall  stability,  kN 

1.3 

15 

not  reported1 

Marshall  flow,  mm 

0.30 

15 

not  reported1 

Indirect  tensile  strength,  kPa 

85 

15 

not  reported1 

Indirect  static  modulus. 

80-480 

55 

not  reported1 

MPa 

Indirect  poissons  ratio 

N/A 

40 

not  reported1 

Dynamic  modulus,  MPa 

not  reported1 

-  4SC 

N/A 

1 5 

-  21’C 

N/A 

30 

-  38*0 

N/A 

20 

Flexural  stiffness,  MPa 

not  reported1 

-  4-C 

520-1020 

20 

-  21°C 

150-280 

20 

Fatique  life  at  21  °C,  cycles 

N/A 

65 

not  reported1 

'  Reports  on  variability  did  not  address  normality. 

II  N/A  Limited  data  or  not  reported. 
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Table  4-9 

Summary  of  Reported  Variabilities  for  Field  Measurements  on 


Property 

Standard 

Deviation 

Coefficient  of 
Variation,  Percent 

Normality 

Mat  thickness,  mm 

8 

10 

not  reported’ 

Benkelman  beam 
deflections,  mm 

2 

15-70 

not  normal2 

Falling  weight  deflectometer 
(FWD) 

Deflections 

-  under  load,  fjm 

-  at  0.9  m  offset,  /jm 

30-340 

8-30 

15 

15 

normal2 

Falling  weight  deflectometer 
(FWD) 

Modulus  calculations 

-  AC  modulus,  MPa 

-  base  modulus,  MPa 

-  subgrade  modulus,  MPa 

210-3450 

20-240 

2-115 

50 

40 

15 

not  normal2 

Initial  serviceability  index 

0.35 

8 

normal3 

Reports  on  variability  did  not  address  normality. 

Based  on  the  shapes  of  histograms. 

3  According  to  the  chi-square  "qoodness-of-fit”  test 

Coefficients  of  variation  for  load  transfer  ranged  from  20  to  55  percent 
Tied  contraction  joints  had  the  lowest  variability.  Keyed  construction  joints 
and  doweled  contraction  joints  had  the  highest  variabilities.  Variabilities  in 
falling  weight  deflectometer  (FWD)  deflections  and  FWD  moduli  were  simi¬ 
lar  to  those  for  flexible  pavements. 


Frequency  distributions  for  most  of  the  fresh  and  hardened  concrete  prop- 
erhes  have  been  reported  to  be  normal  in  shape,  as  shown  in  Tables  4-10  and 
4-11.  The  measured  air  content  of  air-entrained  fresh  concrete  was  an 
exception;  its  distribution  was  reported  to  have  a  negative  skew.  The  FWD 
moduli  were  also  an  exception;  in  some  cases  they  exhibited  positive  skew. 
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Table  4-10 

Summary  of  Reported  Variabilities  for  Portland  Cement  Concrete 
Materials 

Standard 

Coefficient  of 

Property 

Deviation 

Variation,  Percent 

Normality 

Air  content,  percent 

0.85 

15 

not  normal1 

Slump,  mm 

8-25 

30 

normal1 

Particles  finer,  percent 

normal1 

(coarse  aggregates) 

-  25-mm  sieve 

3 

2-4 

-  1 9-mm  sieve 

10 

4-20 

-  1  3.2-mm  sieve 

10 

15-35 

-  No.  4  sieve 

1 

45-90 

-  No.  8  sieve 

1 

65 

Particles  finer,  percent 

normal1 

(fine  aggregates) 

-  No.  4  sieve 

1 

1 

-  No.  8  sieve 

3 

3 

-  No.  1 6  sieve 

4 

6 

-  No.  30  sieve 

4 

9 

-  No.  50  sieve 

5 

1  5-40 

-  No.  100  sieve 

1 

20-60 

Fineness  modulus 

0.10 

4 

normal1 

1  Based  on  the  shapes  of  histograms. _ _ 
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Table  4-1 1  ~~j 

Summary  of  Reported  Variabilities  for  Hardened  Portland  Cement 
|  Concrete 

Property 

Standard 

Deviation 

Coefficient  of 
Variation,  Percent 

Normality 

Density  in  hardened  state, 
kg/m3 

N/A 

2 

not  reported1 

Compressive  strenqth,  MPa 

4 

15 

normal2 

Compressive  modulus,  MPa 

N/A 

30 

not  reported1 

Poisson's  ratio 

N/A 

15 

not  reported1 

Indirect  tensile  strength,  kPa 

570 

15 

normal2 

Indirect  tensile  modulus, 

MPa 

8720 

35 

normal2 

Flexural  strenqth,  kPa 

350 

7 

normal2 

Slab  thickness,  mm 

9 

3 

normal2 

Load  transfer,  percent 

-  doweled  expansion 

-  doweled  construction 

-  keyed  construction 

-  doweled  contraction 

-  tied  contraction 

-  plain  contraction 

7 

5 

10 

4 

5 

4 

25 

25 

40 

4-55 

20 

25 

conflicts3 

Falling  weight  deflectometer 
(FWD) 

Deflections 
-  under  load,  pm 

3-120 

15 

normal2 

Falling  weight  deflectometer 
(FWD) 

Modulus  calculations 
*  PCC  modulus,  MPa 
-  subgrade  modulus,  MPa 

2830-15740 

15-30 

30 

15 

not  reported1 

Serviceability  index 

0.15 

3 

normal2 

Reports  on  variability  did  not  address  normality. 

Based  on  the  shapes  of  histograms. 

3  Reports  concerning  normality  conflicted. 

N/A  Limited  data  or  not  reported. 
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5  Summary,  Conclusions, 
and  Recommendations 


Summary  and  Conclusions 


The  paving  industry  is  moving  forward  with  performance-related  specifi¬ 
cations  that  incorporate  statistical  concepts  for  the  development  of  acceptance 
plans,  including  quality  control  and  quality  assurance  testing.  Personnel 
involved  in  pavement  construction  need  to  understand  these  statistical  con¬ 
cepts  and  they  need  to  know  the  levels  of  variability  that  should  be  expected 
in  their  test  results. 

This  report  provides  information  on  statistical  methods  that  can  be  used  by 
pavement  engineers  and  contractors  to  analyze  the  variability  of  construction 
materials.  In  addition,  it  summarizes  selected  published  data  on  these  vari¬ 
abilities.  These  data  were  measured  within  “lots”  of  material,  where  a  lot  is 
defined  as  a  quantity  of  material  that  can  be  considered  uniform  for  sampling 
and  payment  purposes. 

The  methods  for  analyzing  variability  will  be  useful  for  those  who  need  to 
quantify  the  relative  magnitudes  of  different  sources  of  variation.  If  a  con¬ 
tractor’s  product  is  highly  variable,  these  methods  can  be  used  to  determine 
if  the  production  of  materials  is  out  of  control  or  if  the  methods  for  sampling 
and/or  testing  are  inconsistent.  By  keeping  track  of  the  trends  in  test  results, 
the  contractor  will  be  able  to  adjust  production  and/or  construction  methods 
as  necessary  in  order  to  provide  a  high-quality,  consistent  product. 

The  data  pertaining  to  published  variabilities  provide  the  pavement 
designer  and  specification  writer  with  tools  for  determining  what  can  be  rea¬ 
sonably  expected  of  the  contractor.  A  contractor  cannot  be  expected  to  keep 
material  variation  well  below  the  published  norm.  Some  sources  of  variabil¬ 
ity,  such  as  those  inherent  to  a  particular  test  method,  are  out  of  the  contrac¬ 
tor’s  control.  Knowledge  of  expected  variability  will  also  help  the  owner 
control  two  important  aspects  of  specifications:  the  risks  for  contractual 
parties  and  the  reliability  of  the  final  product. 
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Recommendations 


The  purpose  of  this  report  was  to  review  statistical  concepts  and  compile 
data  related  to  the  variability  of  materials  used  for  the  construction  of  pave¬ 
ments.  Based  on  the  concepts  presented  and  data  reviewed,  the  following 
actions  are  recommended. 

a.  The  current  Corps  of  Engineers  specification  requirements  for  pave¬ 
ment  construction  should  be  reviewed  to  determine  the  appropriateness 
of  the  statistical  concepts  incorporated. 

b.  Where  necessary,  the  Corps  of  Engineers  specifications  should  be 
adjusted  to  reflect  reasonable  expectations  of  the  contractor  so  that  the 
risks  assumed  by  the  contractor  and  owner  are  equitable. 


c.  Corps  of  Engineer  specifications  that  do  not  contain  statistical  concepts 
for  acceptance  testing  should  be  modified. 
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A1 


Table  A1 

Probability  of  Obtaining  a  Random  Value  of  Z 
Margins  (after  Steel  and  Torrie  1980) 


.4602  .4562  .4522  .4483  .4443 


.4207  .4168  .4129  .4090  .4052 


.3821  .3783  .3745  .3707  .3669 


.3446  .3409  .3372  .3336  .3300 


.3085  .3050  .3015  .2981  .2946 


.2743  .2709  .2676  .2643  .2611 


.2420  .2389  .2358  .2327  .2296 


0.8  .2119  .2090  I  .2061  I  .2033  |  .2005 


0.9  .1841  .1814  .1788  .1762  .1736 


.1587  .1562  .1539  .1515 


Greater  Than  the  Values  Shown  in  the 


1.5 


1.6 


1.7 


1.8 


1.9 


2.0 


2.1 


2.2 


2.3 


2. 


2.5 


2.6 


.1515 

.1492 

.1469 

.1446 

.1423 

.1401 

.1379 

.1292 

.1271 

.1251 

.1230 

.1210 

.1190 

.1170 

.1093 

.1075 

.1056 

.1038 

.1020 

.1003 

.0985 

.0918 

.0901 

.0885 

.0869 

.0853 

.0838 

.0823 

.0764 

.0749 

.0735 

.0721 

.0708 

.0694 

.0681 

.0630 

.0618 

.0606 

.0594 

.0582 

.0571 

.0559 

.0516 

.0505 

.0495 

.0485 

.0475 

.0465 

.0455 

.0418 

.0409 

.0401 

.0392 

.0384 

.0375 

.0367 

.0336 

.0329 

.0322 

.0314 

.0307 

.0301 

.0294 

.0268 

.0262 

.0256 

.0250 

.0244 

.0239 

.0233 

.0212 

.0207 

.0202 

.0197 

.0192 

.0188 

.0183 

.0166 

.0162 

.0158 

.0154 

.0150 

.0146 

.0143 

.0129 

.0125 

.0122 

.0119 

.0116 

.0113 

.0110 

.0099 

.0096 

.0094 

.0091 

.0089 

.0087 

.0084 

.0075 

.0073 

.0071 

.0069 

.0068 

.0066 

.0064 

.0057 

.0055 

.0054 

.0052 

.0051 

.0049 

.0048 

.0043 

.0041 

.0040 

.0039 

.0038 

.0037 

.0036 

.0032 

.0031 

.0030 

.0029 

.0028 

.0027 

.0026 

.0023 

.0023 

.0022 

.0021 

.0021 

.0020 

.0019 
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Table  A1  (Concluded) 


2 

.00 

.01 

.02 

.03 

.04 

.05 

.06 

.07 

08 

2.9 

.0019 

.0018 

.0018 

.0017 

.0016 

.0016 

.0015 

.0015 

0014 

3.0 

.0013 

.0013 

.0013 

.0012 

.0012 

.001 1 

.0011 

.001 1 

0010 

3.1 

.0010 

.0009 

.0009 

.0009 

.0008 

.0008 

.0008 

.0008 

0007 

3.2 

.0007 

.0007 

.0006 

.0006 

.0006 

.0006 

.0006 

.0005 

0005 

3.3 

.0005 

.0005 

.0005 

.0004 

.0004 

.0004 

.0004 

.0004 

.0004 

.0003 

3.4 

.0003 

.0003 

.0003 

.0003 

.0003 

.0003 

.0003 

.0003 

.0003 

.0002 

S 

.0002 

.0002 

.0001 

.0001 

.0001 

.0001 

.0001 

.0001 

.0001 

.0001 
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A3 


Table  A2 

Values  of  t  (after  Fisher  and  Yates  1949) 


Probability  of  a  Numerically  Larger  Value  of  t  (Two-Tail  Test) 


0.05 


12.706 


4.303 


3.182 


2.776 


0.02 


31.821 


6.965 


4.541 


3.747 


63.657  636.619 


31.598 


12.941 


8.610 


5 

2.015 

2.571 

3.365 

4.032 

6.859 

6 

1.943 

2.447 

3.143 

3.707 

5.959 

7 

1.895 

2.365 

2.998 

3.499 

5.405 

8 

1.860 

2.306 

2.896 

3.355 

5.041 

9 

1.833 

2.262 

2.821 

3.250 

4.781 

10 

1.812 

2.228 

2.764 

3.169 

4.587 

1 1 

1.796 

2.201 

2.718 

3.106 

4.437 

12 

1.782 

2.179 

2.681 

3.055 

4.318 

13 

1.771 

2.160 

2.650 

3.012 

4.221 

14 

1.761 

2.145 

2.624 

2.977 

4.140 

15 

1.753 

2.131 

2.602 

2.947 

4.073 

16 

1.746 

2.120 

2.583 

2.921 

4.015 

17 

1.740 

2.110 

2.567 

2.898 

3.965 

18 

1.734 

2.101 

2.552 

2.878 

3.922 

19 

1.729 

2.093 

2.539 

2.861 

3.883 

20 

1.725 

2.086 

2.528 

2.845 

3.850 

df 

0.05 

0.025 

0.01 

0.005 

0.0005 

Probability  of  a  Larger  Positive  Value  of  t  (One-Tail  Test) 

|  (Continued)  \ 
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Table  A4 

Values  of  F  (after  Steel  and  Torrie  1980) 


1  . — 

Denominator 

df 

Numerator  df 

Problem  of 
a  Larger  F 

1 

2 

4 

6 

1 

.050 

161.4 

199.5 

224.6 

234.0 

.025 

647.8 

799.5 

899.6 

937.1 

.010 

4,052 

5,000 

5,625 

5,859 

2 

.050 

18.51 

19.00 

19.25 

19.33 

.025 

38.51 

39.00 

39.25 

39.33 

.010 

98.50 

99.00 

99.25 

99.33 

4 

.050 

7.71 

6.94 

6.39 

6.16 

.025 

12.22 

10.65 

9.60 

9.20 

.010 

21.20 

18.00 

15.98 

15.21 

6 

.050 

5.99 

5.14 

4.53 

4.28 

.025 

8.81 

7.26 

6.23 

5.82 

.010 

13.75 

10.92 

9.15 

8.47 

8 

.050 

5.32 

4.46 

3.84 

3.58 

.025 

7.57 

6.06 

5.05 

4.65 

.010 

11.26 

8.65 

7.01 

6.37 

10 

.050 

4.96 

4.10 

3.48 

3.22 

.025 

6.94 

5.46 

4.47 

4.07 

.010 

10.04 

7.56 

5.99 

5.39 

20 

.050 

4.35 

3.49 

2.87 

2.60 

.025 

5.87 

4.46 

3.51 

3.13 

.010 

8.10 

5.85 

4.43 

3.87 

OO 

.050 

3.84 

3.00 

2.37 

2.10 

.025 

5.02 

3.69 

2.79 

2.41 

.010 

6.63 

4.61 

3.32 

2.80 

1 

(Continued) 
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Table  A5 

Critical  Values  for  T  When  Standard  Deviation 
the  Sample  [after  ASTM  (1995b)  E  178] _ 


is  Calculated  from 


0.1% 

_  J  . 

0.5% 

1 .0% 

2.5% 

5% 

10% 

1.155 

1.155 

1.155 

1.155 

1.153 

1.148 

1.499 

1.496 

1.492 

1.481 

1.463 

1.425 

1.780 

1.764 

1.749 

1.715 

1.672 

1.602 

2.01 1 

1.973 

1.944 

1.887 

1.822 

1.729 

2.201 

2.139 

2.097 

2.020 

1.938 

1.828 

2.358 

2.274 

2.221 

2.126 

2.032 

1.909 

2.492 

2.387 

2.323 

2.215 

2.110 

1.977 

2.606 

2.482 

2.410 

2.290 

2.176 

2.036 

2.997 

2.806 

2.705 

2.549 

2.409 

2.247 

3.230 

3.001 

2.884 

2.709 

2.557 

2.385 

3.389 

3.135 

3.009 

2.822 

2.663 

2.486 

3.507 

3.236 

3.103 

2.908 

2.745 

2.563 

3.673 

3.381 

3.240 

3.036 

2.866 

2.682 

3.789 

3.483 

3.336 

3.128 

2.956 

2.768 

3.874 

3.560 

3.411 

3.199 

3.025 

2.837 

3.942 

3.622 

3.471 

3.257 

3.082 

2.893 

3.998 

3.673 

3.521 

3.305 

3.130 

2.940 

3.998 

3.673 

3.521 

3.305 

3.130 

2.940 

4.084 

3.754 

3.600 

3.383 

3.207 

3.017 

4.044 

3.716 

3.563 

3.347 

3.171 

2.981 

4.150 

3.817 

3.662 

3.444 

3.267 

3.078 

4.178 

3.843 

3.688 

3.470 

3.294 

3.104 

4.203 

3.867 

3.712 

3.493 

3.318 

3.129 
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Table  A7 

I  Values  of  Fmax  (after  Ott  1977) 

df 

each 

sample 

Problem 

Number  of  Independent  Mean  Squares  1! 

of  a 

Larger  Fm„ 

2 

3 

4 

5 

mm 

2 

.050 

39.0 

87.5 

142 

202 

266 

.010 

199 

448 

729 

1,036 

1,362 

3 

.050 

15.4 

27.8 

39.2 

50.7 

62.0 

.010 

47.5 

85 

120 

151 

184 

4 

.050 

9.60 

15.5 

20.6 

26.2 

29.5 

.010 

23.2 

37 

49 

59 

69 

5 

.050 

7.15 

10.3 

13.7 

16.3 

18.7 

.010 

14.9 

22 

28 

33 

38 

6 

.050 

5.82 

8.38 

10.4 

12.1 

13.7 

.010 

11.1 

15.5 

19.1 

22 

25 

7 

.050 

4.99 

6.94 

8.44 

9.70 

10.8 

.010 

8.89 

12.1 

14.5 

16.5 

18.4 

8 

.050 

4.43 

6.00 

7.18 

8.12 

9.03 

.010 

7.50 

9.9 

11.7 

13.2 

14.5 

9 

.050 

4.03 

5.34 

6.31 

7.11 

7.80 

.010 

6.54 

8.5 

9.9 

11.1 

12.1 

10 

.050 

3.72 

4.85 

5.67 

6.34 

6.92 

.010 

5.85 

7.4 

8.6 

9.6 

10.4 

12 

.050 

3.28 

4.16 

4.79 

5.30 

5.72 

.010 

4.91 

6.1 

6.9 

CD 

8.2 

15 

.050 

2.86 

3.54 

4.01 

4.37 

4.68 

.010 

4.07 

4.9 

5.5 

6.0 

6.4 

30 

.050 

2.07 

2.40 

2.61 

2.78 

2.91 

.010 

2.63 

3.0 

3.3 

3.4 

3.6 

00 

.050 

1.00 

1.00 

1.00 

1.00 

1.00 

.010 

1.00 

1.00 

1.00 

1.00 

1.00 

(Continued) 
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Table  A7  (Concluded) 


Problem 
of  a 

e  Larqer  F 


.050 


.010 


.050 


.010 


.050 


.010 


.050 


.010 


.050 


.010 


.050 


.010 


.050 


.010 


.050 


.010 


.050 


.010 


.050 


.010 


.050 


.010 


.050 


.010 


.050 


.010 


333 


1,705 


72.9 


216 


Number  of  Independent  Mean  Squares 


403 


2,063 


83.5 


249 


33.6 

37.5 

44.6 

51.4 

58.1 

79 

89 

106 

120 

137 

20.8 

22.9 

26.5 

29.9 

32.6 

42 

46 

54 

60 

66 

15.0 

16.3 

18.6 

20.7 

22.1 

27 

30 

34 

37 

39 

11.8 

12.7 

14.3 

15.8 

16.6 

20 

22 

24 

27 

32 

9.8 

10.5 

11.7 

12.7 

13.1 

15.8 

16.9 

18.9 

21 

24 

8.41 

8.95 

9.91 

10.7 

10.9 

13.1 

13.9 

15.3 

16.6 

18.2 

7.42 

7.87 

8.66 

9.34 

9.45 

11.1 

11.8 

12.9 

13.9 

15.1 

6.09 

6.42 

7.00 

7.48 

8.03 
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Appendix  B  Residual  Fine-Grained  Soil  Deposits 


B1 


n 

Table  B1 

Specific  Gravity  _  _ 

Mean 

Standard 

Deviation 

Coefficient  of 
Variation,  Percent 

Reference 

2.86 

0.03 

1.0 

Schultze  (1 972) 1 

1  N/A 

N/A 

2.0 

Padilla  and  Vanmarcke  (1974) 

1  Normality 
|  N/A  Data  no 

was  not  rejected  by  the  chi  square  "goodness-of-fit“  test, 
t  reported.  -  - - — 

Table  B2 

Natural  Water  Content,  Percent _ — ==J 

Mean 

Standard 

Deviation 

Coefficient  of 
Variation,  Percent 

Reference 

14 

3.5 

26 

Krahn  and  Fredlund  (1983) 

21 

4.7 

23 

Schultze  (1 972)1 

22 

2.8 

13 

27 

6.0 

22 

Fredlund  and  Dahlman  (1972)2 

29 

5.8 

20 

Krahn  and  Fredlund  (1983) 

29 

5.1 

18 

Fredlund  and  Dahlman  (1972)2 

32 

4.4 

14 

Krahn  and  Fredlund  (1983) 

33 

4.1 

12 

Fredlund  and  Dahlman  (1972)2 

35 

3.5 

10 

35 

4.6 

13 

36 

4.3 

12 

1  Normality  was  rejected  by  the  chi  square  "goodness-of-fit  test. 

2  Frequency  distributions  appeared  to  _te_normaj_in_shape. _ -—==—== 

B2 
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Table  B3 
Porosity 


Mean 

Standard 

Deviation 

Coefficient  of 

Variation,  Percent 

0.33 

0.06 

19 

0.38 

0.05 

13 

0.39 

0.04 

10 

0.40 

0.03 

7.5 

0.48 

0.05 

10 

Reference _ 

(Schultze  1972)’ 


Normality  was  not  rejected  by  the  chi  square  "qoodness-of-fit"  test. 


Table  B4 
Void  Ratio 


Mean 


0.42 


0.51 


0.53 


0.56 


0.57 


0.66 


0.75 


0.90 


0.91 


0.91 


0.95 


Standard 

Deviation 


Coefficient  of 

Variation,  Percent  Reference 


(Krahn  and  Fredlund  1983) 
(Schultze  1972)' 


(Fredlund  and  Dahlman  1972) 


u.  M  _LJ4  MKrahnandFredlund  1983) 

Normality  was  not  rejected  by  the  chi  square  "goodness-of-fit”  test. 
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Table  B5 

Dearee  of  Saturation,  Percent 

Mean 

Standard 

Deviation 

Coefficient  of 
Variation,  Percent 

Reference 

78.7 

11.4 

15 

Krahn  and  Fredlund  (1983) 

81.3 

10.6 

13 

85 

16 

19 

Schultze  (1 972)1 

91.5 

6.0 

6.6 

Krahn  and  Fredlund  (1983) 

91.9 

12.1 

13 

Fredlund  and  Dahlman  (1972) 

93.3 

8.0 

8.5 

94.1 

10.0 

11 

Krahn  and  Fredlund  (1983) 

95.6 

8.6 

9.0 

Fredlund  and  Dahlman  (1972) 

97.5 

8.4 

8.6 

1  Normality  was  rejected  by  the  ch 

i  square  "qoodness-of-fit"  test. 

Table  B6 

Drv  Density,  kq/m3  - ____ 

Mean 

Standard 

Deviation 

Coefficient  of 
Variation,  Percent 

Reference  | 

1,430 

99 

6.9 

(Krahn  and  Fredlund  1983) 

1,880 

106 

5.6 

Table  B7 

Coefficient  of  Permeability,  Percent 

1  Mean 

Standard 

Deviation 

Coefficient  of 
Variation,  Percent 

Comment 

Reference 

N/A 

N/A 

90 

80  percent  saturated 

Nielson  et  al. 

N/A 

N/A 

700 

1 00  percent  saturated 

(1973) 

N/A  Data  not  reported. 
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Table 
1  Liquic 

B8 

Limit,  Percent 

|  Mean 

Standard 

Deviation 

Coefficient  of 
Variation,  Percent 

Reference 

16 

9.7 

62 

Mitchell  et  at.  (1977) 

20 

4.0 

20 

27 

5.4 

20 

27 

6.0 

22 

Krahn  and  Fredlund  (1983) 

33 

6.3 

19 

Kennedy,  Hudson,  and  McCullouah  (197RI 

43 

11 

26 

Wahls  and  Futrell  (1966) 

53 

11 

21 

Fredlund  and  Dahiman  (1972)1 

55 

12 

22 

Schultze  (1972)2 

57 

13 

24 

Krahn  and  Fredlund  (1983) 

59 

13 

22 

Fredlund  and  Dahiman  (1972)1 

62 

12 

19 

Krahn  and  Fredlund  (1983) 

63 

14 

22 

Wahls  and  Futrell  (1966) 

63 

11 

18 

Fredlund  and  Dahiman  (1972)' 

64 

11 

17 

73 

13 

18 

Krahn  and  Fredlund  (1983) 

===============================J 

’  Frequency  distributions  appeared  normal  in  shape. 

Normality  was  not  rejected  by  the  chi  square  ”goodness-of-fit"  test. 

Note:  Sampling  error  and  experimental  error  (testing  and  inherent)  accounted  for 

40  percent  of  the  total  variance  (Hampton.  Yoder,  and  Burr  1 9621 
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Table  B9 

Plastic  Limit.  Percent _ _ _ = 

Mean 

Standard 

Deviation 

Coefficient  of 

Variation,  Percent 

Reference 

14 

2.2 

16 

Krahn  and  Fredlund  (1983) 

21 

3.4 

16 

Schultze  (1972)' 

23 

3.8 

17 

Krahn  and  Fredlund  (1983) 

24 

3.2 

13 

Fredlund  and  Dahlman  (1972)2 

25 

3.8 

15 

25 

3.2 

13 

25 

3.4 

13 

Krahn  and  Fredlund  (1983) 

26 

3.4 

13 

Fredlund  and  Dahlman  ( 1 97 2)2 

27 

4.8 

19 

Krahn  and  Fredlund  (1983) 

’  Normality  was  rejected  by  the  chi  square  "goodness-of-fit"  test. 

2  Frequency  distributions  appeared  normal  in  shape. 

Note:  Sampling  error  and  experimental  error  (testing  and  inherent)  accounted  tor 

35  nercent  of  the  total  variance  (Hampton,  Yoder,  and  Burr  1962). - _l 

Table  BIO 

Plasticitv  Index.  Percent 

Mean 

Standard 

Deviation 

Coefficient  of 
Variation.  Percent 

Reference 

4.3 

3.4 

80 

Mitchell  et  al.  (1977) 

10 

4.8 

48 

Kennedy,  Hudson,  and 

McCullouqh  (1975) 

12 

7.2 

59 

Wahls  and  Futrell  (1966) 

14 

5.9 

42 

Kennedy,  Hudson,  and 

McCullough  (1975) 

16 

2.6 

16 

27 

10.2 

39 

Wahls  and  Futrell  (1966) 

33.9 

9.6 

28 

Schultze  (1972)’ 

’  Normality  was  rejected  by  the  chi  square  goodness-of-fit  test. 

Notes:  (1)  Frequency  distributions  appeared  normal  in  shape  (Ingles  1  972). 

(2)  Sampling  error  and  experimental  error  (testing  and  inherent)  accounted  for 

_ 60_£ercentofthetotalvarianceJHamg^n^Yoder^nd_Burr_^962^__________J 
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Table  B1 1 

Particle  Size  Distribution,  Percent  Finer 

Sieve  Size 

Mean 

Standard 

Deviation 

Coefficient  of 
Variation,  Percent 

Reference 

No.  40 

79 

16 

21 

Wahls  and  Futrell  (1966) 

91 

9.4 

10 

No.  200 

14 

2.4 

17 

Auff  and  Choummanivong 
(1994) 

56 

20 

36 

Wahls  and  Futrell  (1966) 

77 

7.0 

9.1 

Hampton,  Yoder,  and  Burr 
(1962)' 

79 

16 

20 

Wahls  and  Futrell  (1966) 

2  pm 

17 

11 

62 

Wahls  and  Futrell  (1966) 

27 

3.5 

13 

Hampton,  Yoder,  and  Burr 
(1962)2 

41 

12 

30 

Wahls  and  Futrell  (1966) 

of  total  variance1  eXperimental  error  (testin9  and  inherent)  accounted  for  50  percent 

oftXtal  variance!  eXperimental  err0r  (testin9  and  inherent)  accounted  for  65  percent 

Table  B12 
Potential  ft 

>r  Volume  Change 

Measure 

Mean 

Standard 

Deviation 

Coefficient  of 
Variation,  Percent 

Reference 

Compression 

index 

0.108 

0.032 

30 

Krahn  and  Fredlund 
(1983) 

0.110 

0.052 

47 

Fredlund  and  Dahlman 
(1972) 

0.155 

0.049 

33 

Krahn  and  Fredlund 
(1983) 

0.159 

0.048 

30 

Fredlund  and  Dahlman 
(1972) 

0.167 

0.048 

29 

0.184 

0.047 

26 

0.205 

0.066 

32 

0.265 

0.070 

26 

(1983) 

Swell  index 

0.032 

0.017 

53 

Fredlund  and  Dahlman 

0.065 

0.034 

52 

(1972) 

0.065 

0.020 

31 

0.065 

0.017 

26 
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Table  B13 

Moisture-Density  Relationship 

Measure 

Mean 

Standard 

Deviation 

Coefficient  of 

Variation, 

Percent 

Reference 

Optimum  moisture 
content,  percent 

19.8 

3.5 

18 

Willenbrock  (1 974) 1 

20.3 

5.6 

28 

Wahls  and  Futrell 
(1966) 

26.4 

5.4  , 

20 

Maximum  dry 
density,  kg/m3 

1,495 

132 

8.8 

Wahls  and  Futrell 
(1966) 

1,655 

164 

9.9 

1,690 

98 

5.8 

Willenbrock  (1974)2 

'  Frequency  distributions  appeared  normal  in  shape. 

2  Frequency  distributions  appeared  to  have  negative  skew. 

Note:  Sampling  error  and  experimental  error  (testing  and  inherent)  accounted  for 

40  percent  and  25  percent  of  the  total  variance  for  optimum  moisture  content  and 

1  maximum  density,  respectively  (Hampton^Yoder^jind^BinMgSg^ 

Table  B14 

Unconfined  Compressive  Strength  (kPa) 


— 

Mean 

Standard 

Deviation 

Coefficient  of 
Variation,  Percent 

Reference 

21 

7.0 

33 

Ladd,  Moe,  and  Gifford  (1972) 

26 

14 

53 

Schultze  (1 972)1 

40 

21 

52 

58 

23 

39 

125 

59 

48 

Fredlund  and  Dahlman  (1972)2 

143 

57 

40 

161 

66 

41 

167 

75 

45 

Krahn  and  Fredlund  (1983) 

199 

98 

49 

Fredlund  and  Dahlman  ( 1 97 2)2 

220 

122 

55 

Krahn  and  Fredlund  (1983) 

253 

179 

71 

331 

86 

26 

Morse  (1 972) 

370 

188 

51 

Krahn  and  Fredlund  (1983) 

1  Normality  was  rejected  by  the  chi  square  ”goodness-of-fit”  test. 

2  Frequency  distribution  appeared  to  have  a  positive  skew. 

Note:  Frequency  distributions  appeared  normaMn_shape_(lr^les_T£72)1 
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Table  B15 

Measures  of  Shear  Strength 

Measure 

Mean 

Standard 

Deviation 

Coefficient  of 

Variation, 

Percent 

Reference 

California 
bearing  ratio 

1 0  percent 

2.7  percent 

27 

Hampton,  Yoder,  and 

Burr  (1962)’ 

Field  vane 

26  kPa 

7.6  kPa 

29 

(19  to  35) 

Ladd,  Moe,  and  Gifford 
(1972) 

Undrained 
triaxial  test,  c 

N/A 

N/A 

1 7  to  432 

Kay  and  Krizek  (1972) 

N/A 

N/A 

19 

Lumb  (1972) 

Undrained 
triaxial  test, 
tan<(> 

N/A 

N/A 

22 

Drained 
triaxial  test,  c 

N/A 

N/A 

14 

Drained 
triaxial  test, 
tan<(> 

N/A 

N/A 

1.6 

Drained  shear 
box  test,  c 

N/A 

N/A 

95  to  1033-4 

Drained  shear 
box  test,  tan((> 

N/A 

N/A 

1 8  to  463,5 

-  - 

the  total' varfance3"0*  eXperimental  error  (testinS  a"d  inherent)  accounted  for  80  percent  of 

2  Depends  on  sampling  technique. 

3  Depends  on  type  of  soil. 

a.nd  ®xperimentai  err°r  (testing  and  inherent)  accounted  for  20  to 

50  percent  of  the  total  variance. 

"  S^P'int9  e.rr°r  and  experimental  error  (testing  and  inherent)  accounted  for  10  to 

25  percent  of  the  total  variance. 

N/A  Data  not  reported. 
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Cl 


Table  Cl 

Controlled  Moisture  Content,  Percent 

Mean 

Standard 

Deviation 

Coefficient  of 
Variation,  Percent 

Reference 

8.3 

2.6 

31 

Auff  and  Choummanivong  (1994) 

9.1 

1.3 

14 

Seiiq  (1966) 

9.2 

1.0 

10 

Yeo  and  Auff  (1995) 

11 

1.1 

11 

12 

2.7 

23 

Auff  and  Choummanivong  (1994) 

12 

2.0 

16 

Seiiq  (1966) 

15 

3.6 

24 

David  (1967)’ 

1  Frequency  distribution  appeared  to  have  positive  skew. 

Notp-  Freauencv  distributions  appeared  normal  in  shape  (Baecher  1987). 
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Table  C2 
1  Field  Cot 

| 

iipaction.  Percent  Relative  to  Laboratory  Density 

|  Mean 

Standard 

Deviation 

Coefficient  of 

Variation,  Percent 

Reference 

87.0 

2.5 

2.8 

Yeo  and  Auff  {1  995) 

87.2 

6.1 

7.0 

Jorgensen  (1968)1 

88.7 

4.5 

5.1 

88.9 

4.1 

4.6 

89.9 

8.8 

9.8 

90.5 

3.1 

3.4 

Sherman,  Watkins,  and  Prysock 
(1966) 

91.5 

3.8 

4.2 

Yeo  and  Auff  (1995) 

92.8 

2.9 

3.1 

Auff  and  Choummanivong 
(1994) 

92.8 

4.8 

5.2 

92.9 

2.4 

2.6 

Sherman,  Watkins,  and  Prysock 
(1966)2 

93.6 

5.5 

5.9 

94.5 

3.0 

3.2 

Brown  (1975) 

96.6 

2.9 

2.9 

Mitchell  et  al.  (1977) 

96.8 

2.5 

2.6 

96.8 

5.7 

5.9 

Williamson  and  Yoder  (1967)3 

97.3 

3.7 

3.8 

David  (1967)3 

97.8 

4.9 

5.0 

Jorgensen  (1968) 

98.2 

4.5 

4.6 

Williamson  and  Yoder  (1967) 

98.2 

3.2 

3.3 

Mitchell  et  al.  (1977) 

99.0 

4.6 

4.6 

Nielson  (1 967)“ 

99.1 

4.5 

4.5 

Van  Houten  (1 967)1 

100.6 

5.3 

5.3 

Williamson  and  Yoder  I1.967)2 

of  5aoTiir~nT°/rd,ef,erimentai  6rr0r  (t6Stin9  and  inherent)  accounted  for  an  average 
of  50  percent  of  the  total  variance  when  tested  by  water  balloon.  9 

3  Frequency  distributions  did  not  appear  to  be  normal  in  shape. 

4  Frequency  distributions  appeared  to  be  normal  in  shape 

of10T:noLe;r0frahndef,erimental  err0r  (testi"9  and  inherent)  accounted  for  an  average 
of  10  percent  of  the  total  variance  when  tested  bv  sand  cone  9 
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Table  C3  1 

Resilient  Modulus,  MPa  _ _ _  -J 

Mean 

Standard 

Deviation 

Coefficient  of 
Variation,  Percent 

Reference 

N/A 

N/A 

20 

Inqles  (1972) 

83 

24 

29 

Kennedy,  Hudson, and 
McCullough  (1975) 

110 

24 

22 

131 

29 

22 

N/A  Data  not  reported.  - — - - - - ' 

Table  C4 

Measures  of  Stiffness  and  Strength 


■ - 

Measure 

Mean 

Standard 

Deviation 

Coefficient  of 
Variation,  Percent 

Reference 

California 

7 

1.5 

22 

Yoder  and  Witczak 

bearing 

8 

1.4 

18 

(1975) 

12 

4.1 

34 

Yeo  and  Auff 
(1  995)1 

18 

4.7 

26 

Yoder  and  Witczak 
(1975) 

21 

6.7 

32 

Yeo  and  Auff 
(1995) 

43 

15 

35 

Mitchell  et  al. 

(1977) 

Modulus  of 

N/A 

N/A 

33 

Yoder  and  Witczak 
(1975) 

MPa/mm 

N/A 

N/A 

35 

Highway  Research 
Board  (1 962b) 

N/A 

N/A 

41 

Kennedy,  Hudson, 
and  McCullough 
(1975) 

Benkelman  beam 
deflections,  mm 

1.4 

0.35 

25 

Yeo  and  Auff 
(1995) 

Dynaflect  tests, 

MPa 

N/A 

N/A 

34 

Kennedy,  Hudson, 
and  McCullough 
(1975) 

1!  - 

1  Estimated  by  dynamic  cone  penetrometer  (DCP)  tests. 
N/A  Data  not  reported. _ _ 
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Table  D2 

Atterberg  Limits,  Percent 


Property 

Mean 

Standard 

Deviation 

Coefficient  of 
Variation,  Percent 

Reference 

Liquid  limit 

12 

9.8 

81 

Mitchell  et  al.  (1977) 

16 

20 

128 

17 

2.2 

13 

18 

6.4 

37 

20 

1.8 

9.3 

Plasticity 

2.9 

1.0 

33 

Mitchell  et  al.  (1977) 

index 

3.0 

1.3 

44 

3.5 

2.3 

65 

3.9 

2.6 

67 

4.0 

1.0 

26 

Ingles  (1972)' 

5.9 

1.4 

24 

Mitchell  et  al.  (1977) 

Frequency  distribution  appeared  tohave  positive  skew. 


Table  D 
Coeffici 

3 

ent  of  Peri 

neability.  Percent 

Mean 

Standard 

Deviation 

Coefficient  of 
Variation,  Percent 

Comment 

Reference  | 

I  N/A 

N/A 

Silty  sand 

Lumb  (1972)'  1 

thp  tntP  9  K  expenmental  error  (testing  and  inherent)  accounted  for  1  percent  of 

the  total  variance  when  tested  by  constant  head  permeameter 

N/A  Data  not  reported.  | 
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Table  D4 

Particle  Size  Distribution,  Percent  Finer 


Sieve  Size 
1  9-mm 


1 3.2-mm 


9.5-mm 


Mean 

Standard 

Deviation 

Coefficient  of 

Variation, 

Percent 

Reference 

20 

0.9 

4.5 

Yeo  and  Auff  (1995) 

22 

2.4 

11 

22 

2.2 

10 

22 

2.5 

11 

70 

7.3 

10 

Willenbrock  (1974b) 

80 

3.9 

4.9 

Kelley  (1969)’ 

90 

2.5 

2.8 

85 

2.3 

2.7 

Auff  and  Yeo  (1992) 

86 

1.8 

2.1 

42 

8.1 

19 

Willenbrock  (1974b) 

50 

5.2 

10 

Kelley  (1969)2 

70 

2.9 

4.1 

Auff  and  Yeo  (1992) 

71 

2.2 

3.1 

30 

6.8 

22 

Willenbrock  (1974b) 

34 

4.3 

13 

Kelley  (1969)1 

45 

6.6 

15 

State  of  California  (1967) 

50 

4.3 

8.7 

51 

3.1 

6.1 

53 

5.7 

11 

53 

2.9 

5.5 

Auff  and  Yeo  (1992) 

54 

4.1 

7.6 

Sherman  (197 1)2 

55 

4.9 

8.9 

Kelley  (1969) 

55 

2.0 

3.6 

Auff  and  Yeo  (1992) 

56 

5.8 

10 

Sherman  ( 1 97 1  )3 

58 

2.8 

4.8 

State  of  California  (1967) 

73 

6.5 

9.0 

1  Sampling  error  and  experimental  error  (testing  and  inherent) 
of  the  total  variance. 

2  Sampling  error  and  experimental  error  (testing  and  inherent) 
of  the  total  variance. 

3  Sampling  error  and  experimental  error  (testing  and  inherent) 

of  the  total  variance.  — 


_  (Sheet  1  of  3)  \ 

accounted  for  40  percent 
accounted  for  35  percent 
accounted  for  20  percent 
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1  Table  D4  (Concluded)  _ _ _  - 

Sieve  Size 

Mean 

Standard 

Deviation 

Coefficient  of 

Variation, 

Percent 

Reference 

No.  200 
(Cont.) 

7.4 

1.9 

26 

Auff  and  Choummanivong 
(1994) 

7.8 

1.3 

17 

State  of  California  (1967) 

7.9 

1.1 

14 

8.6 

1.7 

20 

8.8 

1.6 

18 

Sherman  ( 1 97 1  )2 

9.0 

0.6 

6.7 

Auff  and  Yeo  (1992) 

9.0 

0.8 

9.0 

9.0 

0.6 

6.7 

Auff  and  Laksmanto 

10 

1.1 

11 

(1994) 

10 

1.8 

18 

State  of  California  (1967) 

13 

1.2 

9.2 

Auff  and  Laksmanto 
(1994) 

13 

0.9 

7.0 

Auff  and  Choummanivong 

16 

1.8 

11 

(1994) 

16 

2.9 

18 

Kelley  (1969) 

20 

0.9 

4.5 

Yeo  and  Auff  (1995) 

21 

2.4 

11 

22 

2.2 

10 

22 

2.5 

11 

1  (Sheet  3  of  3)  \ 

2  Sampling  error  and  experimental  error  (testing  and  inherent)  accounted  for  35  percent 
|  of  the  total  variance^^ ======  1 
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Table  D7 

Lift  Thickness,  mm 


j  L.1 1  L  i  moiuici 

Mean 

►  O/  «■«■■■ 

Standard 

Deviation 

Coefficient  of 
Variation.  Percent 

Comment 

Reference 

134 

8.4 

6.3 

Subbase1 

Mitchell  et  al.  (1977) 

307 

35 

1 1 

Auff  and 

343 

54 

16 

Choummanivong 

(1994) 

93 

17 

18 

Base1 

Auff  and  Laksmanto 
(1994) 

97 

5.4 

5.6 

Auff  and  Yeo  (1992) 

98 

13 

14 

Auff  and  Laksmanto 
(1994) 

103 

6.2 

6.0 

Auff  and  Yeo  (1992) 

104 

23 

22 

Auff  and  Laksmanto 
(1994) 

114 

11 

9.6 

Mitchell  et  al.  (1977) 

131 

11 

8.3 

131 

18 

14 

Auff  and 

133 

19 

14 

(1994) 

145 

29 

20 

Auff  and  Laksmanto 

147 

27 

18 

(1  yyo) 

1  Both  unstabilized  and  stabilized. _  _ — 
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Table  DIO 

Measures  of  Stiffness  and  Strength 


|  Measure 

Mean 

—  ^ 

Standard 

Deviation 

Coefficient  of 

Variation, 

Percent 

Reference 

Angle  of  internal 
friction,  <j) 

N/A 

N/A 

5  to  12 

Schultze 

(1972)' 

Unconfined 
compressive 
strength,  kPa 

N/A 

N/A 

1 1 

(3  to  27) 

Liu  and 

Thompson 

(1966) 

California  bearing 
ratio,  percent 

26 

8.3 

32 

Yoder  and 

Witczak 

(1975)2 

59 

13 

22 

Mitchell  et  al. 
(1977)2 

94 

36 

38 

Yoder  and 

Witczak 

(1975)2 

N/A 

N/A 

15 

(7  to  26) 

Liu  and 

Thompson 

(1966)3 

Benkelman  beam 

0.9 

0.2 

24 

Yeo  and  Auff 
(1995) 

deflections,  mm 

1.0 

0.2 

18 

1.1 

0.3 

24 

1.3 

0.3 

22 

Normality  was  rejected  by  the  chi  square  "goodness-of-fit"  test,  howevc 
was  not  rejected  for  cot<}>  . 

2  Performed  in  the  field. 

3  Performed  in  the  laboratory. 

N/A  Data  not  reported. 

“r,  normality  | 
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Table  Dll 

Measures  of  Stiffness  and  Strength  for  Lime-  and  Cement- 
1  Stabilized  Subbase  and  Base  Course  Materials 

Measure 

Mean 

Standard 

Deviation 

Coefficient  of 

Variation, 

Percent 

Reference 

Unconfined 
compressive 
strength,  kPa 

5,792 

841 

15 

Highway  Research  Board 
(1962a) 

N/A 

N/A 

15 

14  to  24) 

Lui  and  Thompson 
(1966) 

Compression 
modulus,  MPa 

4,000 

2,400 

60 

Marshall  and  Kennedy 
(1974) 

Indirect  tensile 
strenqth,  kPa 

938 

338 

36 

Kennedy,  Hudson,  and 
McCullough  (1975)1 

Tensile 

modulus,  MPa 

621 

421 

68 

California 
bearing  ratio, 
percent 

216 

81 

38 

Mitchell  et  al.  (1 977)2 

307 

95 

31 

N/A 

N/A 

15 

(3  to  31) 

Lui  and  Thompson 
(1966)3 

Plate-load 
tests,  MPa/mm 

N/A 

N/A 

71 

Kennedy,  Hudson,  and 
McCullough  (1 975) 

Dynaflect  tests, 
MPa 

N/A 

N/A 

22 

1  Frequency  distributions  appeared  to  be  normal  in  shape. 

2  Performed  in  the  field. 

3  Performed  in  the  laboratory. 

N/A  Data  not  reported-^^^^  - 
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Table  D12 

Measures  of  Stiffness  and  Strength  for  Asphalt-Stabilized  Base 
Course  Materials 

Measure 

Mean 

Standard 

Deviation 

Coefficient  of 

Variation, 

Percent 

Reference 

Indirect  tensile 
strength,  kPa 

724 

152 

21 

Marshall  and 

Tensile  modulus, 
MPa 

401 

144 

36 

Indirect  poisson's 
ratio 

0.27 

0.13 

48 

Flexural  stiffness, 
MPa 

821 

131 

1  61 

Finn  (1967) 

1,007 

255 

252 

1,014 

283 

282 

Monismith  et  al. 

(1967) 

3,034 

607 

202 

3,613 

1,193 

332 

Finn  (1967) 

3,661 

1,027 

28' 

Fatigue  life  of 
beams 

N/A 

N/A 

30  to  753 

Moore  and 

Kennedy  (1971) 

1  Laboratory-molded  specimens. 

2  Specimens  obtained  from  the  field. 

Standard  deviation  increased  linearly  with  increasing  fatigue  life 

(Ke^nnedv^'Hudson^arKf'Mcc'uUougih'l's^s I**118'**  “  b»  "»P* 

N/A  Data  not  reported. 
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El 


Table  El 

ProDerties  of  Extracted  Asphalt  Cement 

Property 

Mean 

Standard 

Deviation 

Coefficient  of 
Variation,  Percent 

Reference 

Penetration  at 

25°C,  0.1  mm 

21 

2.3 

11 

Sherman  (1971) 

32 

9.6 

30 

35 

4.3 

12 

47 

18 

38 

48 

14 

30 

Viscosity  at  60°C, 
Pa-s 

910 

560 

62 

Sherman  (1971) 

940 

240 

26 

1,160 

1,120 

97 

1,900 

1,270 

67 

4,680 

790 

17 

Table  E2 

1  Mixture  Temperature  in 

the  Field 

Mean 

Standard 

Deviation 

Coefficient  of 
Variation,  Percent 

Reference 

135 

7.4 

5.5 

(2.2  to  8.0) 

Kilpatrick  and  McQuate  (1967) 

153 

9.0 

5.9 

Oqlio  and  Zenewitz  (1965)1 

1  Sampling  error  and  experimental  error  (testing  and  inherent)  accounted  for  2  percent  of 
the  total  variance.  _  _ 
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Table  E4 

Particle  Size  Distribution  for  Extracted  Aggregates,  Percent  Finer 

Sieve  Size 

Mean 

Standard 

Deviation 

Coefficient  of 
Variation,  Percent 

Reference 

25  mm 

94 

4.5 

4.8 

SCSHD  (1966) 

19  mm 

93 

1.4 

1.5 

Granley  (1969) 

92 

3.2 

3.4 

SCSHD  (1966) 

P" 

99 

0.9 

0.9 

Mitchell  et  al.  (1977) 

9.5  mm 

98 

1.7 

1.8 

Oqlio  and  Zenewitz  (1965) 

86 

2.5 

2.9 

Granley  (1969) 

=====  1 

(Sheet  1  of  4) 


Note:  Sampling  error  and  experimental  error  (testing  and  inherent)  accounted  for  the 
following  percentages  of  total  variance. 

-  25-mm  sieve:  85  percent  (SCSHD  1966). 

-  19-mm  sieve:  75  percent  (Granley  1969). 

-  13-mm  sieve:  65  percent  (SCSHD  1966). 

-  9.5-mm  sieve:  60  percent  (Granley  1969). 


E4 


Appendix  E  Asphalt  Concrete 


Table  E4  (Continued)  _ _ _ ====== 

Sieve  Size 

Mean 

Standard 

Deviation 

Coefficient  of 
Variation,  Percent 

Reference 

36 

2.9 

8.0 

Mitchell  et  al.  (1977) 

42 

2.8 

6.7 

Granley  (1969) 

2  0  mm 

35 

3.7 

1 1 

SCSHD  (1966) 

39 

2.1 

5.5 

Oqlio  and  Zenewitz  (1965) 

55 

4.0 

7.4 

SCSHD  (1966) 

0.85  mm 

26 

2.1 

7.8 

Oalio  and  Zenewitz  (1965) 

0.6  mm 

26 

1.7 

6.7 

Granley  (1 969) 

0  425 

18 

1.7 

9.2 

Oqlio  and  Zenewitz  (1965) 

mm 

20 

3.6 

18 

SCSHD  (1966) 

29 

2.1 

7.2 

15 

1.6 

11 

Mitchell  et  al.  (1977) 

15 

1.4 

9.1 

Granley  (1969) 

0.1 8  mm 

9.5 

1.5 

15 

Oqlio  and  Zenewitz  (1965) 

0.15  mm 

9.0 

1.0 

11 

Granley  (1969) 

- - - 

(Sheet  3  of  4) 


Note:  Sampling  error  and  experimental  error  (testing  and  inherent)  accounted  for  the 
following  percentages  of  total  variance. 

-  2.36-mm  sieve:  25  percent  (Granley  1969). 

-  2.0-mm  sieve:  50  percent  (SCSHD  1966). 

-  0.6-mm  sieve:  30  percent  (Granley  1969). 

-  0.425-mm  sieve:  50  percent  (SCSHD  1966). 

-  0.30-mm  sieve:  35  percent  (Granley  1969). 

-  0.1  5-mm  sieve:  30  percent  (Granley  1  969). 

-  0.075-mm  sieve:  50  percent  (SCSHD  1 966). _ _ — _____ 
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Table  E4  (Concluded) 


Sieve  Size 

Mean 

Standard 

Deviation 

Coefficient  of 
Variation,  Percent 

Reference 

No.  200 

2.6 

0.5 

19 

Hode-Keyser  and  Wade  (1963) 

4.1 

0.6 

13 

Shook  (1966) 

4.3 

0.5 

11 

Highway  Research  Board  (1962a) 

4.3 

0.7 

16 

SCSHD  (1966) 

4.9 

1.0 

20 

Shook  (1966) 

5.6 

1.2 

21 

Oglio  and  Zenewitz  (1965) 

5.6 

1.2 

21 

Shook  (1966) 

5.8 

1.1 

20 

5.9 

1.2 

20 

Highway  Research  Board  (1962a) 

6.0 

0.9 

16 

Granley  (1969) 

1.0 

17 

Shook  (1966) 

ism 

1.0 

16 

SCSHD  (1966) 

7.0 

1.1 

16 

Adams  and  Shah  (1965) 

7.1 

0.8 

11 

Mitchell  et  al.  (1977) 

7.1 

1.0 

14 

Adams  and  Shah  (1965) 

7.3 

1.0 

13 

Shook  (1966) 

7.5 

1.3 

17 

Adams  and  Shah  (1965) 

7.9 

0.7 

9.3 

Shook  (1966) 

8.0 

1.3 

16 

Adams  and  Shah  (1965) 

9.6 

0.9 

9.0 

Shook  (1966) 

10 

1.0 

9.5 

10 

0.8 

8.1 

Hode-Keyser  and  Wade  (1963) 

(Sheet  4  of  4)  \ 


aN°'l°n’he  ^era9®'  Sami>li"9  error  and  experimental  error  (testing  and  inherent) 
accounted_for_5Q  percent  of  thejtotalvariancelRranlev  1969;  SCSHD  1966). 
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Table  E5 

Laboratory  Density  and  Voids  Analyses _ _ _  -J 

— 

Property 

Mean 

Standard 

Deviation 

Coefficient  of 
Variation, 

Percent 

Reference 

Bulk  specific 

2.353 

0.038 

1.6 

Highway  Research 

Board  (1962a) 

gravity 

2.388 

0.026 

1.1 

Theo.  Max. 

specific 

qravity 

2.460 

0.011 

0.4 

Granley  (1 969)1 

Voids  total 
mix,  percent 

1.8 

0.8 

44 

Sherman  (1971) 

2.3 

1.5 

65 

2.4 

0.4 

17 

2.5 

1.5 

60 

4.3 

1.0 

23 

Granley  (1 969)2 

4.9 

0.8 

17 

Mitchell  et  al.  (1977) 

5.7 

0.9 

16 

Sherman  (1971) 

6.5 

1.5 

23 

Highway  Research 

Board  (1962a) 

7.7 

1.0 

13 

Voids  filled, 
percent 

56.5 

3.5 

6.2 

Highway  Research 

65.4 

5.3 

8.1 

Board  (1962a) 

70.0 

3.7 

5.3 

Sherman  (1971) 

85.6 

2.2 

2.6 

86.8 

7.5 

8.6 

87.4 

6.8 

7.8 

89.4 

4.3 

4.8 

1  Sampling  error  and  experimental  error  (testing  and  inherent)  accounted  for  30  percent 

of  the  total  variance.  _ . 

1  Sampling  error  and  experimental  error  (testing  and  inherent)  accounted  for  20  percent 

of  the  total  variance. _ _  -1 

E8 


Appendix  E  Asphalt  Concrete 


Table  E6 

Field  Compaction,  Percent  Relative  to 


Laboratory  Density 


Standard 

Deviation 


Coefficient  of 
Variation.  Percent 


Reference 


1  Sampling  error  and  experimental  error  (testing  and 
of  the  total  variance. 

2  Sampling  error  and  experimental  error  (testing  and 
of  the  total  variance. 

3  Airfield. 

4  Highway. 

N/A  Data  not  reported 


Highway  Research  Board  (1962a) 


Kennedy,  Hudson,  and  McCullouqh 
(1975) 


inherent)  accounted  for  35  percent 
inherent)  accounted  for  40  percent 


Table  E7  ! 

Measures  of  Stiffness  and  Strength  — 1 

Measure 

Mean 

Standard 

Deviation 

Coefficient  of 
Variation,  Percent 

Reference 

Marshall  stability, 
kN 

9.9 

1.3 

13 

Mitchell  et  al. 

(1977) 

10.3 

1.3 

12 

Granley  (1969)1 

Marshall  flow,  mm 

2.2 

0.3 

15 

Granley  (1969)* 

2.3 

0.3 

15 

Mitchell  et  al. 

(1977) 

Indirect  tensile 
strength,  kPa 

531 

85 

16 

Marshall  and 

Kennedy  (1974) 

Indirect 

static  modulus, 

MPa 

289 

84 

29 

Marshall  and 

Kennedy  (1 974) 

397 

258 

65 

Kennedy,  Hudson, 
and  McCullough 
(1975) 

654 

477 

73 

Indirect  poisson’s 
ratio 

0.40 

0.11 

27 

Marshall  and 

Kennedy  (1 974) 

N/A 

N/A 

52 

(38  to  73) 

Kennedy,  Hudson, 
and  McCullough 
(1975) 

Dynamic  modulus 
at  4°C,  MPa 

N/A 

N/A 

13 

(9.2  to  16) 

The  Asphalt 

Institute  (1974) 

Dynamic  modulus 
at  21  °C,  MPa 

N/A 

N/A 

nn 

The  Asphalt 

Institute  (1974) 

N/A 

N/A 

40 

(24  to  62) 

Kennedy,  Hudson, 
and  McCullough 
(1975) 

Dynamic  modulus 
at  38°C,  MPa 

N/A 

N/A 

The  Asphalt 

Institute  (1974) 

1  Sampling  error  and  experimental  error  (testing  and  inherent)  accounted  for  60  percent  of 

the  total  variance.  .  _c  _  . 

2  Sampling  error  and  experimental  error  (testing  and  inherent)  accounted  for  75  percent  ot 

the  total  variance. 

N/A  Data  not  reported.  . .  ■ 
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Table  E9 

Mat  Thickness 

Mean,  mm 

Standard 

Deviation 

Coefficient  of 
Variation,  Percent 

Reference 

66 

5.6 

8.4 

Attoh-Okine  and  Roddis  (1994) 

71 

6.6 

9.3 

79 

12 

15 

Kevser  and  Waell  (1968)1 

112 

5.6 

5.0 

Attoh-Okine  and  Roddis  (1994) 

N/A 

6.6 

N/A 

Granley  (1969) 

N/A 

10 

N/A 

Huculak  (1  968) 

N/A 

10 

N/A 

Sherman  (1971) 

N/A 

N/A 

2.7  to  5.9 

Voder  and  Witczak  (1975)2 

N/A 

N/A 

3.5  to  19.2 

Yoder  and  Witczak  (1975)3 

N/A 

N/A 

22 

Voder  and  Witczak  (1975)4 

'  Frequency  distributions  appeared  to  be  normal  in  shape. 

2  New  airfield. 

3  New  highway. 

4  Overlay. 

N/A  Data  not  reported. _ _  _ 

Table  E10 

Benkelman  Beam  Deflections,  mm 

Mean 

Standard 

Deviation 

Coefficient  of 
Variation,  Percent 

Reference 

2.6 

1.8 

69 

Kerr  and  Parkes  (1966)1 

3.2 

1.3 

41 

3.6 

1.7 

47 

N/A 

N/A 

15 

C.G.R.A.  (1962)2 

N/A 

N/A 

22 

C.G.R.A.  (1962)3 

1  Frequency 
3  WASHO  R 
3  State  high 
N/A  Data  no 

distributions  appeared  to  have  a  positive  skew. 

Dad  Test, 
ways. 

t  reported. _ _ _ _ ^=====j 
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Table  Ell 
Falling  Weight 


Measure 


Deflectometer  Results 


Deflections  under 
load,  fjm 


Deflections  at 
0.9  m  offset,  //m 


Backcalculated 
moduli  for  asphalt 
concrete,  MPa 


Backcalculated 
moduli  for  base 
course,  Mpa 


Coefficient  of 
Standard  Variation, 

Deviation  Percent 


Reference 


Siddharthan,  Sebaaly, 
and  Javaregowda 
(1992)' 


Grogan  (1991) 


Siddharthan,  Sebaaly, 
and  Javaregowda 
(1992)1 


Grogan  (1991) 


Siddharthan,  Sebaaly, 
and  Javaregowda 
(1992)' 


Grogan  (1991) 


Siddharthan,  Sebaaly, 
and  Javaregowda 
(1992)2 


Grogan  (1991) 


Siddharthan,  Sebaaly, 
and  Javaregowda 
(1992)2 


Grogan  (1991) 


/Continued) 


’  Frequancy  distributions  appeared  to  be  normal  in  shape. 
_ EH£5H^P^y_ffistributions  appeared  to  have  a  positive  skew. 
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1  Table  Ell  (Concluded)  _ - 

Measure 

Mean 

Standard 

Deviation 

Coefficient  of 

Variation, 

Percent 

Reference 

Backcalculated 

505 

123 

24 

Siddharthan,  Sebaaly, 
and  Javaregowda 
(1992)2 

moduli  for  base 
course,  Mpa 

539 

171 

32 

Backcalculated 

37 

2.3 

6.3 

Siddharthan,  Sebaaly, 
and  Javaregowda 
(1992)2 

moduli  for 
subgrade,  MPa 

38 

2.1 

5.3 

110 

12 

11 

119 

17 

14 

Grogan  (1  991) 

142 

16 

11 

146 

14 

10 

Siddharthan,  Sebaaly, 
and  Javaregowda 
(1992)2 

435 

35 

8.0 

Grogan  (1991) 

483 

116 

24 

2  Frequency  distributions  appeared  tojTaye_a_positivejy<ew:M^_>_M^_<a_g 

Table  El 2 

Initial  Serviceability  Index 

1  Mean 

Standard 

Deviation 

Coefficient  of 
Variation,  Percent 

Reference 

I  4,2 

0.33 

7.9 

Darter,  Hudson,  and  Brown  (1973)' 

I  1  Normality 

was  not  rejected  by  thechi-souare^goodness-of-fit"  test. _  1 
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FI 


Table  FI 

Air  Content  (Percent)  of  Concrete  in  its  Plastic  State 

Mean 

Standard 

Deviation 

Coefficient  of 
Variation,  Percent 

Reference 

4.4 

0.7 

16 

Hanna,  McLauqhlin,  and  Lott  (1967) 

4.6 

0.8 

17 

Neamen  and  Laquros  (1967) 

4.6 

0.8 

18 

Newlon  (1966) 

4.8 

0.7 

14 

4.8 

0.6 

13 

Brown  (1975) 

5.0 

0.8 

16 

Baker  and  McMahon  (1969) 

5.4 

0.8 

15 

Willenbrock  (1974b)1 

5.4 

1.0 

19 

SRCWV  (1968) 

5.5 

0.8 

15 

Baker  and  McMahon  (1969) 

6.5 

1.4 

22 

1  Frequency  distribution  appeared  to  have  negative  skew. 

Note:  On  the  average,  sampling  error  and  experimental  error  (testing  and  inherent) 
accounted  for  25  percent  of  the  total  variance  (Baker  and  McMahon  1 969;  Hanna, 
McLauahlin.  and  Lott  1967;  SRCWV  1968;  Neamen  and  Laquros  1967). 

Table  F2 

Slump  (mm)  of  Concrete  in  its  Plastic  State 

Mean 

Standard 

Deviation 

Coefficient  of 
Variation,  Percent 

Reference 

38 

8 

20 

Willenbrock  (1979b)1 

38 

20 

53 

Neamen  and  Laquros  (1967) 

51 

13 

25 

Kennedy,  Hudson,  and  McCullough  (1975) 

51 

15 

30 

Brown  (1975) 

53 

15 

28 

Baker  and  McMahon  (1969) 

56 

18 

32 

56 

18 

31 

Newlon  (1966) 

61 

20 

33 

SRCWV  (1968) 

76 

25 

33 

Newlon  (1966) 

76 

25 

33 

Hanna,  McLauqhlin,  and  Lott  (1967)  H 

1  Freque 
Note:  Or 
accounte 
SRCWV 

ncy  distribution  appeared  to  be  normal  in  shape. 

the  average,  sampling  error  and  experimental  error  (testing  and  inherent) 
d  for  20  percent  of  the  total  variance  (Hanna,  McLaughlin,  and  Lott  1967; 

1968;  Neamen  and  Laquros  1967;  Baker  and  McMahon  1969). - 
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Table  F3 

Particle  Size  Distribution  for  Concrete  Aggregates,  Percent  Finer 


Sieve  Size 
25  mm 


Mean 

95.6 


Standard 

Deviation 


Coefficient  of 

Variation, 

Percent 


Reference _ 

Louisiana  Department  of 
Highways  (1966) _ 

Willenbrock  (1974b)' _ 

Baker  and  McMahon  (1969) 


5-4  10-8  14  Louisiana  Department  of 

- - - - - Highways  (1966)  _ 

- LL8 - M _ 3J* _ Baker  and  McMahon  (1969) 

^  mm  35.5  12.7  36  Louisiana  Department  of 

- - - -  Highways  (1966) _ 

- - - - LZ _ 1Z _ _  Willenbrock  (1974b)1 _ 

4.75  mm  1.3  1.2  92  Louisiana  Department  of 

- - - Highways  (1966) _ 

— — - U _ 46 _ Willenbrock  (1974b)1 _ 

- ^2^ - LI - _ ®4_ _ Willenbrock  (1974b)1 _ 

4.75  mm  - 96^2 -  1-1  1-1 _ Willenbrock  (1974b)' _ 

97-8  15  1  -5  Louisiana  Department  of 

- ■ - Highways  (1966) _ 

- — mm - Z£^ - ZJj _ 3J _ Willenbrock  (1974b)1 _ 

- -  8  mm - 64^ - 3J3 _ 6.0 _ Willenbrock  (1974b)' _ 

- mm - HLJ - 1^2 _ _ Willenbrock  (1974b)1  _ 

0.30  mm  15.9  6.5  41  Louisiana  Department  of 

- - - Highways  (1966)  _ 

- — - 22J8 - 2^7 _ 13 _ Willenbrock  (1974b)1 _ 

0.15  mm  2.1  1.3  52  Louisiana  Department  of 

- _ - Highways  (1966)  _ 

- - 07  T8 _ Willenbrock  (1974b)1 _ 

Frequency  distribution  appeared  to  be  normal  in  shape. 

Note:  Sampling  error  and  experimental  error  (testing  and  inherent)  accounted  for  the 
following  percentages  of  total  variance. 

-  25-mm  sieve:  30  percent  (Louisiana  Department  of  Highways  1966) 

"  3  mm  !!fVe:  «  PSrCent  !,Baker  and  McMahon  1969;  LA  Department  of  Highways  1966). 

-13-mm  sieve:  85  percent  (Louisiana  Department  of  Highways  1966). 

-4.75-mm  sieve:  50  percent  (Louisiana  Department  of  Highways  1966) 

Note:  Sampling  error  and  experimental  error  (testing  and  inherent)  accounted  for  the 
following  percentages  of  total  variance  (Louisiana  Department  of  Highways  1 966) 

-  4.75-mm  sieve:  15  percent. 

-  1 .1 8-mm  sieve:  8  percent. 

-  0.1 5-mm  sieve:  20  percent. _ 
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Table  F4 

Fineness  Modulus  of  Concrete  Aggregai 

tes 

im 

Standard 

Deviation 

Coefficient  of 
Variation,  Percent 

Reference  1 

2.86 

0.12 

4.2 

Willenbrock  (1 974b)1  1 

1 

H  1  Frequency 

distribution  appeared  to  be  normal  in_sha^g_: _ — ^^==1 

Table  F5 

Density  of  Hardened  Concrete,  kg 

/m3 

Mean 

Standard 

Deviation 

Coefficient  of 
Variation,  Percent 

Reference 

2,371 

14.3 

0.6 

Hanna,  McLauqhlin,  and  Lott  (1967) 

N/A 

N/A 

2.1 

(1.3  to  2.4) 

Kennedy,  Hudson,  and  McCullough 
(1975) 

N/A 

N/A 

1.7 

(1.1  to  4.7) 

Marshall  and  Kennedy  (1974) 

N/A  Datanotreported. _  _ — 
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Table 

I  Comp 

F6 

ressive  Strength  of  Hardened  Concrete,  MPa 

|  Mean 

Standard 

Deviation 

Coefficient  of 
Variation,  Percent 

Reference 

20 

5.5 

28 

Highway  Research  Board  (1962a) 

26 

5.0 

19 

Neamen  and  Laauros  (1967) 

26 

2.0 

7.6 

Highway  Research  Board  (1962a) 

26 

4.1 

16 

SRCWV  (1966) 

26 

2.8 

11 

27 

4.9 

18 

28 

3.9 

14 

(1965) 

28 

4.2 

15 

29 

4.6 

16 

29 

3.8 

13 

Newlon  (1966) 

29 

2.5 

8.6 

Highway  Research  Board  (1962a) 

31 

3.5 

12 

Newlon  (1966) 

31 

4.6 

15 

32 

3.9 

12 

SRCWV  (1966) 

32 

5.5 

17 

33 

3.3 

10 

33 

2.7 

8.0 

Willenbrock  (1974a)' 

33 

4.0 

12 

State  Road  Depaartment  of  Florida 
(1965) 

34 

4.7 

14 

34 

3.7 

11 

Highway  Research  Board  (1962al 

35 

4.2 

12 

State  Road  Depaartment  of  Florida 
(1965) 

35 

4.6 

13 

40 

2.8 

6.9 

Kennedy,  Hudson,  and  McCullouqh 
(1975) 

42 

2.6 

6.1 

Highway  Research  Board  (19fi?at 

Frequency  distribution  appeared  to  be  normal  in  shape. 

Note:  On  the  average,  sampling  error  and  experimental  error  (testing  and  inherent) 

1966!  sRCWV  1966r)Cent  ^  ***  ^  VarianCe  (Neamen  and  Laquros  1967;  Newton 
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Table  F7  j  ^ 

ComDressive  Modulus  and  Poisson's  Ratio  for  Hardened  Concrete  | 

Property 

Mean 

Standard 

Deviation 

Coefficient  of 
Variation,  Percent 

Reference 

Compression 

25,787 

6,189 

24 

Marshall  and 

Kennedy  (1  974) 

N/A 

N/A 

34 

(21  to  49) 

Kennedy,  Hudson, 
and  McCullough 
(1975) 

Poisson's  ratio 

N/A 

N/A 

14 

(9.4  to  20) 

Gibeaut  (1960) 

N/A  Data  not  re 

ported. _ — — =— 

Table  F8 

Measures  of  Tensile  Stiffness  and  Strength 

Measure 

Mean 

Standard 

Deviation 

Coefficient  of 
Variation, 

Percent 

Reference 

Indirect  tensile 
strength,  kPa 

3,248 

648 

20 

Marshall  and  Kennedy 
(1974) 

3,806 

269 

7.0 

Kennedy,  Hudson,  and 
McCullough  (1 975)1 

4,399 

793 

18 

Indirect  tensile 
modulus,  MPa 

23,788 

8,067 

34 

Kennedy,  Hudson,  and 
McCullouqh  (1975)’ 

27,511 

9,377 

34 

Marshall  and  Kennedy 
(1974) 

Flexural 
strength,  kPa 

3,792 

255 

6.7 

Highway  Research 

Board  (1962a) 

4,344 

241 

5.6 

4,461 

359 

8.0 

Kennedy,  Hudson,  and 
McCullouqh  (1975) 

4,564 

434 

9.5 

Brown  (1  975) 

4,895 

365 

7.5 

Highway  Research 

Board  (1962a) 

5,282 

421 

8.0 

Brown  (1975) 

6,068 

365 

6.0 

Highway  Research 

Board  (1962a) 

’  Frequency  distributions  appeared  to  be  normal  in  shape. 

Note:  Frequency  distributions  for  flexural  strength  appeared  normal  in  shape  (Kher  and 
Darter  1973).  _ _ _ ^^==^=====^  1 
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Table  F10 

Load  Transfer  (Percent) 


Joint  Type 

Mean 

Standard 

Deviation 

Coefficient  of 

Variation, 

Percent 

Reference 

Doweled 

expansion 

30.5 

7.4 

24 

Rollings  (1987) 

Doweled 

12.2’ 

3.6 

29 

Hammons,  Pittman,  and 

construction 

15.5' 

5.8 

38 

Mathews  o  yyo) 

16.32 

4.1 

25 

19.3' 

1.8 

9.3 

28. 52 

4.6 

16 

30.6 

1 1.6 

38 

Rollinqs  (1 987) 

Keyed 

construction 

25.4 

10.5 

41 

Rollings  (1 987) 

Doweled 

13.0’ 

6.6 

50 

Hammons,  Pittman,  and 

contraction 

16.3' 

9.1 

56 

Mathews  (1 9yb) 

17.9' 

3.9 

22 

21.2’ 

2.1 

9.9 

26. 62 

3.1 

12 

30.72 

1.1 

3.5 

Tied 

24. 72 

8.6 

35 

Hammons,  Pittman,  and 

contraction 

25.5’ 

oo 

6.9 

Mathews  (1 99o) 

27.0' 

4.3 

16 

Plain 

6.4' 

2.9 

46 

Hammons,  Pittman,  and 

contraction 

9.0' 

7.4 

82 

Mathews  (199o) 

1 0.62 

3.2 

30 

10. 82 

3.2 

30 

11.9' 

2.0 

17 

11.9' 

3.4 

29 

12.82 

1.6 

13 

13.62 

1.6 

12 

14. 22 

3.7 

26 

17.4' 

6.6 

38 

(Continued) 


'  Cold  testing  conditions. 

2  Warm  testing  conditions. 

3  Frequency  distribution  appeared  to  have  positive  skew. 

4  Frequency  distribution  appeared  tc^bejTorrnaM^shape^ 
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Table  F10  (Concluded) 

Joint  Type 

Mean 

Standard 

Deviation 

Coefficient  of 

Variation, 

Percent 

Reference 

Plain 

17.42 

3.4 

19 

Hammons,  Pittman,  and 

Mathews  (1 995)3 

contraction 

(Cont.) 

17.5’ 

4.3 

25 

17.62 

2.2 

13 

17.6’ 

5.3 

30 

20. 02 

3.1 

16 

20.32 

5.5 

27 

20.4' 

6.7 

33 

20.8’ 

2.6 

12 

25. 22 

2.2 

8.7 

37.2 

7.1 

19 

Rollings  (1987) 

'  Cold  testing  conditions. 

Warm  testing  conditions. 

Table  F1 1 

Falling  Weight  Deflectometer  Results 


Measure 

Mean 

Standard 

Deviation 

Coefficient  of 
Variation,  Percent 

_  ..  1 

Reference 

Deflections  under  load. 

81 

8.9 

11 

'  1 

Grogan 

102 

2.9 

2.8 

(1991) 

396 

45 

11 

4,53 

121 

27 

Backcalculated  moduli 

25,690 

2,830 

11 

concrete,  MPa 

29,150 

15,740 

54 

Backcalculated  moduli 

80 

14 

17 

109 

17 

16 

167 

32 

19 

Table  FI  2 

Serviceability  Index  for  Rigid  Pavements 


Mean 


Standard 

Deviation 


Coefficient  of 
Variation,  Percent 


Reference 


4.7 


0.14 


3.0 


Fordyce  and  Teske  (1963) 


J__Frgquency  distribution  appeared  to  be  normal  in  shape. 
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